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2 Be An SE 1. Ziffernmäßig arbeitende Rechenautomaten 

h 6 _ Wenn wir eine Zahl, z. B. 673, in der üblichen Weise schreiben, so heißen die Symbole 6, 7, 
B 3die Ziffern der Zahl, und man nennt eine mechanische oder elektrische Einrichtung, welche 
EN die einzelnen Ziffern jeder Zahl bearbeitet und anzeigt, oft Ziffern-Rechenmaschine. 
e- - — Inden vergangenen 12 Jahren hat eine beachtliche Entwicklung dieser Maschinen stattgefunden, 
so daß sie zwei wichtige Eigenheiten gewonnen haben: Erstens können sie langwierige und ver- 
4 wickelte Rechnungen gänzlich automatisch ausführen, wenn der Rechengang in geeigneter, Form 
er eingegeben worden ist, und zweitens sind sie sehr anpassungsfähig, so daß man mit einundderselben 
Maschine die verschiedensten Arten von Rechnungen durchführen kann, z.B. die Ermittlung von 


Funktionswerten aus der Pötenzreihe, die Lösung umfangreicher linearer Gleichungssysteme, die 
Bestimmung der Eigenwerte von Matrizen, die schrittweise Integration gewöhnlicher Differential- 
gleichungen. Um diese beiden Eigenheiten hervorzuheben, nennt man diese Maschinen manch- 
mal universell verwendbare automatische Ziffern-Rechenmaschinen, kürzer 
Rechenautomaten. Meine Absicht ist es, über einige Hauptzüge solcher Maschinen zu berichten, 
sowohl hinsichtlich der Maschinen selbst, als auch hinsichtlich der Art ihrer Anwendung. Frei- 
lich kann ich manche Dinge nur streifen, denn ein einzelner Vortrag ist viel zu kurz, um den 
Gegenstand erschöpfend zu behandeln. (Für einen umfassenden allgemeinen Überblick vgl. 
- Literaturhinweis [10], für genauere Berichte über einzelne Maschinen vgl. [1, 4, 8, 9, 11, 13, 14 
Bund 19)). 
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2. Die „Analytieal Engine‘ von Babbage 


B Obgleich erst vor etwa 8 Jahren die erste derartige Maschine in Betrieb genommen wurde, 
ist der Grundgedanke nicht neu. Er ist vielmehr schon ungefähr 115 Jahre alt und stammt im 
wesentlichen von Charles Babbage, dessen geplante Maschine, die er analytical engine nannte, 
ein Rechenautomat werden sollte, automatisch bis auf eine Ausnahme, nämlich bis auf den Ge- 
_ brauch von Funktionentafeln, der eine Bedienungsperson erforderte. Die analytical engine sollte 
aus drei Hauptteilen bestehen: in dem einen, den Babbage store (Speicherwerk) nannte, konnten 
_ numerische Daten in einer Reihe von Speichern niedergelegt werden; in einem zweiten, den er 
 mäll (Rechenwerk) nannte, konnten arithmetische Operationen ausgeführt werden an Zahlen, 
die ihm vom Speicherwerk zugeführt wurden; ein dritter nicht besonders benannter Teil wählte 
- die Speicher aus, von denen.Zahlen entnommen werden sollten, gab die auszuführende Operation 
an und bestimmte den Speicher, in welchen das Ergebnis überführt werden sollte. 

Diese Steuerung der Maschine sollte durch eine Folge von Lochkarten geschehen. Fühler, 
die durch die Löcher in den Karten hindurchgriffen, sollten die Speicher für die Übertragungen 
_ zum und vom Rechenwerk und die auszuführende Rechenoperation auswählen. 


3. Allgemeine Gesichtspunkte für Rechenautomaten 


= ER Die Betrachtung von Aufbau und Ablauf einer Rechnung, die von einem Menschen mit 

- Hilfe einer gewöhnlichen Rechenmaschine ausgeführt wird, soll zeigen, welche Einrichtungen 
rgesehen werden müssen, um die gleiche Rechnung automatisch auszuführen. Die Schluß- 
lgerungen lassen sich an Aufbau und Wirkungsweise von Babbages analytical engine ver- 
- Bild 1a zeigt schematisch Aufbau und Ablauf einer ‚„Handrechnung‘. Der Rechner ver- 

über eine gewöhnliche Rechenmaschine für die arithmetischen Operationen, ferner über 
Arbeitsblatt, um Ausgangsdaten, Zwischenergebnisse und die erforderliche Reihenfolge der 
jerationen festzuhalten; er selbst steuert die Übertragung von Zahlen in beiden Rich- 
wischen der Rechenmaschine, dem Arbeitsblatt und den Funktionentafeln. Weiter 
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steuert er das Arbeiten der Rechenmaschine und rain den Arbeitsblatt die Vorschiiien 
für die auszuführenden Übertragungen und Operationen. 

Eine automatische Maschine kann man sich von ganz ähnlicher Arbeitsweise vorstellen 
(vgl. Bild 1b). Sie muß ein Rechenwerk haben (die mill in Babbages analytical engine), 
um die arithmetischen Operationen auszuführen, ein Speicherwerk für Ausgangsdaten, 
Zwischenergebnisse und Rechenvorschriften,,und ein Kommandowerk an Stelle des 
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Bildi1. a) Aufbau und Ablauf einer Handrechnung. b) Aufbau und Arbeitsweise eines Rechenautomaten ( 
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Rechners, der bei einer Handrechnung die Operationsfolge steuert. Weiter benötigt die Maschine 
| eine Eingabe- und eine Auslieferungs-Einrichtung, um von außen Rechenvorschriften und 
r zahlenmäßige Ausgangswerte aufnehmen und Ergebnisse herausgeben zu können. Das Speicher- 
E werk kann so eingerichtet werden, daß es die Aufgaben des Arbeitsblattes und der Funktionen- 
tafeln in der Handrechnung mit übernimmt und es kann zweckmäßig in zwei Teile aufgespalten 
sein, nämlich in ein Hauptspeicherwerk mit direkter Verbindungsmöglichkeit zum Rechenwerk 
und ein Hilfsspeicherwerk, in welches Gruppen von Zahlen und Rechenvorschriften überführt 
werden können, wenn sie im Laufe der Rechnung längere Zeit nicht gebraucht werden. 


4. Die Darstellung von Zahlen in der Maschine 


Oft ist es zweckmäßig, die Ziffern einer Zahl durch technische Elemente darzustellen, die 
nur zweier sich gegenseitig ausschließender Anzeigen fähig sind. Beispielsweise können wir eine 
‚Reihe von Orten auf einer Karte benutzen, an denen jeweils entweder ein Loch gestanzt ist oder 
nicht, oder eine Reihe von Relais, deren jedes entweder ein- oder ausgeschaltet ist, oder eine 
‘Folge von Signalen, deren jedes aus der Anwesenheit oder Abwesenheit eines elektrischen Im- 
pulses besteht. Der Gebrauch solcher Elemente legt die Darstellung der Zahlen in der Maschine 
im Dualsystem nahe, ohne es allerdings zwingend zu erfordern; denn man kann auch Dezimal- 
zahlen in verschiedener Verschlüsselung mit Hilfe solcher Elemente darstellen. Verschiedene 
Maschinen arbeiten innerhalb der Maschine mit Dualzahlen, fast alle nehmen jedoch die ein- 
gehenden Zahlenwerte im Dezimalsystem auf und liefern die Ergebnisse ebenfalls im Dezimal- l 

. system; die Dezimal-Dual- und die Dual-Dezimal-Umwandlung wird durch die Maschine als 
Teil des Eingabe- bzw. Auslieferungsvorganges erledigt. 


2. 
e 

4 
Kr 
ne 
y 


7 
1 
D 
0 
o 


Bild2. Darstellung von Zahlen durch Impulse. a) Paralleldarstellung der Ziffern. b) Seriendarstellung der Zittern. 
c) Parallele Dualdarstellung jeder Dezimalzifter, is or der ee i 


Die Ziffern einer Zahl können dargestellt werden entweder durch die Ka en Zu ist 
einer Anzahl von Elementen, oder durch die zeitlich aufeinanderfolgenden Zustände eines einz 
Elementes. Beispielsweise können die Ziffern einer Dualzahl durch ‚gle ge zii 
dargestellt werden, nämlich durch das Vorhandensein oder Fehlen ein ses. 
Leitung für jede Dualstelle in einem S ystem von Leitungen, die z 
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zigenLeitung (vgl. Bild 2b). Die erste Darstellung heißt Parallel-, die zweite Serien- Darstellung, 
und entsprechend können wir von Parallel- oder Serienbetrieb sprechen, je nachdem die Ziffern 
einer Zahl gleichzeitig oder nacheinander bearbeit werden. Auch eine gemischte Serien-Parallel- 
Darstellung ist möglich. Zum Beispiel kann man in einer Dezimal-Maschine jede einzelne Dezimal- 
ziffer dual im Parallelbetrieb darstellen, die einzelnen Dezimalstellen jedoch im Serien- 
betrieb behandeln (vgl. Bild 2e). 

i Da bei der Addition der Übertrag stets von der niedrigeren in die nächsthöhere Stelle er- 
folgt, ist es bei einer Serien-Maschine am zweckmäßigsten, die Zahlen so darzustellen, daß die 
niedrigste Stelle zeitlich zuerst kommt. 


5. Rechenoperationen 


Bei einer Dezimalmaschine kann man die Addition auf zwei verschiedene Arten ver- 
wirklichen: die eine beruht auf dem Zählen, d.h. auf fortgesetzter Addition der Einheit, die 
andere benutzt eine ‚‚Additionstafel‘, mit der man z. B. bei der Addition von 3 zu 4 direkt zum 
Ergebnis 7 geht, ohne durch 5 und 6 als Zwischenergebnisse zu gehen. Bei einer Dualmaschine 
besteht zwischen beiden Methoden kein Unterschied; beide umfassen nur die 4 Beziehungen 

0+0=0 1+0=1 

0+1=1 1+1=10 | 
(die letzte muß natürlich im Dualsystem gelesen werden, worin 10 ‚eine Zwei + keine Eins“ be- 
deutet). In einer Serien-Dualmaschine wird in keinem Augenblick mehr getan als die Addition 
einer einzelnen Ziffer O oder 1 zu einer anderen einzelnen O0 oder 1; die Tragweite und die Ge- 
schwindigkeit der Maschine beruhen auf der Möglichkeit, diese elementaren Operationen sehr 
schnell — etwa in der Größenordnung von einer Million je Sekunde — und auf eine Weise auszu- 
führen, daß sie bequem angeordnet und gesteuert werden können. 

Die Multiplikation kann entweder durch fortgesetzte Addition oder mit Hilfe einer Multi- 
plikationstafel erfolgen; diese Unterscheidung ist wiederum nur bei Dezimalmaschinen von 
Bedeutung. 

Die Behandlung der Vorzeichen und die Stellung des Dezimal- bzw. Dualkommas kann 
ich hier nur als weitere zu beachtende Gesichtspunkte erwähnen. Es gibt verschiedene Möglich- 
keiten dafür, jedoch würde deren Besprechung zu viel Zeit erfordern (für eine kurze Besprechung 
vgl. [10], Seite 64). 

| 6. Die Form der Rechenvorschriften 


In der analytical engine von Babbage wurden die Rechenvorschriften in ganz anderer Art 
wiedergegeben als die Zahlen und demgemäß auch in einem besonderen Speichersystem nieder- 
gelegt. Die Zahlen nämlich wurden dargestellt durch die Stellungen einer Reihe von Zählrädern, 
die Rechenvorschriften hingegen wurden in Lochkarten niedergelegt. Eineähnliche Unterscheidung 
in der Darstellung und Speicherung von Zahlen und Rechenvorschriften findet man noch bei den 
ersten tatsächlich gebauten Maschinen, von denen zwei im 7. Abschnitt behandelt werden. Aber 
bei neueren Maschinen, darunter wohl allen, die im vergangenen Jahr fertiggestellt wurden oder 
die noch im Bau sind, werden Rechenvorschriften in derselben Art verschlüsselt und gespeichert 
wie Zahlen; ein Unterschied besteht lediglich in der Art ihrer Verwendung. So werden z.B. in 
einer Serienmaschine, mit Darstellung durch elektrische Impulszüge die den Rechenvorschriften 
entsprechenden Züge vom Speicher zum Kommandowerk und umgekehrt übertragen, die den 
Zahlen entsprechenden Züge vom Speicher zum Rechenwerk; diese Übertragungen selbst werden 
wieder vom Kommandowerk gesteuert. 

Wie auch das Speicherwerk physikalisch beschaffen sein mag, es muß aus einer Anzahl 
identischer Einzelspeicher bestehen, die man sich zweckmäßig numeriert denkt. Die Nummer 
jedes Einzelspeichers heißt gewöhnlich seine Adresse oder die Adresse seines Inhalts. So 
bezeichnet man den Inhalt des Speichers Nr. n mit O (n). 

Man unterscheidet zwei Hauptarten von Rechenvorschriften (Befehlen), die an einem 
Beispiel gezeigt werden sollen: Nehmen wir an, es soll die Summe der in den Speichern n, und n, 
enthaltenen Zahlen gebildet und das Ergebnis nach Speicher n, gebracht werden, so kann das 
erstens auf ‚Grund eines einzigen Befehls geschehen, der sich schreiben läßt: 
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Die are Möglichkeit ist folgende: Das Rechenwerk möge einen sogenanntes Sammelwerk (ab- 


gekürzt: Sm.) haben, d.h. ein Gerät, das eingegebene Zahlen solange aufaddiert, bis es gelöscht 


Do wird. ‚Dann kann die gewünschte Rechenoperation auf Grund von drei Befehlen erfolgen: 


C (n,) > Sm 
il 6 Anl. 
C(Sm)—n; 3 
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Jeder Befehl von der Form (1) gibt drei Adressen im Speicher an, jeder von der Form (2) nur 
eine einzige. Daher bezeichnet man ihn auch als Drei-Adressen-Befehl bzw. Ein- 
Adressen-Befehl. 

Außer den einzelnen Befehlen selbst muß man aber auch die Reihenfolge festlegen, in der 
sie ausgeführt werden sollen. Bei einer Maschine, in der die Befehle in derselben Form wie die 
Zahlen dargestellt werden, kann das auf zwei Arten geschehen. Entweder man nimmt in jeden 
Befehl diejenige Adresse auf, von welcher der nächste Befehl entnommen werden soll; dadurch 
wird aus einem Drei-Adressen-Befehl ein Vier-Adressen-Befehl. Oder man speichert für nor- 
malen Rechengang die einzelnen Befehle mit Adressen, die in derselben Reihenfolge numeriert 
sind, in der sie zeitlich nacheinander ausgeführt werden sollen. Die Adresse des gerade ausge- 
geführten Befehls wird in einem Speicher festgehalten, dessen Inhalt normalerweise nach Ab- 
arbeitung dieses Befehls um eins vermehrt wird, und der Inhalt dieses Speichers steuert die 
Auswahl des nächsten Befehls. Will man von der Reihenfolge, in der die Befehle im Speicher- 
werk untergebracht sind, abweichen, so bedarf es eines besonderen Befehls. 

Das Steuerungssystem einer Maschine hängt davon ab, welche allgemeine Form man für 
die Befehle angenommen hat und in welcher Weise man den jeweils nächsten Befehl auswählen 
will. Ist die Maschine einmal gebaut, so müssen die Befehle dem Typus entsprechen, für dessen 
Verarbeitung das Steuerungssystem entworfen worden ist. 


7. Der Harvard Mark I Caleulator und der ENIAC 


Der erste tatsächlich gebaute Rechenautomat war der Harvard Mark ICalcu- 
lator[l, 11], der durch Prof. H.H. Aiken von Harvard und die IBM Company entwickelt 
wurde und jetzt im Computation Laboratory der Harvard-Universität aufgestellt ist. Diese 
Maschine, die als „Verwirklichung von Babbages Traum‘‘ begrüßt wurde, war schon vor dem 
Krieg geplant, wurde aber erst in ihm fertiggestellt. 

Sie ist eine Parallel-Dezimal-Maschine und hat eine Anzahl mechanischer Zählwerke, die 
sowohl als Addiergeräte als auch als Speicher dienen. Sie werden von einer gleichmäßig- um- 


laufenden Welle über elektromagnetische Kupplungen angetrieben, die durch Relais gesteuert 


Bild3. Harvard Mark I Oaleulator. Links: Schalter zum-Eingeben konstanter Werte, Mitte: Rechenwerke. 
Rechts: Interpolierwerke und Kommandogerät 3 


werden. Diese Relais wiederum werden durch Schaltkreise mit Kontakten anderer Relais oder 
Zähler gesteuert oder durch Löcher in Lochkarten oder Lochstreifen. Die Reihenfolge der 


Operationen wird durch einen Lochstreifen gesteuert, auf dem die Befehle in verschlüsselter 


Form in aufeinanderfolgende Zeilen gelocht sind. Stromkreise, die durch die Löcher geschlossen 


werden, betätigen Relais, welche die durch die Löcher kommandierte Operation herbeiführen. ” 
Bild 3 zeigt eine Gesamtansicht der Maschine, Bild 4 in einem Ausschnitt ein Steuerungsgerät 


für den Rechenablauf mit eingesetztem Lochstreifen.. & & 
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Die Verwendung elektronischer Rechen- und Schaltkreise erhöht die Arbeitsgeschwindigkeit 
einer Maschine beträchtlich. Die erste Maschine, die elektronische Elemente benutzte, war der 
ENIAC[S, 9] von Dr. J. W. Mauchly und Dr. J. P. Eckert. In dieser Maschine werden elektrische 


ru 


Bild4. Kommandogerät des Harvard MarkI Calculator 


Impulse mit einer Frequenz von 200 kHz durch elektronische Schaltkreise zu elektronischen 
Rechenkreisen geleitet. Sie ist eine Parallel-Dezimal-Maschine, und besteht aus einer Reihe von 
Elementen, die durch drei Gruppen von elektrischen Übertragungsleitungen miteinander ver- 


Ä 
Bild 5. Gesamtansicht des ENIAC 


bunden sind. Die erste Leitungsgruppe versorgt alle Elemente mit einer Grundauswahl von 
“ Impulsen, die zweite („Ziffern-Leitungen‘) überträgt Impulse zur Darstellung von Zahlen 
zwischen den Elementen, die dritte (‚‚Programm-Leitungen“) überträgt die Steuerimpulse, welche 
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die Recheneinheiten in Betrieb setzen und so den Ablauf der von der Maschine ausgeführten 
Operationen steuern. 

Ausgangswerte können in die Maschine eingegeben werden an Tafeln mit Handschaltern 
(für Tabellenwerte benutzt) oder durch Lochkarten, deren Inhalt in Relais-Speicher überführt 
wird. Zwischenergebnisse werden in elektronischen Zählwerken, die gleichzeitig als Addier- 
elemente dienen, gespeichert. 

Die Speicherung der Befehle besteht einerseits in den Verbindungen zwischen den einzelnen 
Rechenelementen und den Ziffern- und Programmleitungen, andererseits in der Einstellung von 
verschiedenen Steuerungsschaltern an den einzelnen Elementen, was alles vor Inbetriebnahme 
der Maschine von Hand erledigt werden muß. 

Bild 5 gibt eine Gesamtansicht des ENIAC, Bild 6 zeigt im Ausschnitt zwei Rechenelemente 
mit ihren Verbindungen zu den Ziffernleitungen (in der Mitte) und zu den- Programmleitungen 
(unten). 

Über beide Maschinen sind Berichte veröffentlicht worden, auf die für Einzelheiten ver-. 
wıesen sel. 


Bild 6. ENIAC. Zwei Addierwerke und Verbindungsleitungen 


8. Neuere Entwicklungen 


Die neueren Maschinen unterscheiden sich von den soeben kurz erwähnten in zwei wesent- 
lichen Eigenheiten: 

Erstens haben sie ein weit größeres Speichervermögen, während sie gleichzeitig mechanisch 
und elektrisch einfacher geworden sind. Dies beruht auf der Entwicklung verschiedener. Arten 
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von Speichern, die große Speicherkapazität ohne großen technischen Aufwand geben. Drei 
Arten von Speichern sind für diesen Zweck entwickelt worden: 

i 

| 


1. akustische Verzögerungsstrecken, in denen Informationen (Zahlen oder Rechenbefehle) 

als Züge akustischer Impulse in einer Quecksilbersäule gespeichert werden; diese liegt zwischen 
| zwei piezoelektrischen Kristallen, die elektrische Impulse in akustische umwandeln und um- 
f gekehrt [19], 
: 2. magnetisches Material, in welchem die Informationen durch den magnetischen Zustand 
.. des Materials festgehalten werden, 
; 3. ein isolierender Schirm, auf dem die Informationen als elektrische Ladungsverteilung 
| gespeichert werden. 


Die erste Art wird in mehreren Maschinen benutzt, die zweite als Hauptspeicher im Har- 
vard Mark III, einer kürzlich fertiggestellten Maschine und als Hilfsspeicher in mehreren 
anderen Maschinen. Ein Speicher der dritten Art wurde von Professor Williams an der Uni- 
versität Manchester [21] entwickelt, wo eine Maschine mit einer solchen Speicherung gebaut 
worden ist [14]. 

Die zweite Eigenheit der neuesten Maschinen wurde schon in Abschnitt 6 erwähnt, nämlich 
die Speicherung von Befehlen und Zahlen in derselben Form. Dies hat die wichtige Folge, daß 
Befehle genau so wie Zahlen behandelt und übertragen und mit Hilfe des Rechenwerks aufgebaut 
oder abgeändert werden können. Wenn die Einzelspeicher durch Nummern unterschieden sind, 
so wird ein Befehl gewöhnlich wenigstens eine Nummer als Adresse enthalten, an die er gerichtet 

ist. Es ist also eine ganz naturgemäße Anwendung des Rechenwerkes, einen Befehl herzunehmen 
und beispielsweise eine Eins zu der darin enthaltenen Adresse zu addieren. Diese Freiheit, Be- 
fehle im Laufe der Rechnung abzuändern, sogar, wenn nötig, in einer Weise, die von Ergebnissen 
der Rechnung selbst abhängt, gestattet große Freizügigkeit beim Entwerfen des Rechenprogramms, 
das sogar auf Rechnungen zugeschnitten werden kann, die eine Menge von Entscheidungen und 
die Auswahl unter verschiedenen Vorgehensweisen erfordern. 


9. Der EDSAC 


s Statt den Versuch zu machen, die verschiedenen Eigenheiten einer Anzahl von Maschinen 
aufzuzählen, was bedeuten würde, daß ich über jede einzelne nur wenig sagen könnte, will ich 
mich darauf beschränken, über die Maschine zu berichten, mit derich am meisten vertraut bin 
nämlich den EDSAC (Electronic Delay-Storage Automatic Calculator) im Mathematical Labora- 
tory der Universität Cambridge [18, 19]. 


y Bild 7, EDSAC (im Mathematical Laboratory Cambridge/England). Rechts: Die Lochstreifenabtaster des 
- Eingabewerks und Fernschreibmaschinen des Auslieferungswerks 


Er ist eine Serien-Dual-Maschine mit Ultraschallspeichern. Er benutzt Ein-Adressen- 
Befehle und Serienspeicherung für sie. Eine Ansicht gibt Bild 7, der allgemeine Aufbau geht in 
etwas vereinfachter Form aus Bild8 hervor. Die Kreuze bezeichnen Schalter, die vom Kommando- 
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werk bedient werden und die Übertragung von Impulszügen zwischen den verschiedenen Ele- 


menten der Maschine steuern. 


Es gibt eine Speichereinheit zum Aufnehmen des laufenden Befehls (,,Speicher für laufen- 
den Befehl‘ in Bild 8), und eine weitere zum Aufnehmen der Adresse, von der der laufende Be- 
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Bild 8. Allgemeiner Aufbau des EDSAC. Vom Kommandowerk 
betätigte Schaltersind mit x bezeichnet 


zum ande . 


zum Kommando- 
werk 


fehl stammt (,,Speicher für Rechenablauf- 
Steuerung‘). Der Inhalt der zweiten Spei- 
chereinheit wird normalerweise nach Ausfüh- 
rung jedes Befehls um eins vermehrt, aber es 
ist auch möglich, in sie die in dem laufenden 
Befehl enthaltene Adresse selbst zu übertra- 
gen und so von der Reihenfolge, in der die 
Befehle gespeichert sind, abzuweichen. 


Die Maschine führt abwechselnd zwei 
Arbeitsgänge aus: Im ersten wird ein Be- 
fehl ausgewählt entsprechend der Nummer, 
die im Speicher für Rechenablaufsteuerung 
steht und in den Speicher für laufenden Be- 
fehl gebracht. Im zweiten Arbeitsgang 
überträgt sie diesen Befehl in das Kom- 
mandowerk auf einen statischen Speicher, 
von dem aus dann die einzelnen Schaltele- 
mente betätigt werden, die den Ablauf der 
geforderten Operation bestimmen. 


Zahlen und Befehle werden in die Ma- 
schine in verschlüsselter Form auf Loch- 
streifen einer Fernschreibmaschine eingege- 
ben und die Ergebnisse von einer etwas ab- 
geänderten Fernschreibmaschine ausgelie- 
fert. Zahlen werden als Dezimalzahlen ein- 


gegeben und ausgeliefert; die Dezimal-Dual- 
Umwandlung bei der Eingabe und die Dual-Dezimal-Umwandlung bei der Auslieferung wer- 
den von der Maschine nach den entsprechenden Befehlen ausgeführt. 


10. Programmgestaltung und Verschlüsselung | 


Dieser und der folgende Abschnitt beschäftigen sich mit der Gestaltung der Rechnung bei 
einem Rechenautomaten. Diese Arbeit hat zwei Stufen: Die erste zuweilen als Programm- 
gestaltung bezeichnete Stufe besteht darin, die ganze Rechnung in eine Folge elementarer von 
der Maschine ausführbarer Operationen aufzuspalten wie Addition, Multiplikation, Übertragung 
und Auswahl des nächsten Befehls. Ein ähnliches Vorgehen ist bei einer Handrechnung er- 
forderlich; bevor man mit irgendeiner Rechnung anfangen kann, muß man sich darüber klar 
werden, wie man vorgehen will. Im Prinzip ist diese Arbeitsstufe weitgehend unabhängig von der 


Maschine, die wir für die Rechnung gebrauchen wollen. In Einzelheiten hängt sie freilich von 


der Maschine ab, da nicht alle Maschinen die gleichen elementaren Operationen zur Verfügung 


haben. Bei manchen kann z. B. |x| aus x durch eine einzige Operation gefunden werden, während 


man bei anderen ein Programm dafür aufstellen muß. 


Die zweite, als Verschlüsselung bezeichnete Arbeitsstufe besteht darin, das Schema der 
Befehle in die Form zu bringen, in welcher sie in die Maschine eingegeben werden sollen. Dies 
hängt von der jeweiligen Maschine ab; weniger von der Art der Speicher oder der Art, in der die 
Zahlen innerhalb der Maschine wiedergegeben und verarbeitet werden, vielmehr hauptsächlich 
von den beiden erwähnten Eigenheiten, nämlich der für die Befehle gewählten Grundform und 
der Art, wie man den jeweils nachfolgenden Befehl auswählen läßt. i 


11. Beispiel für die Verschlüsselung ” 


{ 


Die Unterschiede in der Form der Befehle und in der Art, den nächsten Befehl auszuwählen, 
erschweren es, auf kurzem Raume die verschiedenen Möglichkeiten allgemein und anschaulich ’ 
darzustellen. Ich will mich deshalb wieder auf die Maschine beschränken, die mir am ver- 


trautesten ist. Ihre für unsere Betrachtung wichtigsten Eigenheiten sind die Verwendung d 
Ein-Adressen-Form für Befehle mit Serienspeicherung und die Tatsache, daß ihre Spei 
Inhalt eines Einzelspeichers so lange unverändert enthalten, bis man eine andere Zahl 
anderen Befehl in ihn einführt. Im folgenden sind einige Befehle in abgekürzter Form zu 
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gestellt, die für ein Beispiel, das ich im Anschluß besprechen will, benötigt werden (wegen einer 
vollständigen Liste der elementaren Operationen und der Form der zugehörigen Befehle vgl. [17]): 


An: Addiere C(n) zum Inhalt des Sammelwerks. 
Sn: Addiere —C(n) zum Inhalt des Sammelwerks. 


T n: Übertrage den Inhalt des Sammelwerks in den Speicher Nr. n und lösche das Sammel- 
werk. 


N 
Un: Übertrage den Inhalt des Sammelwerks in den Speicher Nr.n, lösche aber das Sammel- 
: werk nicht 


En: Prüfe den Inhalt des Sammelwerks; wenn er null oder positiv ist, nimm ( (n) als 
nächsten Befehl; ist er negativ, gehe in der normalen Reihenfolge weiter, d. h., wenn 
dieser Befehl En im Speicher m stand, so entnimm den nächsten dem Speicher 


m+1. 


Betrachten wir die folgerde Aufgabe (die man sich als Teil einer größeren Rechnung denken 
mag): Addiere die Inhalte der Speicher 100 bis 140 und setze die Summe in den Speicher 170; 
symbolisch geschrieben: 

2 40 
ee 


# j= 


von Befehlen A 100, A 101, A 102, ..., 4 140, T 170. Das setzt voraus, daß das Sammelwerk 
zu Anfang infolge eines vorausgegangenen Befehls gelöscht ist. Ist dies nicht der Fall, so muß 
es erst durch den Befehl 7’ 0 gelöscht werden; es kann nur durch einen 7-Befehl gelöscht werden, 
und Einzelspeicher Nr. 0 wird vereinbarungsgemäß als ‚„„‚Müllgrube‘‘ benutzt, in welche ein nicht 
mehr gebrauchter Inhalt des Sammelwerks zum Zwecke des Löschens übertragen wird. 


Diese Folge von Befehlen beansprucht außergewöhnlich viel Speicherraum und die Rech- 
nung kann folgendermaßen viel sparsamer ausgeführt werden: Wir benutzen die drei Befehle 
A 170, An, T 170, wodurch jeweils 0 (170) um CO m) vermehrt wird, und zwar immer wieder, 
wobei n der Reihe nach die Werte 100, 101, 120, ..., 140 haben soll. Das bedeutet, daß wir 
den Wert von n, der die Adresse im Befehl A n ist, "jedesmal um eins vermehren müssen, wenn - 
der Befehl ausgeführt ist. 


Das kann geschehen, indem man den Befehl A n SR ins Rechenwerk übernimmt und zu 
der Stelle, die der Einheit in der Adresse n entspricht, jeweils eine Eins addiert. Die Zahl, die 
einer Eins in der Adresse n entspricht, wird zu diesem Zweck im Speicher 2 bereitgehalten. Vor- 
her muß noch dafür gesorgt werden, daß das Sammelwerk und der Speicher 170 gelöscht sind. 
Damit haben wir die folgenden Befehle, mit Adressen, die beispielsweise mit m beginnen: 


1 Der einfachste und direkteste Weg der Programmgestaltung und Verschlüsselung wäre die Folge 
| 
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Spe = Br Inhalt e N Erklärung 


m To Zwecks Löschung des Sammelwerks | BALL, 

ee AR ade hai) Zwecks Löschung des Speichers 170 | Vorbereitung 
m+2 4170 | Der eigentliche Additionsvorgang; der Befehl im Speicher (m + 3) 
m+3 .. 4100 | wird der Reihe nach zu A100, 4101, 4102, ..., 4140 infolge 


m + Ar 7.170 der Befehle in den Speichern (m + 5) bis (m + 7) 
'M-+5 Alm-+3) ' \ 
"-m+6b "42 . + Adresse im Befehl des Speichers (m + 3) wird um eins erhöht. 


m-+7 Zum 3) 


k Bei hear Wiederholung dieser Folge von Operationen wird die Adresse im Befehl des 
Y “Speichers (m + 3) um eins vermehrt. Vor jeder Wiederholung muß die Maschine jedoch fest- 
‚stellen, daß die Addition noch nicht fertig ist; war die letzte addierte Zahl C (m) so steht, nach- 
dem der Befehl im Speicher (m -+ 6) ausgeführt ist, im Sammelwerk A(n + 1). Ersetzt man 
ie ‚den Befehl im Speicher (m + 7) durch U(m +3), so wird das Sammelwerk nach Übertragung 
seines Ergebnisses nach Speicher (m + 3) nicht gelöscht, es bleibt A(n + 1) in ihm stehen. 
M Ko Wenn man nun A 141 vom Inhalt des Sammelwerks subtrahiert, so erhält man (n — 140), was 
En  megativist, wenn der Additionsvorgang noch nicht fertig ist, null, wenn er fertigist. Es ist zweck- 
. mäßi den Befehl A 141 im Speicher (m + 12) zu haben und den nächsten Befehl, der nach 
 Beendigur cn da ch Ben Ben soll, im Speicher (m + a Dann dienen die len 
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ED Inhalt “ %s.. Erklärung © 

m -+7 U(m +3) tritt an Stelle von 7 (m + 3) der früheren Tabelle, 

m 8 S (m + 12) 

m-+ 9 E (m -+-13) Wenn die Additionsaufgabe beendet ist, dann ist ( (Sm) = 0 
und der nächste Befehl wird aus dem Speicher (m-+ 13) ent- 


nommen. Ist die Additionsaufgabe nicht beendet, so ist C' (Sm) 


negativ und der nächste Befehl wird aus dem nächsten Speicher, 
also (m + 10) entnommen. 


m+10 TO Löschung des Sammelwerks. 


m+11  E(m-+2) Zu diesem Befehl kommt man nur, solange die Alldıore aufgabe 
nicht beendet ist. Es ist C (Sm)=0 als Folge des vorher- 
gehenden Befehls. Deshalb wird der nächste Befehl dem Speicher 
(m — 2) entnommen, um den Additionsprozeß zu. wiederholen. 
* Dieser Sprung, zurück zum Befehl in Speicher % —- 2) geschieht 

bei gelöschtem Sammelwerk. 

m-+12 A14 
m-+- 13 enthält den Befehl, der nach age der a 
ausgeührt werden soll. 


Der Befehl A 141 im Speicher (m + 12) wird hier tatsächlich niemals durchgeführt, er ist 
nur nötig für den von der Maschine ausgeführten Vergleich mit dem Befehl A n, um tenpzu E13 
ob die Additionsaufgabe schon beendet ist oder noch nicht. 

Bei dieser Art der Programmgestaltung und Verschlüsselung der Addition benötigt man 
zur Ausführung der Aufgabe (3) nur 13 Befehle. Das gilt für jede ähnliche Addition, wieviele 
Zahlen auch immer addiert werden sollen. Dieses Beispiel zeigt, daß die Anzahl der Befehle 
nicht unbedingt mit der Anzahl der auszuführenden Operationen wachsen muß. 


Beim Gebrauch dieser Befehle ging die Zahl 100, die Adresse des ersten Summanden, ver- 


‚loren, da der Inhalt des Speichers (m + 3) im Laufe der Rechnung abgeändert wurde. Bilden, 


aber diese Befehle einen Teil einer größeren Gruppe, die im Laufe einer Rechnung auch noch 


wiederholt werden muß, so muß der Befehl A 100 wieder in den Speicher (m + 3) hineingegeben _ 


werden; noch besser speichert man ihn irgendwo anders und gibt ihn durch einen besonderen 


Vorbereitungsbefehl vor Beginn jedes Additionsprozesses an die entsprechende Adresse. 


Man sieht, daß sich einzelne der Befehle des obigen Schemas auf Speichernummern be- 
ziehen, die mit dem Wert m verknüpft sind, d.h. sie hängen von der Stellung im Speicherwerk 
ab, in der diese Befehle selbst gespeichert sind, andere Befehle hingegen sind davon unabhängig. 
Es ist jedoch zweckmäßig, beim Entwurf des Programms und bei seiner Niederlegung im Loch- 


streifen die Befehle so abzufassen, daß sie den Wert m nicht selbst explizit enthalten. Deshalb 


fügt man zu jedem Befehl, der in den Lochstreifen gestanzt wird, einen Schlüsselbuchstaben 


hinzu, und zwar F bei einem Befehl, in dem die Zahl, so wie sie dasteht, die Adresse angibt, auf - 


die sie sich bezieht und © bei einem Befehl, in dem zu der angegebenen Zahl erst noch die Zahl m 
addiert werden muß, um die Adresse zu Arhaften: auf die sich dieser Befehl bezieht. So erhält 


' man für die oben besprochenen Befehle schließlich die folgende Form: nt u 


Speicher Nr. Befehl selbst seine Form auf dem Lochstreifien x 
. @ aiRas br x 1 P 
m PLN TORE NE NE 
m + 1 T 170 a PORN 2 We 
m 2 A 170 AZOREN RE 
m + 3 4 100 A100 ig BEN 
m + 4 T 170 ie 217 T10OF VIEL DEE 
m +5 A (m + 5) rg 
m -+ 6 42 
m + 7 U (m +5) 
m-+ 8 S (m + 12) 
m + 9 E (m + 15) 
m -+- 10 TO 
m + 11 (Em -+2) 
m+ 12,840 4141 
m -- 13 er 
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_ Die Schlüsselbuchstaben @K, F und © erscheinen in den Befehlen nicht in der Form, die sie im 
Speicher haben; sie sind Befehle für die Maschine hinsichtlich der Ausdeutung der Lochungen 
des Streifens und der Bildung der entsprechenden Befehle. Der Eingang @K auf dem Loch- 
streifen befiehlt der Maschine, daß sie beim Ablesen des Streifens und der Bildung der zu speichern- 
den Befehle die Nummer des Speichers, in dem der nächste Befehl niedergelegt werden soll, 
festhält. Das bezweckt, daß die Nummer zur Bildung der Adresse solcher Befehle verfügbar ist, 
die auf dem Lochstreifen durch den Buchstaben © gekennzeichnet sind. (Außerdem dienen 
die Buchstaben F und © zur Markierung des Endes der Zahl, die die Adresse des jeweiligen Be- 
fehls angibt.) 

Sind diese Befehle einmal in dieser Form verschlüsselt, so ist es gleichgültig, an welcher 
Stelle des Speichers sie niedergelegt werden sollen. Wenn sie einmal aufgestellt, gelocht und 
bestätigt sind, können sie kopiert werden, da sie für jede gleiche Rechnung ein für allemal be- 
stehen bleiben. Ein Hilfsgerät erledigt dieses Kopieren mechanisch und prüft die Kopie. 

Es gibt eine wichtige Ausdehnung dieses Gebrauches von Schlüsselbuchstaben. Sie ver- 
langen die Ausführung von Operationen an gelochten Befehlen oder Zahlen, schon bevor diese in 
den Speicher gebracht werden. Damit ist es möglich, den gleichen Satz von einmal ge- 
lochten Befehlen für je de’Rechnung des Typs | 


BLlEr 


zu benutzen, gleichgültig, welche Adressen a, b und c in einem speziellen Beispiel vorkommen. 

Natürlich müssen die speziellen Werte für die Parameter a, b und c der Maschine in jedem Einzel- 

falle mitgeteilt werden. Dies geschieht durch besondere Eintragungen im Lochstreifen, die den 
eigentlichen Befehlen vorausgehen. . 


G 
P- 
> 
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| " Das gesamte Arbeiten der Maschine einschließlich Lesen und Übersetzen der Eintragungen 
auf dem Eingabestreifen kann nur durch Befehle gesteuert werden, die aus dem Speicherwerk 
zu entnehmen sind. Wenn die Maschine eingeschaltet wird, sind die Speicher leer. Deshalb ist 
ein für allemal eine Grundgruppe von Anfangsbefehlen in einer Reihe von Drehwählern (mehr- 
polige Dreh-Schrittschalter) fest verdrahtet. Wenn nun der Startknopf gedrückt wird, werden 
diese Anfangsbefehle in den Speicher übertragen. Sie befähigen die Maschine, die Eintragungen 
des Eingabe-Lochstreifens zu lesen und zu übersetzen. 

Zum Schluß noch eine Bemerkung: Man hat diese Maschinen ‚Elektronische Gehirne“ 
genannt. Diese Bezeichnung verleitet zur Annahme, daß sie selbständig denken könnten, was 
sie jedoch nicht tun können. Sie können nur ganz blind die'Folge von Befehlen ausführen, welche 
für sie ausgedacht worden ist. Wie Byrons Tochter, Lady Lovelace, vor über hundert Jahren 
von Babbages Analytical Engine sagte, ‚‚kann die Maschine keinen Anspruch erheben, irgend- 
etwas zuerschaffen. Sie kann nur ausführen, was wirihraufzutragenver- 
stehen.‘ Und das gilt für die bisher gebauten Maschinen genau so; wie damals, als es ge- 
schrieben wurde, "” 


Ä 
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he eine Methode der Störungsrechnung und ihre 
: Anwendung auf Schwingungsprobleme 3 
Von Rudolf Haag in Stuttgart!) 


Das gewöhnliche Verfahren der Störungsrechnung versagt bei manchen Problemen in gewissen kritischen 
Bereichen. In der vorliegenden Arbeit wird gezeigt, wie in diesen Fällen die Störungsrechnung mit Erfolg 
gehandhabt werden kann. Die Methode wird angewandt auf die Berechnung von Mitnahme- und Kipp- 
schwingungen. Man erhält in verhältnismäßig einfacher Weise die Lage der Mitinahme- bzw. Kippbereiche, 
sowie bereits bei erster Näherung die wesentlichen darin auftretenden Phänomene. f 


The usual method of little disturbances fails in certain critical areas of some problems. It is shown how 
to handle the method of little disturbances effectively in these cases. T’he method is applied to the computation 
of driving and relawation vibrations. The domains, where driving and relaxalion phenomena occur, are j 
obtained in a relatively simple way, and already the first approximation gives the essential effects. 

En tant de probl&mes le procede ordinaire du calcul de perturbation manque dans certaines domaines 
critiques. Dans le travail pr sent on fait voir, comme, en se cas lü, le calcul de pertwrbation peut ätre 
employ& avec succes. La methode est appliquee pour calculer les oscillations d’entra'nement .et de bas- 
culement. On obtient, d’une maniere relativement simple, les domaines d’entrainement et de basculement. 
On obtient, d’une maniere relativiment simple, les domaines d’entrainement et de basculement et, deja 
dans la premitre approximation, les phenomenes essentiels, que s’y prösentent. 


IIpıu pemennun HeKoTopbIXx 3anay OÖBIYHBIH MeToA BOSMyINeHnÜ HENPHTONeH B USBECTHBIX 
KpHuTuyeckux O0nacTax. B HacTonmei paboTe TIOKABEIBAeTCH, KAK B ITUX CHYYAAX MO3KHO 
© YCHeXOM HONB3OBATBCA METONOM BOsmyIennä. Jlo 3Tomy MerTony PacuuTsIBamTea arech 
KOJMeOAHHUA SAXBATBIBAHHA U OTPORHAEIBAHHA. IlyTrem CpaBHeTHIBHO IPOCTOrO pacyeTa HAaXO- 
AUTCA INONOSKEHNE ODNACTEH BAXBATBIBAHHA MIU ONIPOKHABIBAHUA H ya m3 NePBOorO IPuÖan- 
3KCHHA BEIACHAIOTCH CAMBIE CYINECTBEHHBIE ABIICHUS B ITUX OÖNACTAX. Br 
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1. Allgemeines 

Auf allen Gebieten der Physik ist die Zahl der Probleme, die einer EOnlossenen En 
matischen Lösung zugeführt werden können, relativ gering. Bei den ‚übrigen Problemen hilft 
man sich im allgemeinen dadurch, daß man durch Wegstreichen gewisser ‚‚Störungsglieder“ den x 
Anschluß an ein idealisiertes Problem gewinnt, das in einfacher Weise streng lösbar ist undvon - 
da aus Störungsrechnung treibt. Ihre Methode läßt sich ganz allgemein so formulieren: das } 
idealisierte Problem gehe aus dem ursprünglichen dadurch hervor, daß ein Parameter A gegen 
Null strebt. Dann setzen wir die Lösung als Potenzreihe in A an und berechnen rekursiv ihre 
Koeffizienten. Es sind nun häufig gerade solche Fälle von Interesse, in denen’ größer ist als 
der Konvergenzradius r der Potenzreihe, denn erst dann sind wesentlich andere Phänomene als 
bei dem idealisierten Problem zu erwarten. Vom Be der Funktion ntheorie aus handelt 


Fällen hilft dagegen eine gewisse nrendeas des ER, der Be gsrechnung. zum ER 
folg, das durch folgende Überlegung nahe gelegt wird: Der Konvergenzrad en ab vo h me RR 
Lage der übrigen Parameter des Problems — die ich BR durch ‚pchar isieren 1 a) t 
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für die Störungsrechnung erhält, wenn man statt p den Parameter q einführt, der die beiden 
„kleinen Größen“ (p— pr) und A koppelt: 


P—-NM)=g'X ee a (LE 


Die Lage dieser kritischen Stellen ist meist leicht anzugeben, da sie eine unmittelbare physi- 
kalische Bedeutung besitzen. Bei Schwingungsproblemen sind es die Resonanzstellen. 

Die hier skizzierte Methode eignet sich vor allem für Schwingungsprobleme und einige 
Fälle der quantenmechanischen Störungsrechnung (z.B. bei der Berechnung der Elektronen- 
zustände im Kristallgitter oder bei Fast-Entartung). Sie soll im folgenden auf die Berechnung 
von Mitnahme- und. Kippschwingungen angewandt werden. 


2. Schwingungen in einem Freiheitsgrad. Überblick und Ansatz der Störungsrechnung 


Das idealisierte Problem, von dem wir im folgenden stets ausgehen, ist die Gleichung der 
harmonischen Schwingungen 


BF DBEL.DSBZDK, 608 (iO). ee er) 


Was ist zu erwarten, wenn sich eine der Konstanten mit dem Ausschlag oder mit der Zeit ändert ? 
Es ist für die Art der auftretenden Phänomene unwesentlich, welcher Parameter in Gl. (3) ver- 
änderlich ist (C. Schäfer, Ann.d.Phys. Bd. 33, 1910), wir können uns dabei auf den Fall 
veränderlicher Eigenfrequenz beschränken. Bei der Zeitabhängigkeit interessiert vor allem 
der Fall der periodischen Schwankung, bei der «-Abhängigkeit beschränken wir uns auf die 
Berücksichtigung des nächsten Gliedes der Potenzreihe. Das ‚‚gestörte‘‘ Problem lautet also 


+ Pat oc +0oco(pt), HA —= 3 K,cos (wt—9Q) - : (4). 


In nullter Näherung erhalten wir eine Überlagerung harmonischer Schwingungen mit den 
Erregerfrequenzen w; und der Eigenfrequenz o des Oszillators. Die Parameterschwan- 
kung bewirkt, daß wenn in m-ter Näherung eine Frequenz 2 vorkommt, in (m + 1)-ter 
Näherung die Frequenzen Q +» auftreten. Die Nichtlinearität erzeugt zu je zwei 
Frequenzen 2, und 2, in nächsthöherer Näherung die Frequenz Q,—+ Q2,. Um zu einem ein- 
fachen und konsequenten Verfahren der Störungsrechnung zu kommen, empfiehlt sich der 
folgende Ansatz: 

1. An Stelle der goniometrischen Funktionen verwenden wir durchweg die komplexen 
Exponentialfunktionen. Dies hat — abgesehen von der größeren Durchsichtigkeit der Rechnung 
— noch den Vorteil, daß instätionäre, exponentiell an- oder abklingende Schwingungen einfach 
dureh komplexe Frequenzen dargestellt werden können. 


2. Wir setzen als Lösung an 


BeuheH ... 2.2.2.0 n2eneß) 


f: wobei die 2; diejenigen Frequenzen sind, die aus den Grundfrequenzen ® und w; des idealisierten 
Problems durch die Näherungschritte erzeugt werden. Setzt man diesen Ansatz in (4) ein, so 
erhält man ein (unendliches) Gleichungssystem für die Ar. 


3. Wir wollen nicht den Anschluß an ein spezielles Anfangswertproblem, sondern das 
allgemeine Integral von (4) wiedergeben. Daher lassen wir A, in (5) als willkürliche komplexe 
” Integrationskonstante und führen die übrigen Amplituden darauf zurück. Dafür müssen wir 
zulassen, daß die Eigenfrequenz 2, = w gegenüber dem ungestörten Problem verändert, also 
ebenfalls zunächst unbekannt ist. 


4. Wir setzen alle Az (k #0), sowie w als Potenzreihen i in den Störungsparametern A und 
(oder) « an. Im Fall «= 0 etwa 


i Er KeSAPR; SE BPELE (6% 


Damit läßt sich das Gesamtsystem rekursiv berechnen. 

Zu dieser Gebrauchsanweisung ist für den Fall der nichtlinearen Schwingungen noch 
eine Bemerkung zu machen. Während bei linearen Problemen (A= 0) Real- und Imaginärteil 
einer komplexen Lösung je selbst eine Lösung darstellen, ist dies bekanntlich bei nichtlinearen 

Problemen nicht mehr der Fall. 

Man muß daher für sie von vornhereineine reelle Lösung anstreben, in dem man in (3) zu 
jedem ‚Glied A,ei®st das konjugiert komplexe ebenfalls aufnimmt. Sind alle Frequenzen 2; 

.. reell, so bedeutet dies keine nennenswerte Komplikation. Dagegen ist für die nichtstationären, 

nichtlinearen Schwingungen (komplexe Eigenfrequenz) der obige Rahmen ein Prokrustesbett. 

Der Barde. Ban dafür ist, Sy bei nichtlinearen Schwingungen allgemein die Eigen- 


wie E x en 
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frequenz amplitudenabhängig ist. Bei instationären inne ändert sich aber die Am- 
plitude mit der Zeit, so daß auch die Eigenfrequenz zeitabhängig wird. Der Ansatz einer kon- 
stanten Eigenfrequenz ist daher dem Problem nicht angemessen. 


3. Lineare, freie Mitnahmesehwingungen r 
Wir diskutieren hier Gl. (4) für den Fall A=0, k;—= 0, also 


But £+B&-+o} ct+o(erite-w)a—=0 . DER RER EDER HR 
Re Im allgemeinen erhält man bei kleinem & einen, gemäß der Dämpfung ß abklingenden, nahezu a 
BR periodischen Vorgang, der nur durch das schwache Auftreten der erwähnten ‚‚Variations- ie; 
Be: [requenzen“ &-+ kp gegenüber dem der rein harmonischen Schwingung etwas modifiziert ist. 
| Bekanntlich stellt aber die periodische Änderung der Streife unter geeigneten Bedingungen eine a 
h dauernde Energiezufuhr dar, durch die in diesen sogenannten „Mitnahmebereichen“ total andere 
‚FOR Verhältnisse geschaffen werden. Wie spiegelt sich dieser Sachverhalt in unserer ‚Störungs- ig 
BB rechnung ? £ 3 
IR Wir machen den Ansatz von Gl. (5), wobei wir nach dem oben gesagten alle Frequenzen B- 
2 
A=o+k'p Et) si 
zulassen müssen. Gl. (7) reduziert sich dann auf ER 4 
BR, Ar + a (A-ı + Arıı)=0 RE je & Erge n eline (9), r 
wobei R; der Resonanznenner der Frequenz 2; ist. RE Ir pi 
| B-otPBR—M 2.2... (0). 2 
Normale Störungsrechnung u 
Machen wir den Ansatz (6) (mit « statt A) und beschränken uns au PeStrE Potenzen von &, Sr 
so folgt zunächst aus (9) für die nullte Potenz vonabik=0 
2 RO=0. ER TS R BR 
d.h. i ER a a 
e' (0) — HR *4 FE 2 ß? Ey 3 ® 
[0) of, mit 2 Aare EM . a2). 5 
f® Welches Vorzeichen wir rählen, ist gleichgültig, wir bentitzen das positive, Ren e 0) und | 
(12) erhält man für Rz: 
Rı po) 2 Ha HEN .@s) N 
Damit liefert Gl. (9) für die einzelnen Potenzen von &: FR RR 3 
Potenzeny; | kp@atkp) AP 0. . Er BR Bert, 4 
Da die Klammer nie verschwindet, hat man A j2 REN : END 15 Y E 
BE . en An Or ee AN = —A bel a  \ 
P.ote'n zenüat: | ee PR a EN 
.hr70: kp@a+ kp) AD = AP. a ah BETEN s: 
ADS ar ei Sen 
” ' A h a 7 - EA AN 
„ R am = Tr Zterr 4% = a N = re 
| 0 P@o+p) TINO | 
Bu: | ERRER 


200 WA=0; = ie: 


Potenzen.a: 


[4 


BE | +0: kp(20+ kp) AP 
Er i 2-0 für 


Bor ni, 200m4- AHA 


> 
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ar 


f. 


Man sieht, daß sich ein reeller Zusatz zur Eigenfrequenz ergibt. Die Entwicklung läuft reibungs- 
los so weiter. Sie versagt jedoch, wenn für irgendein % (es sei für k—= — m) der Resonanz- 
nenner sehr klein wird, wenn also 


m DERO 20212.10.(40) 


wird. In diesem Fall werden die Glieder in der Reihe für A, mit jedem Näherungsschritt größer, 


Abgewandelte Störungsrechnung 
rn den Divergenzgebieten der normalen Störungsrechnung springt die abgewandelte ein. 
Wir wollen sie für.den ersten kritischen Bereich durchführen. Es sei also 28 — p=0. Setzen 
wir nach dem Muster von (1) 

N. 2. ne. (20) 
so wird 
/ Rr=—4o:k(l-+k)— oa(o@V! -+kg)2w(1+2%) | 

— fo 2ol +2) Hd Hk) +... 2:2... (21). 


‚Es wird also R_ı —od (das Zeichen soll die Größenordnung andeuten). Da A, wo ist (als 


Integrationskonstante) folgt, daß auch A_,— a’ wird. Die übrigen Amplituden sind wieder 
von kleinerer Größenordnung, so daß wir keine negativen Potenzen von & benötigen. Im einzel- 
nen liefert dann Gl. (9): 


Potenzen oo: AD—=0 für k#0,; :.k#—l. 
Dotenzen co; 

k=0: . . 200) A= A, 

k=—1: („U —)2A, = —A. | 


. . . . . (0) Be . 0 . 
Dies sind zwei Gleichungen zur Bestimmung von o"!) und A"), . Elimination von A, ergibt 


00 late ;) EA Eh (22). 


Die Störungsrechnung läuft in dieser Weise schematisch weiter. Sie ist etwas mühsamer als im 
Normalfall, da immer gelegentlich Gleichungen mit zwei Unbekannten auftreten. Aus der 
une anten Näherung läßt sich jedoch alles Wesentliche ablesen. Man erkennt, daß Ir 


l<- — also RE | 
De A . (23) 


ein imaginärer Zusatz zur Eigenfrequenz entsteht (mögliche Enddämpfung). Die Schwingung 
enthält in diesem Bereich zwei Frequenzen stark: und @— p. Der Realteil dieser Frequenzen 
ist in erster Näherung: / 


Re (0) = Re(o + a0) = 5 EN NE ee AR (24). 


(Es läßt sich übrigens zeigen, daß die höheren Näherungen an dieser Gleichung nichts mehr 
ändern, da nur rein imaginäre Zusätze zur Eigenfrequenz auftreten.) 


Man hat also folgendes Bild: In dem durch Gl. (23) gegebenen Bereich ist eine der beiden 
Epektallösnngen entdämpft mit der ne (Gl. (22)) 


SEN Var u_, A ER On (25); 
für i 
en SEN ER (26) 


klingt die Schwingung zunächst exponentiell an, bis die Amplitude durch irgendeinen, nicht in 
GI. (7) berücksichtigten, Faktor begrenzt wird. Von da ab haben wir eine stationäre Schwingung 


KR 


_ mit der durch (24) gegebenen Frequenz. Die Schwingung ist sogar rein periodisch, da die 


RN =: $£ + kp höhere harmonische von ® = r sind. Ihre Grundfrequenz ist mit der Frequenz 


de Paramelerschwankung: Bene im Verhältnis 1:2. Daher die Bessiehuung „Mitnahme- 
v BayBung i 


. 


Diehöheren Resonanzen 
Nach genau demselben Muster lassen sich die Divergenzen des normalen Verfahrens bei 


2@ mp behandeln. Man muß dabei (bei jedem m) ansetzen: 2 
20 1‘ 
Doz ir ZT a ER ENN e - (27). 


Es ergibt sich, daß für m>2 keine Enddämpfung mehr auftritt. Dagegen erhält man 
bei m = 2 einen Mitnahmebereich ö 


Hi, 8 4 1 / 
und Instabilität für: % } | FA 
a2 w® & 112 *2 

= 1-4, @-0+3| SB. ee / 


Bemerkungen zumAnsatzvonGl. (l, Konvergenz 

‚In Gl. (1) war der Exponent y noch unbestimmt gelassen worden. Welche Gesichtspunkte 
Ha hat man für seine Wahl? Es gibt zwei: Erstens gibt es einen Wert y’ derart, daß für y> y’ keine 
[ES rekursive Berechnung des Gleichungssystems mebr möglich ist. Dieser Wert ist leicht zu 
ermitteln, indem man mit dem Ansatz (1) in der Ausgangsgleichung (in unserem Fall Gl. (9)) 
R die Größenordnungen der Anfangsglieder abschätzt, wie es etwa auf Seite 15 geschah. Mit 
y==y’ erhalten wir dann die gewünschte abgewandelte Störungsrechnung. Wenn es möglich ist, 
die Konvergenz des normalen Verfahrens abzuschätzen, so gibt es einen zweiten Weg. Es sei 
etwa: 

"Agrenz D% (p vr Pr)” e 
Dann führt der Ansatz ! 4 
?—- M=1' RE EN ENG art (30) 
zum Erfolg 1). In dem behandelten Fall ist eine solche grobe Abschätzung durchführbar. Schreibt 
man (9) in der Form 
PAARE L.) Arie, 

so sieht man, daß ein Amplitudenglied der Schwingung Q; nach zwei lei wieder 
ein Glied dieser Schwingung erzeugt, durch Kopplung mit den beiden Nachbarschwingungen 
Q4 1 sowie durch den Ausdruck RP AU”. en wir nur diesen unmittelbaren ‚Zusammen- 
hang der Amplitudenkoeffizienten A’ und AÜ”” ins Auge fassen, so wird 


AN RD RR 1 RE 
ee + RD Rn, oz RO Rf RO, © . . . ea WA EN. 


Die Entwicklung wird divergieren, wenn irgendein R; klein ist, also Gl. (19) bestähtt Um ‚die 


Breite des kritischen Gebietes abzuschätzen setzen wir Br, ‚n|= v Sn | ==, 


bi 


RW 2 me RE Kat Aue Se 
RV ER La 9 | R® 1 0) . : un 2 ER: 
—m+ı = (m IP m—ı © (m+1)p? (> Rdn) pa TR SE 


Dann wird en GERN ER N 


1 m+1 1% rel Re 
Pre 
— mM % TR # t nn, a ne D 
Man hat somit als rohes Konvergenzkriterium beimzi Bf i 
Alm A pr a? : a2 FR 
na " ormm— EM 


a? 
Egrenz I Dan En op m(m 1) ie 


1) &, bzw. a, eiP ist eine singuläre | ‚Stelle der Lösung, di A 
(n negativ oder gebrochen). Mit der Substitution in (30) wir: 
rückt und die Trennung der Bee &-Ebene in zwei ERE 
Ben. wegfällt. u... 


de! En Er | Ko, r 
. 2. ew.M: th. Mech. ö > fe: 24 
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Bei m= 1 fällt wegen AV) = 0 das dritte Glied auf der rechten Seite von Gl. (31) fort und das 
erste Glied ist das entscheidende. Hier ist nun: 


h FR 
2 1 An 1 
Rn 2 ae also era, 2 
direkt 


40 —p? ep ep® 
& E 
Egröng 5 ES 5 Dem le u a er (33). 


Man sieht, daß diese Abschätzungen die Größe der Mitnahmebereiche — von Zahlfaktoren ab- 
gesehen — richtig wiedergeben und zu den Ansätzen (20) bzw. (27) führen. 


4. Erzwungene niehtlineare Schwingungen 
Wir diskutieren die Gleichung 

stBstosti®=KlevtiLeriwm) ....n... 1 2 
3 und beschränken uns dabei auf den stationären Anteil, sehen also ab von der Eigenschwingung 
und ihren möglichen Instabilitäten. Dann läuft Gl. (5) auf den einfachen Fourieransatz hinaus: 


4 a Sue (weckte) mi = A a... (80) 
und Gl. (34) liefert 
BA SA rt) wobei dm für 17 
mit - i 
Keen -Findwnn Ne lern (37). 


Normale Entwicklung 


Die in R; vorkommenden Größen sind alle DekaunE Mit dem Ansatz (6) für A, erhält 
man RER \ 


2 (0) x 
| dr = für k#-+1 
ä AD — Be I 
R 

I ER N EP (BBR 
K A für 20, +2 
k Br, | 24. wn p) K:? a) ae _ ae K:2 

Say ER TEN Sets Ar AR ee DE 
; a BISLS NEBEN: R, RiR, 
Wiederum versagt die Entwicklung, wenn einer der Resonanznenner klein wird, d.h. wenn 
ö Be nun, BE NER (39). 
67 Eine Kon abachätzung zeigt, daß dies bei Resonanz mit der Grundschwingung (m = 1) 
R der Fall ist, wenn ER oG@ ) 
E- ‚BR, 0 @—u) Mae Be De 4 
“ k Az K? K: (40) 
| ‚bei Resonanz mit der ersten Oberschwingung (m = 2), wenn 
E: ER i Ne BETEN, NT . (41) 


ER Ru. Ki 
RE (Von Zahlfaktoren ist abgesehen und die Dämpfung ist als hinreichend schwach voraus- 
gesetzt. ) 


u Entwicklung für dieerste Resonanzstelle 
Br Wir setzen 
RE: ee Kar ik! We ee auäen. . (42), 
Y = sowie u pn 6 wegen «) R, »42!3 werden soll) und entwickeln nach Potenzen von 
% iz 2 An 
O1. O0) schreib sich dann: 
7 An Brise Sa, Ari = Kl +. sa TPIR VE VPE (43) 


Er, BARS 7 
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Wie man sieht, erhält A,ı einen zweifachen Pol in A’ —0. Damit erhalten A, und As ein- 
fache Pole. Für die Anfangsglieder der Entwicklung bei k = 1, 0, 2liefert Gl. (43) die Gleichungen 


2 
AD H2IDAD-O | RE LAH): 
3010 + (AD) =0 | 


Daraus lassen sich die drei Amplitudenkoeffizenten AT®, ASP, AV) berechnen, wenn man 
bedenkt, daß A_; = 4£ ist. Wir setzen: 


— 209" JA9 + 240 An +24 VAR | 


AN = Re NR ae ae ee u 
und erhalten nach Elimination der übrigen Amplituden 
10 ” iß" 
art2alg)r-—Kei N 2 aa nn 


Multiplikation jeder Seite mit ihrem konjugiert Komplexen führt auf die Gleichung dritten 
Grades für r?: | 


(rt) tagte ro...) 


Zu jeder positiv reellen Lösung von (48) gibt es eine von (47) und damit auch eine Lösung des 
Ausgangsproblems. Zunächst können wir die Grenze des Bereiches bestimmen, indem (48) drei 
reelle Lösungen besitzt. Auf dieser Grenze müssen zwei (reelle) Lösungen von (48) zusammen- 
fallen. Die Bedingung dafür ist bekanntlich das Verschwinden der Diskriminante. 


Mit der Abkürzung 
; 5 k242\1/8 
ale) . EN) 


erhält man nach einiger Zwischenrechnung für diese Bedingung: 


oe ler er leer ner 
N ee | Sin re —)| 1 2 0.2.. } 
Be Urea. El bau] 0) 
= & 2 als Ordinate, 2 als Ab- 
szisse. Aus (48) erkennt man, daß für 8 = 0 (#’—=0) und © — w= 0 (g= 0) nur eine reelle 
Lösung existiert. Der Bereich mit drei reellen Lösungen ist also der Innenbereich des durch die 
beiden Grenzkurvenzweige gebildeten Keils. Ferner sind 
die drei reellen Lösungen in dem Keilgebiet sämtlich posi- 
tiv, so daß ihnen auch drei Lösungen unseres Ausgangs- 
problems entsprechen. Die bekannte Tatsache, daß bei 
nichtlinearen Schwingungen in der Resonanzgegend ver- 


Die Gestalt dieser Grenzkurve ist in Bild 1 aufgezeichnet mit 


2b 
20 EERTIIIERTITIINIITIT) „ 


28 02 04 08 08 mW 12 4 71 78 


BAER 77 > WM, ö ns TR k 
Abb.1. Die Grenze er Abb.2. Typische Resonanzkurve der Schwingung (34) Fe? 
nach G1.(50) 2:75,55 8.74 a 
A N 1: 


schiedene stationäre Schwingungsformen möglich sind, ergibt sich also zwangslos aus unserer 
Störungsrechnung. Eine typische Resonanzkurve der Schwingung (34) zeigt Bild2. Sie gibt de 

Amplitude der Grundschwingung in Abhängigkeit von der Erregerfrequenz bei festem o, ß, Au. 
Aus Bild 1 können wir nun bei beliebiger Lage aller Parameter den Mehrdeutigkeitsbereich a 


.. ) Für w<w<w, sind drei Lösungen vorhanden, von denen allerdings die mittlere nicht stabil Pe i 
Bei Anderung der Erregerfrequenz springt die Amplitude bei w, vom unteren Ast zum oberen, bei w, ? 
oberen zum unteren (Kippschwingung). & Be N DE ENTER 


el 


i 
| 
F 


”. 
ge 
& 
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lesen. Man sieht, daß es einen kritischen Wert der Dämpfung gibt, jenseits dessen keine Mehr- 
deutigkeit mehr auftritt. Bei fehlender Dämpfung erstreckt sich der Mehrdeutigkeitsbereich ins 
Unendliche, die Resonanzzacke in Bild 2 ist offen. 


Entwicklungfür die zweite Resonanzstelle 
Wenn © = 2 w, also R, »- Oist, müssen wir auf Grund des Ergebnisses (41) der Konvergenz- 
betrachtung den Ansatz R, — A, also 
vw—@o-+ qAR; p=Pß=P'i4R; Be NEE REN 


als Ausgangspunkt nehmen. 
- Damit ergeben sich folgende Größenordnungen 


Beer Ama AS, Am AUS, 


(Es sei übrigens nochmals darauf hingewiesen, daß man die Größenordnungsbetrachtungen auch 
zur Ableitung der Ansätze (51) und (42) benutzen kann.) Es wird 


2 k2 “3 > 
eg ‚4 
Rı.=o® (1 n) 4 (#4 9 ko-+. 


Für die Anfangsglieder der Entwicklung erhält man - 


To84D—K 


aa Ar —=0 
en e 
= (2 ah) 20" HAN +27 +2 MAD =0 


RE a 


Setzt man AU = rei? und eliminiert die übrigen Amplituden, so ergibt sich 


+40 CW 
| Er + q wr= N a DR" sencie R 
also eine zu (47) vollkommen analoge Gleichung/ In ähnlicher Weise lassen sich auch die 
höheren Resonanzstellen durchrechnen. 


im=3 führt der Ansatz R, A? zum Erfolg. 
„Be 3 


Der Einschwingvorgang bei erzwungenen nichtlinearen 
Schwingungen 


Bei nicht zu schwacher Dämpfung läßt sich auch die Eigenschwingung des durch (34) 
gegebenen Vorgangs in.dem hier zugrunde gelegten Rahmen durchrechnen. Man findet, daß 
sie im allgemeinen mit der Dämpfungskonstanten ß abklingt, daß jedoch in den Resonanz- 
gebieten, die hier durch 1» = nw gegeben sind, unter Umständen eine Entdämpfung auftritt, 
wie wir sie bei den Mitnahmeschwingungen diskutierten. Auf eine Durchführung dieser Rechnung 
sei indessen verzichtet, da (wie auf S. 14 bemerkt) der Ansatz einer konstanten Eigenfrequenz 
dem Problem nicht gemäß ist. 


Herrn Prof. E. Fues möchte ich für die Anregung dieser Untersuchung herzlich danken. 


Eingegangen am 15. Mai 1949, 


2* 
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Der rotierende Tragflügel als Strahlungsproblem * 
Von Wilhelm Ernsthausen 
Mitteilung aus dem Helmholtz-Institut, Nußdorf/Inn 


Die auf einen, im er Gas rotierenden Tragflügel ausgeübten Kräfte entsprechen der Rück- 
wirkung seines Schallfeldes. Sie lassen sich daher aus*dem Energie-Inhalt dieses Feldes bestimmen. 


The forces exerted upon an airfoil rotating in a compressible gas correspond to the reaction of it’s field 
of sound. Therefore they are determined from the capacity of energy of this field. 


Les forces op£rant sur une aero-aile tournant.dans un gaz compressible correspondent & la reaction de son 
champ acoustique. Par consequent ils sont determinables de la capacite de l’energie de ce champ. 

Cap, AeaCTByPINHe HA KAKYW-IHÖO HeCyINYP IIOBePXHOCTb, BPAINAWINYIOCH B CHKHMACMOM 
Ta3e, COOTBETCTBYPT PeaklluH 3BYKOBOTO TIOIA, BO3HUKAWINETO IIPN 9TOMM Bpamennn. Benen- 
CTBHe 9TOTO YRASAHHBIe CH’IBI MOTYT ÖBITB ONPpenelleHkl H3 KOAMYeCTBA HHEePTHH HTOTO NOIA. 


Einleitung!) 


Der rotierende Tragflügel?) erzeugt im kompressiblen Gas ein Schallfeld. Das gleiche Feld 
kann man sich von einer geeigneten Anordnung raumfester Strahler erregt denken. Es ist in 
vollkommener Weise nachgebildet, wenn durch die Eigenschaften der Ersatzstrahler alle Wir- 
kungen des bewegten Körpers auf das Gas erfaßt sind. Die Schallfeldkräfte an der Strahlerfläche 
stehen dann in eindeutigem Zusammenhang mit den zeitlich unveränderlichen Kräften, welche 
die kompressible Strömung am Tragflügel selbst hervorruft. Im Bilde der Strömungsphysik 
erklärt sich der Einfluß der Kompressibilität mit der Bewegung des Tragflügels in der eigenen 
Schallatmosphäre. Denn die von seiner Bahn ausgehenden wellenförmigen Störungen bringen 
durch ihre Interferenz am Tragflügelkörper die gesamte befahrene Kreisfläche zur Wirkung. 


Der folgenden Untersuchung liegt daher der Gedanke zugrunde, unmittelbar vom Schallfeld 
aus in die Zusammenhänge von Strömungs- und Wellenvorstellung einzudringen. Sie befaßt 
sich mit der Bestimmung der am Tragflügel angreifenden Strömungskräfte aus der Energie des 
Schallfeldes, wobei die den Auftrieb beschreibende Massenbewegung im Nahfeld und die für den 
Widerstand verantwortlichen Strahlungsverluste im Fernfeld ihren Ausdruck finden. 


-” 


Die Schallfelderregung 

Ein bewegter Körper kann in zweifacher Weise auf ein reibungsfreies Gas einwirken, einmal 

über die durch eine zirkulationsfreie Umströmung hervorgerufene Druckverteilung und ferner 

durch konzentrierte Kräfte, die außerhalb eines endlichen Gebietes um ihren Angriffspunkt 

verschwinden. Diese sind als Auftrieb an eine Zirkulationsströmung gebunden, welche beim 

Tragflügel als Folge seines Anströmwinkels entsteht. Die damit verbundenen Gasbewegungen 

führen zu Strahlungsverlusten und erzeugen Rückwirkungen in Form von Massen- und Feder- 
widerständen. 

Es läßt sich zeigen, daß an jedem Punkt der durch den umlaufenden Flügel beschriebenen 

„ Fläche — der raumfesten Strahlerfläche — die von der zirkulationslosen Umströmung herrührende 


“ Druckänderung durch die Wirkung einer ungerichteten Quelle — eines Nullstrahlers —, die 


konzentrierten Kräfte durch die Wirkung eines strahlenden Dipols dargestellt werden können?). 

Bild 1 soll die Vorgänge an der Nullstrahlerfläche, Bild 2 an der Dipolfläche veranschau- 
lichen. 

In 1a bewege sich längs der &-Achse zirkulationsfrei ein zylindrischer Körper von profil- 
ähnlichem Querschnitt mit der Geschwindigkeit v. Wir stellen uns, wie in 1b dargestellt, vor, 
daß in der Bahnebene am jeweiligen Ort des Körpers sehr kleine, gegenphasig und entsprechend 
der Br hundung zeitlich ausgelenkte Kolbenstrahler wirken. Sie erregen, da sie immer 


er a AT 
EL “ : 


4 FR der Flü 
Zi bahnen 
BL een. b Beben immer mit einem festen Kern, dem Tragflügel, verkettet (Bild 12). Die vom 
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Kein gegen die Wellenlänge des Schalles sind, ein Kugelfeld. Jeder dieser Strahler erfährt den 
gleichen durch die Körperform gegebenen zeitlichen Verlauf d®(t). Ist der Vorgang periodisch, 
läuft also eine regelmäßige Folge von Körpern durch das Gas, so läßt sich die Zeitfunktion der 
Erregung d°(t) in harmonische Kom- 


v—- 
ponenten zerlegen. Das hierbei auf- x 
tretende konstante Glied verschwin- ? 

det bei der Bestimmung der für die Bild 13 | 


Kraftwirkung maßgebenden Ge- 
schwindigkeit. Somit kann die Form BEER, 
und die Zahl der an einem Umlauf be- 


teijligten Flügel auf das Fourier-Spek- AR 
trum des Erregungsvorganges zurück- —ef Eos x 
geführt werden. Die Strahlerphase Budr6 


schreitet mit der Geschwindigkeit v 
fort. Wird die Unterteilung genügend rt 
klein gemacht, so muß die Anordnung 
1b das gleiche Feld wie der Körper 
selbst hervorbringen. 

Es ist ohne weiteres möglich, das R 
Feld der Ausweichströmung (bewegter id Te 
Körper, ruhendes Gas, ruhender Beob- 133° 
achter) um einen solchen Körper, wel- Bild 1. Bewegter Körper ohne Zirkulation 
ches vereinfacht durch das Feld eines 
stationären Dipols darstellbar ist, mit dem Kugelfeld der Ersatzstrahler in Beziehung zu bringen. 
In unserer Darstellung kommt dieses Dipolfeld dadurch zustande, daß infolge der Phase 
zwischen den benachbarten Strahlern jeder einzelne dieser Strahler während seiner Auslenkung 
aus der Ruhelage hintereinander den Quellen- und Senkenzustand durchläuft. Der feste Phasen- 
abstand der Strahler untereinander ergibt dann, da auf der Vorderseite der fortschreitenden 
„„Körperwelle‘‘ der Quellenzustand, auf der Rückseite der Senkenzustand überwiegt, das Dipol- 
strömungsbild. Der durch die in 1b eingetragenen Pfeile angedeutete Geschwindigkeitsverlauf 
ist für die Ober- und Unterseite in Bild lc dargestellt. 

Während die Strahler selbst nur Bewegungen senkrecht zu ihrer Fläche ausführen, bewegen 
sich die Gasteilchen auf den radialen Bahnen ihres Kugelfeldes. Durch phasenrichtige Über- 
lagerung der den Strahlern zugehöri- 


EA) " y. 
al ai 
ZN 


gen Bewegungen erhält man die nor- 2 PER 
malen und tangentialen Bewegungs- var ae 
und ° Geschwindigkeitskomponenten ! x 
der Gasström } En 
an Bild 22 ae 
Um die Vorgänge an der Dipol- R ne 


fläch® zu behandeln, gehen wir nach 
Bild 2a von der Zirkulation J' eines 
‚angestellten, längs der x-Achse be- 
wegten Profils aus. Die aus der Strö- 
mung an seiner Ober- und Unterseite 
resultierende Verteilung der Normal- 
geschwindigkeit, wie sie qualitativ 
auch dem Zirkulationsvorgang nach 
2a entnommen werden kann, zeigt 
Bild 2b. Die gleiche Geschwindigkeits- 
verteilung erhalten wir, wenn wir jetzt 
in den Bahnpunkten des Tragflügels 


U’(x)=U’t) 


. \ 
* * . * . [ 
kleine Dipole (in Form frei schwingen- Bild 2c % SA 
der starrer Scheiben) so arbeiten I. 
lassen, daß sie die einmal ausgebildete Bild 2a—2c. Bewegter Körper mit Zirkulation 


'Zirkulationsströmung mit der Be- 
wegungsphase weitergeben. So-erhalten wir Bild 2c. Der Anströmwinkel a* ist damit in 
die Auslenkung dt Mer elementaren Membranen übergeführt. 

"Den Vorgang d!(f) zerlegt man ebenfalls in harmonische Komponenten. Der zeitlich ver- 
änderlichen Auslenkung entspricht im Bild des angestellten Tragflügels die Verteilung der längs 
geltiefe wirkenden Kraft. Ihre Teilschwingungen verlaufen sinusförmig auf den Kreis- 
der Flügelelemente. Die dabei gebildeten, den Auftrieb ergebenden Zirkulationsströ- 


1 
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Tragflügel in jedem Bahnpunkt ausgelenkten Gasteilchen können ihre kinetische Energie an 
diesen nicht mehr voll zurückgeben, da er weiterläuft. Deshalb entsteht ein axial gerichteter 
Gasstrom, dessen Reaktion der Auftrieb ist. Er steht unmittelbar mit der Dipolamplitude in 
Beziehung und ist bei kleiner Flügelzahl ungleichmäßig über die Strahlerfläche verteilt. 

Zwei Versuche sollen qualitativ unsere Annahmen stützen. Die Versuchsanordnungin Bild 3a 
besteht aus einer drehbaren Scheibe, auf der eine Rippe und eine gleichartige Mulde angebracht 
sind. Mit Hilfe eines vor der Scheibe befindlichen Druckempfängers kann der Druckverlauf 
phasengetreu aufgezeichnet werden®). Markiert man die Zeitpunkte, zu welchen sich Rippe bzw. 
Mulde vor dem Druckmesser befinden, so ergibt sich nach Bild 3b an der Mulde Überdruck und 
an der Rippe Unterdruck. 


Bu Druck 


Sog 


Bild3b 


Druck 
| 
Sog 


Auslenkun 
zum Mikrofin 


I 
Rıppe 


Bild3a Bild3c 
Bild 3a—c. Vergleich der von einem bewegten Körper und einem feststehenden Strahler erregten Druckfelder 


Zum Vergleich wurde der Druckverlauf an der Membran eines Lautsprechers gemessen. 
Die Schallwellenlänge war groß gegen den Membrandurchmesser, das Schallfeld kugelförmig. Die 
unsymmetrische Bewegung der Membran wurde durch einen geeignet verzerrten Strom hervor- 
gerufen und zugleich mit dem Druckverlauf registriert (Bild 3c). Bewegt sich die Membran auf 
das Mikrofon zu, so entsteht an ihr Unterdruck und umgekehrt. Dies ergibt auch die Theorie bei 
gekrümmten Wellenflächen. Danach ist der Druck auf einer schwingenden Kugel, die klein 
gegen Aist, so daß die Massenkräfte die Strahlungskräfte übersteigen, gleich der Beschleunigung 
der Gasteilchen®). Er steht also in Gegenphase zur Wegamplitude der Kugeloberfläche. In 
einem ebenen Schallfeld dagegen treten nur Strahlungskräfte und keine Massenkräfte auf, der 
Druck ist in Phase mit der Geschwindigkeit. Ein Strahler endlicher Abmessung nähert sich mit 
steigender Frequenz diesem Zustand, weil die mitschwingende Gasmasse immer mehr abnimmt. 


„Nach Bild 3 entspricht die Membranbewegung der flächennormalen .Auslenkung der Gas- 
teilchen, welche durch die Strömungsvorgänge an Rippe und Mulde entsteht. Zwischen dem 
Druckverlauf an einem festen Strahler und dem an einem rotierenden Körper kann daher kein 
grundsätzlicher Unterschied bestehen. 


Das Verhalten des zirkulationsfrei bewegten Körpers oder eines örtlich feststehenden Null- 
strahlers wird (in Strenge nur bei sehr großem Scheibenradius) durch die rotierende Rippe be- 
schrieben. Ein dipolähnliches Gebilde oder einen Zirkulationskern erhält man, wenn man Rippe 
und Mulde vereinigt. Dies würde etwa eine Form, wie sie in Bild 2c durch die Anordnung der 
Dipolmembranen entsteht, ergeben. Hier kann sich eine Zirkulationsströmung entwickeln, da 
die Bewegung der Elementarmembranen in der Bahnebene gleich derjenigen der Gasteilchen ist. 


Die bisherigen Überlegungen lassen sich wie folgt zusammenfassen: Ein rotierender Trag- 
flügel kann durch raumfeste Strahlerflächen, Nullstrahler- und Dipolflächen, ersetzt werden. 
Ihr Interferenzfeld und die Verteilung der Normalgeschwindigkeit bestimmen Größe und Verlauf 
der auf ihn wirkenden Kräfte. Die Nullstrahlerfläche führt gegenphasige Bewegungen aus, 
bewirkt also Verdrängung und Verdichtung des Gases. Die Elemente der Dipolfläche machen 
einfache Bewegungen in der Normalen. Diese haben zusätzlich Umströmungen der Flächen- 


4) Die Druckmeßmethoden sind im Anhang II beschrieben. Sie wurden von W. v. Wittern entwickelt, 
dem ich darüber hinaus für interessante Diskussionen zu danken habe. 
5) Rayleigh, Theory of sound II $ 325. 
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nnd zur Folge. Die von den Flügelelementen durchlaufenen Bahnen kann man sich mit 
F punktförmigen Quellen besetzt denken. Dabei entsprechen: 
Die Strahleramplituden der Profildicke und dem Anströmwinkel. 
| Die Frequenz der Zahl der Tragflügel und ihrer Geschwindigkeit. 
Der Phasenwinkel dem Ort des Tragflügels auf seiner Bahn. 
Das Fourier-Spektrum der Form des Flügelquerschnittes und seiner Zirkulationsverteilung. 


Das Potential der Strahlerflächen 

Für die Beschreibung dieser Vorgänge setzen wir voraus, daß die dem Gas erteilten Ge- 
schwindigkeiten von einem Potential abzuleiten sind. Sie sind als die Normalgeschwindigkeiten 
der Strahlerflächen klein gegen dje Schallgeschwindigkeit, selbst wenn die Spitzengeschwindigkeit 
des Tragflügels diese überschreitet. Wir rechnen im folgenden daher mit der linearisierten 
Wellengleichung. Die Gültigkeit dieser Betrachtungen kann erst durch das Auftreten von Ver- 
dichtungsstößen eingeschränkt werden. Dies zeigen auch Untersuchungen über das Schallfeld 
von Luftschrauben!!). Dabei ergab sich bis zu Umfangsgeschwindigkeiten von der Größe der 
Schallgeschwindigkeit eine gute Übereinstimmung der Meßergebnisse mit den auf der oben dar- 
gelegten Grundlage berechneten Werten für den Strahlungsverlust. 

Den Zusammenhang zwischen den Kräften, die auf ein Gasteilchen wirken, stellt die Euler- 
sche Gleichung her. Wir führen die vom Tragflügel verursachten Störkräfte als eine äußere 
Kraft K ein und erhalten, wenn p den Druck auf die Oberfläche eines Teilchens, » seine Ge- 

) schwindigkeit und o seine Dichte bedeuten 


9 
0 & + (vgrad) I BAU D RE N en A (1). 
Die Bedingung der Kontinuität ergibt: 
00 EEE 
4 iv) PER NIE I NE EIER: (2). 


Die Teilchengeschwindigkeit kann von einem Skalar p, dem Geschwindigkeitspotential abgeleitet 
werden (v= — grad p). 
Mit der Vektorbeziehung 


(v grad) v— 5 grad v? — 5 u Be Er Ta (3), 


worin das letzte Glied wegen der Wirbelfreiheit verschwindet, erhält man aus (1) und (3) 
P 
D) 1 ... fd 
a er gran Be (de 
oi 2 0 [ 
% ® 
Auch das quadratische Glied kann wegfallen, da die Teilchengeschwindigkeiten klein gegen die 


Schallgeschwindigkeit sein sollen. Nach Prandtl®) leiten wir die Beschleunigung B = en von 
einem Skalar Y, dem Beschleunigungspotential ab: eg 


Bert) Sea div Be — AP N. Men ns (5). 


Wir können dann für (4) schreiben 


u SE ae 


— — 2 grad p-+ grad je 2 ad On si (6). 


Sind die Geschwindigkeiten klein, so ist der Druck der Dichte proportional. Aus der Adiabaten- 
_ beziehung erhält man 


OPEN DE 0 DSH it. 
ne EEE mt  o=n(l+e) gilt: 
P 
[*?=e fiezemar ee. ER OIETREREENL N (7). 
ge 


Dr L.Prandtl, Theorie des Flugzeugtragflügels im zusammendrückbaren Aieslunt: Luftfahrtforschung 
13, 1936, Lieferung 10. 


e a4 fi } Fi ” £; / . #4 Zi 
24 Ernsthausen, Derrotierende Tragflügel äls Strahlungsproblem ya sm ee 


Durch Integration und mit Beachtung von (7) folgt aus (6) 


— 4804 W-+ oonst = I a ns (8). 
Hieraus durch Differentiation und wegen (2) und (7) 
18 Ir 
rer ne NrTrru TE ne Nahe Le ET (9). 


Da die rotierende Bewegung gleichbedeutend mit einer harmonischen Erregung ist, ergibt sich 


10% [n) ’ 
Ip = — — —; a a a Fr BR Fr ; 
Tr rTEG k > (10) 
(Die Zeitwerte sind dabei wie die Amplituden geschrieben: p = ge‘®t,) Die Lösung der Gl. (9) 
ist nach Rayleigh’) | 
LEE Ri, Fett. 
RR: An/J@a r 


—ikr 


(r—= Abstand zwischen Quelle und Feidpunkt), worin E eine Lösung .der Gleichung 
Ay + k?p= 0 für Kugelkoordinaten darstellt. 

Zunächst betrachten wir die Kräfte, welche auf die zirkulationsfreie Bewegung eines Körpers 
zurückwirken. Der Widerstand, den der Gasraum dieser Bewegung entgegensetzt, ist bestimmt 
durch den auf beiden Seiten der Strahlerfläche wirkenden Druck. Dem symmetrisch umströmten 
Körper entspricht dabei eine gegenphasige Bewegung beider Flächenseiten (Bild 1). 

Das Potential einer derartig erregten Fläche läßt sich, wie Rayleigh®) gezeigt hat, aus (11) 
unter der Annahme gewinnen, daß die Funktion 7 (x, y,2,t) außerhalb eines beliebig kleinen 
Gebietes um den Angriffspunkt der Kraft verschwindet. Dann geht des Volumintegral in ein 
Flächenintegral über. 


Wir integrieren (9) über dV unter der Voraussetzung, daß Y so wächst wie V abnimmt. 


02 
Damit verschwindet f SB dV. Wenn grad gleichfalls so ansteigt wie » abnimmt, ergibt 
E n - > 


4 i 
sich er 
mr, min EraT 
Damit erhalten wir | 
frifen 
und hieraus mit (11) | 
—ikr 3 
9° = — — | grad ®8° Re ER Ur RaeaN 
n e s 
wenn wir grad g durch die Amplitude der Flächengeschwindigkeit UP = — grad ©° ersetzen. 


Hier soll der Zeiger 0 für die Nullstrahlerfläche (wie später der Zeiger 1 für die Dipolfläche) gelten. 


Da das Integralin (13) symmetrisch zu Fist, stellt ®, das auf beiden Seiten gleich verteilte 
Flächenpotentialdar. Es setzt sich stetig durch die Fläche fort, während hierbei die Gschwindig- 
keit um die Phase x springt. Das Strahlergebilde entspricht einer Doppelmembran, die gegen- 
phasige Bewegungen ausführt, also mit Druck belastet ist. Ihr Ersatzbild ist eine mit bz va 
Quellen belegte Fläche (Nullstrahlerfläche). 2 

Die Wirkung eines Tragflügels besteht aber in der Entwicklung des Auftriebes. Diese‘ Ann 
Verteilung von „konzentrierten Kräften“ über seine Fläche stellt den erregenden Vorgang der Ar 
Dipolfläche dar, durch welchen das Gas senkrecht zur Strahlerfläche beschleunigt wird. Den der 
Beschleunigung zugeordneten Druck P an jedem Flächenpunkt bestimmen wir aus der-Poten- AR P 
tialdifferenz ®! von Ober- und Unterseite des angestellten Flügels, die gleich der a sein Br 
muß. Hierzu ziehen wir Gleichung (11) heran, in welcher die Wirkung 3 Kraft be 
das Beschleunigungspotential 7 ausgedrückt ist. Dieses muß, wenn es: sc) ‚wirbel- 
freier Gasbewegungen geeignet sein soll, ebenso wie 4 Gern 


?) Theory of Sound II, $ 277. ve 
8) Theory of Sound II, $ 278. Ban 
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gleichung (4Y-+%k2Y = 0) befriedigen. Wir erhalten, da nach (5) Yund somit auch B periodisch 
in 9 ist 


1 Ik BEN .. 
(11) schreibt sich dann FAR 
| | en Ik BE uv 15 
i PR IV EN EI Ne TEILTE 
v 
F woraus man wegen der Vektorbeziehung) 
} fpdiv BV=— ((Bgradp)dV ........22..(16) 
1 14 Bahn 
e 
das Potential einer periodisch auf das Gas wirkenden, in einem kleinen Gebiet konzentrierten 
h Kraft erhält: 
| aa |» RE ar 17 
= — — oO 

Yp Ser grad — an la ALTE, 
Vv 
| Wenn Yaußerhalb eines endlichen, beliebig kleinen Gebietes von der Breite dn verschwindend 
kleine Werte annimmt und Bin Richtung der Normalen wirkt, dann ist wegen (5) Y — — il Ban. 
| Und wenn wir n so abnehmen und B so wachsen lassen, daß ER Bdn = — Y wird, geht das 


Raumintegral (17)in ein Flächenintegral über. Wir ersetzen daher T Bav durch BdndF = — YdF 
und erhalten 


i % —ikr 
Pa Yan DE Ba RL EEE LÄDT 
% Wir führen. 
‚ vaDdlı. pP 
E re ON 
} 


ein. Damit ergibt sich das Potential einer mit konzentrierten Kräften erregten Fläche. 
l 


) te \ 
Seren re Ar 22.2.2... (0) 


Aus (20) eixtnehmen wir, daß P bzw. Y auf der einen Seite positiv, auf der andern jedoch negativ 
ist, also um die Phase x springt, während seine Ableitung, welche der auf das Gas wirkenden Be- 


> schleunigung gleich ist, auf beiden Flächenseiten gleiche Richtung und Größe hat. Das Potential 
- entsteht hier durch gleichphasige Bewegungen beider Flächenseiten, muß also durch eine Dipol- 
E; belegung dargestellt werden 1°). Da Y bzw. P nach Voraussetzung auf der Dipolfläche vorgegeben 
b ist und außerhalb dieser verschwinden soll, kann mit Hilfe von (20) die resultierende Wirkung 
R im Außenraum ermittelt werden. 


Die Geschwindigkeit der Strahlerflächen 


F: Wie bereits bemerkt, werden Bewegungs- und Zirkulationsverlauf in einem Punkt der 
4 Strahlerfläche bei periodischer Erregung zweckmäßig durch ihr Fourier-Spektrum ersetzt. Die 
Fourier-Koeffizienten sind aus den geometrischen Flügeleigenschaften (Flügelform- und Aul- 
 triebsverteilung über die Flügeltiefe) zu ermitteln, ‘Wir bezeichnen sie für Nullstrahler- und 
 Dipolfläche mit fm und fn- 

2 #, "Die harmonischen Komponenten der zu behandelnden Vorgänge kennzeichnen wir mit dem 
Zeiger m. Wenn wir die Ordnungszahl der Fourierglieder » und die Flügelzahl z einführen, ist 
on vz. Die Flügelzahl ist — gleicher Winkelabstand der Flügel vorausgesetzt — ein ganz- 


Pr zahliger Frequenzfaktr wie die Zahl v. Wir schreiben jetzt für k in (10): k= m — } 
air, ee ) 
Een, 

Fi ese Beziehung gilt nur, wenn der Integrand für sehr große r stärker als 1/r® abnimmt, was durch die 


einer „konzentrierten“ Kraft gegeben ist. 
at: m b (Dynamical Theory of Sound, $ 78) hat gefunden, daB, das Potential einer „konzentrierten 
© ee Irrelmele äquivalent ist. 
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Die Geschwindigkeit der Nullstrahlerfläche ist durch die Dicke d® des Tragflügelprofils 
gegeben. Wir erhalten für eine Harmonische, wenn wir für die Bewegung d,,(t) = d„e'**t und 


2 1 
daher für U), = 7 


d ansetzen, h r 
ds 
U] =ikeT; a OF TR VO DRS Eee (21). 

Das Nullstrahlerpotential (13) besteht aus einer Fourierreihe, deren Anteile durch die Form 
des Flügelkörpers bestimmt sind. 

Die Verteilung der Normalgeschwindigkeit U!(kr) auf der Dipolfläche ist, wenn J'oder P 
gegeben sind, nicht allein durch ihre radiale Abhängigkeit, sondern auch noch durch das Strah- 
lungsfeld bestimmt. Deshalb setzt sie sich aus einem inkompressiblen und einem kompressiblen 
Teil zusammen. Entsprechend dem Vorgehen bei der Nullstrahlerfläche wird der inkompressible 
Teil von U! berechnet. 

Die Potentialdifferenz der Flügelober- und -unterseite ist bei z Flügeln 


Bean Ns a Se 


wo mit /’ die Zirkulation der inkompressiblen Strömung bei vorgegebenem Anströmwinkel a* 
am Element eines unendlich langen Tragflügels bezeichnet ist. 

Der Impuls, den der Tragflügel bei jedem Umlauf dem Gas erteilt, hängt von Zirkulation, 
Blattzahlund Drehzahlab. Erist gleich dem Stoßdruck ®t. Wir schreiben@aher für die zeitliche 
Änderung in einem festen Raumpunkt ®!=zJ'ei®t. Die Kraftamplitude je Flächeneinheit ist | 
dann unter Verwendung der Kutta-Joukowskischen Auftriebsformel 


0 d%y . 
P=- 601 emiat — 
a R e 


Sie beschreibt den Druck, der im inkompressiblen Medium auf das erregende Element allein 
infolge der Massenbewegung wirken würde, und damit die durch die geometrischen Eigenschaften 
erzwungene Erregungsgröße. Wenn wir I'„= I'- fm setzen, können wir für eine Teilschwingung 
EN R Ale RL WL & 24 
a ek Im; De 
schreiben. 


Wir führen im folgenden für die Zirkulation die Beziehung 
T= 3 bo NE N Mt (25) 
mit c,„ als mittlerem Auftriebsbeiwert, b als Tragflügeltiefe, v, als Umfangsgeschwindigkeit 


und b„ gleich 5 - fm ein. Dabei nehmen wir Z’oder P längs des Flügels zunächst als unveränderlich 
an und setzen c, gleich dem Mittelwert der hierdurch bestimmten radialen Verteilung. 


hl i K- 
Aus der Kraftamplitude erhält man wegen P:F= me die „inkompressible“ Ge 8 
schwindigkeit U = | Pdt. Mit F=nR? und der mitschwingenden Dipolmasse für inkom- 
pressible Bewegungen M—-* 5 rs R® nach Fußnote 14, sowie mit (24) wird ; > 
BR | 
BR en da in wären ET Fa ae ae: un 5 . (26) 3 
und daraus mit (25) Br 
Ba 3% 042 2 

Un= ike— b ne =ikedi; da Falle en 


(d,„ sei die Auslenkung der Elementarstrahler). 
Die Beschleunigung ergibt sich dann zu 
DULLHR 
Bi = — k2c 2 F by en 


m: 


et Kir di 


Die Rückwirkung des Schallteldes 
Hierunter verstehen wir die Kräfte, die das Schallfeld auf di 


m 
“ 
% 
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n 
. 


Ba 


Die an das Fernfeld gebundenen Strahlungsverluste wurden in verschiedenen Arbeiten 
behandelt!!). Sie sind identisch mit zusätzlichen Bewegungswiderständen. Diese stellen jedoch 
nur einen Teil des kompressiblen Einflusses dar. Zur Ermittlung der den Auftrieb bestimmenden 
Massenkräfte muß das Nahfeld in die Betrachtung einbezogen werden. 

Wir beschränken uns darauf, daß der Tragflügel auf der Stelle rotiert, womit sich das 
Potential 9 der Fläche harmonisch verhält. 


a) Die Nullstrahlerfläche 


Zunächst berechnen wir die Wirkung, für welche. nach (13) die Normalgeschwindigkeit U, 
zuständig ist. Hierzu haben wir über die Druckverteilung zu summieren. Es ist 


K= [pdF=o/[gdF = ikco [pdF RE Be U 0) 
F F F 


Aus (13) und (29) erhalten wir unter der vereinfachenden Voraussetzung, daß die Erregung längs 
des Tragflügels konstant ist — dies bedeutet konstante Dicke —, für eine harmonische Kraft- 
komponente 


ık U —ikr—imx 
goal dar KR a in. 
4 2 N 


| Im Zähler tritt ein Faktor 2 auf, da die Strahlerfläche bei der Auslenkung um 2 [vgl. (24)] doppelt 
zählt. "x 
Das Doppelintegral behandeln wir an Hand von Bild 4. Zu der räumlichen Phase kr in (30) 
tritt die durch die rotierende Bewegung gegebene zeitlich-örtliche Phase ma. Zunächst wird 
das Potential in ? bestimmt, herrührend von der Kreisfläche zr?. Damit ist auch der Wert 
des Integrals für die Ringfläche 2 ar,dr, bekannt, da es gleich in jedem 
ihrer Punkte ist. Wir erhalten schließlich das gesuchte Integral, indem dr; 
wir das Ringintegral von O0 bis $ nehmen, wenn R die Länge des Trag- (A 
flügels ist. So ergibt sich für die auf die Nullstrahlerfläche wirkende Kraft EN 
der m-ten Komponente 


Ka a RUE > N RR 6) Biker se Koordinaten der 


Nullstrahlerfläche. 


? acda a R 


In (31) ist 
1 ikr—ima 
GH) al) — . Ya dr; dr„da 
R !aq PA 


“ 1 eikr—ima ulguge (32). 
alte] nf I aalantarı 
Ö Ö Ö : 
"Worin nach Bild 4 
> ae (88): 


Dieses Integral läßt sich so weit behandeln, daß es berechenbar wird!). Schreibt man nach 
Weyrich für den Integranden 


E: \oikr 
A EA 
| r R YA — 2 
| so ergibt sich durch die Auftrennung des Integrationsweges & in der YA? — k2-Ebene nach 


Busemann, Zur Gasdynamik des drehenden Schaufelsterns. ZAMM, Bd.18, Heft I, 1938. 
Das räumliche Problem mit Zirkulation: { 
Ernsthausen, Der Einfluß aerodynamischer Eigenschaften auf Schallfeld und Strahlungs- 
leistung einer Luftschraube. Luftfahrtforschung 18, Lieferung 8, 19, Bd.9, H.]. 
Auch die Arbeit von 
Gutin, Über das Schallfeld einer rotierenden Luftschraube. Physikalische Zeitschrift d. Sowjet- 
a Union, Bd.9, H.1, 1936. 
- bezieht sich nur auf die Strahlungsverluste. Sie werden mit Hilfe der als bekannt vorausgesetzten Schub- 
und Drehmomentwerte berechnet. Die Luftschraubenfläche wird durch einen Kıreisring angenähert. Der 
Verfasser hat dabei das gesamte Drehmoment, dem ja in erster Linie die kinetische Energie der Rotation des 
- Schraubenstromes entspricht, für die Entstehung der Schallenergie verantwortlich gemacht. Da es sich aber 
um Luftbewegungen handelt,die vorwiegend tangential zur Strahlerfläche verlaufen, kommen sie für die 
ur Roma Dirht in Betracht. Der durch diesen Drehmomentsbeitrag gegebene Widerstand unterliegt deshalb 
auch nicht dem Einfluß der Kompressibilität in dem hier ausgesprochenen Sinne. 
" 12) Hierfür bin ich R. Weyrich zu großem Dank verpflichtet. 


h 11) Das ebene Problem ohne Zirkulation: 
“ 
14 


Real- und Imaginärteil unter ee umfangreicher Zwischenrechnungen (Rekursion) RE E 


|, [o.am 2 Tre + man = Ei lee 
Durch die Substitution A=kbofu und A=kcosv wird 
G,—n{oy(kR)—iS%(kR)} mit 


= Som (kRCofu)du und Fr aa 
\ Me 
8,70, (kBeosv)dv ; 
6 
4 
r —ikr 
« worin das Minus-Zeichen wegen : : steht. Die Funktionen E 
; AR al SE | 
. Du(AR)= 378 =|, j BREUER Der BERN tus (36) | 
R 
berechnen sich für m= 2,4,6,8 und mit @=AR zu 

is 4 ib BE 
1 8 24 5 Ei 
= J,(2) — Be: "| ra)+ | 
2 12 60 " \ s (37). Br 
DB, — a () Il) +3 J3(2)+5J;(2)] — Bi: +8} R 
Aus (31) und (35) ergibt sich die Amplitude der vom Schallfeld auf seinen Erreger ausgeübten 2 
Kraft ir 
"% 
— Kun = UncoFkRSn(kR)+ DEAD BESSER). 2 =. 4.88). sr 
Das erste Glied ist der Geschwindigkeit U}, proportional, ist also eine Verlustkraft. Das zweite $ 
Glied ist der Beschleunigung By = ikcUn, BIOBOrDONGE und daher eine Massenkraft. Hieraus : 
bestimmen wir den Strahlungswiderstand “a 2 
Ru = cgFkRSy (BB): 22 Ken Ku E 

und den Massenwiderstand ET Ren a Ka 

= 00 Fk REHTKR) 3 er Pr 0 . (40). 


Sm(kR) bezeichnen wir im’ ERNIN als Strahlungsfunktion und &,(k.R) als Massenfunktion. 3 
Sie wurden an Hand von (35) und (37) für m = 2, 4, 6, 8 berechnet (Bild 5aundöb). Diedamitt 
gebildeten Widerstandsfunktionen k RS}, und k R ©, nach (39) und (40) enthält Bild 5eunddd. 
(Die Kurven für verschiedene Werte von m auf ‚der Abszisse gegeneinander versetzt, 


da der Maßstab mit m multipliziert ist: kR = mn ®  Beider Bildung des resultierenden Wertes Ei 


muß daher der Maßstab jeder Kurve durch das Zuger m dividiert er ‚Es rücken dann 2 > 


alle Kurven der Kurve m—= 2 AU 4 NR, AR 


Me b) Ein mechanisches Ersatzbild ” = 
Die Gleichung (38) ist von der Form - TER 

— Kn=U Um dm 3a ER Be 

worin die komplexe Impedanz TH 
| Bet RE en Fire) 
ist. 
Um ein einfacheres Bild von den Strahlungsvorgängen zı zu 
Impedanz Z}, durch diejenige eines mechanischen ’Schwingers mit . 


widerstand darzustellen. Hier kommen uns die Überlegungen v 
über Se zu Hilfe. In der  ; Arbe 
EURER Ye 


N a 2; HE R a v un as SE 
ENTE, ei Re Y R ER AN? 


‚ fr a an ar \ s 
Ri Be Ele an ehr. 191 Ernsthausen, Der rotierende Tragflügel als Strahlungsproblem 29 
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strahler als auch für den Dipol mechanische Ersatzbilder abgeleitet, deren Impedanzen durch 
die Wahl geeigneter Konstanten mit den potentialtheoretisch ermittelten übereinstimmen. Wir 
verallgemeinern diese Gedankengänge, indem wir das mechanische Modell des strahlenden Dipols 


ee 
=| 


2 


+ 8 


h> > Bild 5d. Der Verlauf des Massenwiderstandes der 
E et a  ertluke ne re Nullstrahlerfläche 
F als typisch für den Strahlungsmechanismus überhaupt ansehen, da hiermit indirekt die normalen 
h und die tangentialen Gasbewegungen erfaßt werden. 
Er: Die Impedanz eines solchen, in Bild 6 gezeigten Schwingungsmodells läßt sich dann in 
ä folgender allgemeiner Form schreiben: 
Ye a 1 1 
et! 0m NER.) 
t i - 3 o* —i NEN & 
% 


Kaßel 4, 0*,o die dimensionslosen Koeffizienten von Masse, Widerstand und 
Steife bedeuten und x der Frequenzkoordinate kR gleich ist. Durch Tren- 
nung von Real- und Imaginärteil erhält man 


2 
a 
Eh) e 2,72 0* RO x 
eh 8 zu 0 dh 21] ra EEE ERIRER N EERBEE (44) Bild 6, Das mechanische 
nn 2 + ze Br “.  Brsatzbild eines 
» , c *2 2 ,,2 a Strahlers 
Merian zue € TE 72 2n0 


Nee wieder als Verlustwiderstand und der Imaginärteil als Massenwiderstand 


Str: | 
m) 2.B. ergibt sich aus (44) mit den Konstanten u= 1, e*=1, 5 —=0 E DT pe und mit 
Rn e er 0* = er die Dipolimpedanz. 


ur Ka N > air u 
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sein unbekannter Massenwiderstand näherungsweise ermittelt werden; denn beide sind durch 
4, 0*,o völlig bestimmt, da hierdurch die geometrischen Dimensionen festliegen. 


Wir prüfen diese Aussage an Hand der für die Nullstrahlerfläche streng ermittelten Funk- 
tionen von Bild 5. Dazu nähern wir z.B. 8) möglichst gut an!?) und berechnen danach die 
Funktion ©. Das Ergebnis zeigt Bild7, wo die genauen Kurven ausgezogen und die Näherungs- 
kurven punktiert sind. 

Da sich Strahlungswiderstände wegen des unendlichen Aufpunktabstandes immer viel 
leichter berechnen lassen als die Massenwiderstände, kann man mit Hilfe dieses Verfahrens auch 
bei komplizierten Strahler- 
formen die Vorgänge an der 
Strahlerfläche untersuchen, 
vorausgesetzt, daß sich der 
Verlauf des Strahlungswider- 
standes durch den ersten 
Term von (44) gut genug 
erfassen läßt. Schwierigkei- 
ten können allerdings bei 
höheren Harmonischen auf- 


g 2 3 4 5 6 7 8 9 ”0 " 72 ZRRD? : 
Bild 7a, Die mit Hilfe des mechanischen Ersatzbildes angenäherte Strahlungs- treten, wenn der einfache 
funktion Sy (gestrichelt) im Vergleich mit dem Ergebnis der genauen Rechnung Bau der Näherungsfunktion 


den zu fordernden Ansprü- 
chen nicht mehr genügt. 
Geht man in (44) von kon- 
stanten Koeffizienten u, o*, 
co zu Funktionen wu(e), 
o*(z), o(x) über, so lassen 
sich durch den Realteil auch 
Strahlungswiderstände mit 
verwickelterem Verlauf an- 
nähern. Dies entspricht der 


Einführung frequenzabhän- 
Bild 7b. Die en Hand des mechanischen Ersatzbildes näherungsweise bestimmte giger Elemente in das 

Massenfunktion ©g (gestrichelt) im Vergleich mit dem Ergebnis der genauen Rechnung . . 
Schwingungsmodell (Bild 6). 


Die Gasbewegung kann nach diesen Darlegungen auch so gesehen werden: Werden dem 
Gas (wie Punkt A)nur kleine Geschwindigkeiten erteilt, so bleibt die Feder-starr, es findet eine 
reine Massenverdrängung wie bei der inkompressiblen Bewegung statt. Sie ist in Phase mit der 
aufgedrückten Bewegung, der Widerstand arbeitet nicht, es entstehen keine Strahlungsverluste. 
Erhöhen wir die Geschwindigkeit — in unserem Modell die Frequenz — so wirkt die Feder. Die 
Amplituden übersteigen die des Punktes A bis zu einem Maximum, sie sind dann um 90° phasen- 
verschoben. Bei weiterer Erhöhung der Frequenz nehmen die Massenausschläge ab, da der Feder- 
widerstand kleiner und der Massenwiderstand größer wird. Schließlich führt die Masse keine 
Bewegung mehr aus, die zugeführte Energie wird im Strahlungswiderstand verbraucht. 


c) Die Dipolfläche \ 

Hier müssen wir wie bei der Nullstrahlerfläche ein Zusammenwirken aller Flächenelemente 
zu einem Interferenzfeld annehmen. Dies bedeutet, daß der einzelne Punkt über das gemeinsame 
Feld auf die andern zurückwirkt. Da die Vorgänge im Feld eng mit denen an der Strahlerfläche 
gekoppelt sind, erzeugt jeder Elementardipol einen anderen Schallwiderstand für die übrigen 
Strahler der gleichen Fläche. Sie können nicht ungestört strahlen, obwohl sie nicht direkt von 
fremden Wellen getroffen werden. An der Strahlerfläche muß aber die Schnelle des Gases in die 
Membrangeschwindigkeit übergehen. Wird der erregende Druck erzwungen, so stellt sich eine 
vom Gesamtield abhängige Geschwindigkeit ein, wird diese erzwungen, so ändert sich der Druck 
an der Fläche. Die Berechnung der Schallfeldkräfte an der Dipolfläche mit Hilfe von (20) ver- 
langt deshalb einen bedeutenden Aufwand. 

Die gute Übereinstimmung in Bild 7 berechtigt uns indeß, auch diese Vorgänge mit Hilfe 
des Ersatzbildes in Bild 6 zu beschreiben. Wir wollen deshalb mit (20) nur die Strahlungswider- 
stände und über diese erst die Massenwiderstände ermitteln. | 

Wie angegeben, soll die Strahlungsleistung zunächst für eine längs des Scheibenradius 
konstante J-Verteilung berechnet werden. Da jedoch im Strahlungsfeld höchstens der Energie- er 
betrag für ebene Wellen untergebracht werden kann, läßt sich ein endlicher Zirkulationswert am ° 


'5) Dies gelingt mit den Koeffizienten a, = 0,113 o% = 0,4219 = 2,0896. N EN 2J 
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Tragflügelende physikalisch nicht aufrecht erhalten. Der Strahlungswiderstand kann aber aus 


der Strahlungsleistung immer durch die Bedingung gefunden werden, daß er für AR>w in 
einen konstanten Wert, dem der hier der ebenen Welle zukommt, übergehen muß. Auf das 
gleiche Ergebnis führen strömungsphysikalische Überlegungen wenn man bei der Berechnung der 
Strahlungsleistung die aus der a eiblen Tragflügeltheorie verfügbare radiale Auftriebs- 
funktion einführt. 


Ein Teil der Tragflügelenergie findet sich in dem vom Tragflügelende abgehenden Wirbel 
wieder. Da hierdurch eine Wechselkraft weder auf den Körper noch auf die Strahlerfläche 
entsteht, erweist sich auch der ‚‚induzierte Widerstand‘‘, der diese konstante Rückwirkung dar- 
ehr hier als eine inkompressible Größe. 

- Mit (20) gewinnt man die Strahlungsleistung der Dipolfläche, indem man die Strahlungs- 
intensität auf einer in sehr großem Abstand die Strahlerfläche konzentrisch einhüllenden Kugel 
summiert. Wir führen sphärische Polarkoordinaten an Hand des Bildes 8 ein, und erhalten das 


0) 0 
Potential für sehr große Werte von r und wegen an an ß 
urn ei (kr—mo) F 
1 225 r 
oh = aß Deal rein 232,45). 


Der mit (24) gegebene Druck P,„, nach Voraussetzung über die Fläche 
konstant, steht vor dem Integral. Zu der räumlichen Phase kr tritt die Um- 
laufphase ma. Die zum Aufpunkt führenden Strahler r und r’ sind wegen 


Bild 8. 
Zur Berechnung des 


- des unendlichen Abstandes parallel. Man entnimmt dem Bild die geometrische Strahlungspotentials 


der Dipoliläche 


Beziehung: E 
2 Dee rn Bose . (46). 
Alle konstanten Phasen seien ebenso wie das Zeitglied weggelassen, so daß nur die Änderung 
der Amplitude bleibt. 
Dann ist 
Fe (kr, sin Bcos%— ma) 
= DR: s TUI lee (47) 
oder 
N PS ee 
7 ee EOS era ran D)aran u en a a ers AB 
i - 0 
Mit z= kr,sin ß erhält man 
WAT, 1 
| = (R) Ba a eu (49), 
worin 
Z 
EEE a | male: ENTE RO or RER (50). 


Das Integral kann nach (37) berechnet werden. Der effektive Schalldruck ergibt sich durch 
Ableitung nach der Zeit, 


n UHR 1 t 
et > U k R a IE AA ART & 
Beil 45 Tor %m(P) (51) 
Die Strahlungsleistung 
WW ’ 2 sr’? 1 
NL, = FF I a te. (52) 
führen wir nach den vorausgegangenen Überlegungen in 
U2 
Ni, = > a ur Bye ss ER 3 (53) 


üben, wo für ebene Wellen U == ink, 


m folgt für den Strahlungswiderstand entsprechend (39) 
j NS coF ee N ar a ern (54) 
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mit den Strahlungsiunktionen 
v1 n® Ü 12 . | a 
Sm=zkR Yu (B) sin BABY el. N ever 3 

M) i 


Die Funktionen $4,(k. R) wurden mit (50) bzw. 7) für m—2,4,6,8 berechnet und in Bild 9a 
aufgezeichnet. 

x Die mit $% verknüpften Miissenfunk onen ©, bestimmen wir nach Gl. (44). Dazu nähern 
wir die'mit (54) aus dem Fernfeld berechneten und in Bild 9a ausgezogenen Strahlungsfunk- 


Bild9a. Die Strahlungsfunktionen der Dipolfläche Bild 9b, Der Verlauf des Strahlungswiderstandes der Dipolfläche 


tionen S7 mittels der Koeffizienten u, 0%, 0 
möglichst gut im Gesamtverlauf an. Dies er- 
reicht man mit konstanten Werten für „und 
O7, und den Funktionen 07,(x). Von den letz- 
teren muß man voraussetzen, daß sie grund- 
sätzlich ähnlich den Strahlungswiderständen 
in. Bild 9b verlaufen. Denn der bei N! 
Frequenzen dominierende Einfluß von o* 
nimmt nach tiefen Frequenzen hin stark ab, 
da die Dipolfläche hier ein schlechter Strahler 
ist. Für diein Bild Ya gestrichelten Näherungs- 
kurven ergeben sich die in der Fußnote 1%) 
angeschriebenen w„- und o„-Werte und die  -» 
om -Funktionen in Bild 10. 


Damit schreibt sich die Impedanz der Dipolfläche er. 
= 00 FRISCHE} Dr a Ina 


! 


2 4 6 8 no 12 % 16 18 20 
Bild 10. Die Funktionen en (k R) 


Mit : 
©O, l + 

Bmen ER)— mund E00, EN 

A rn ER. = EURE >. HRS Kö 

om (kR)-+ Re Re ui + AB — 2 Um Im N + NUR» h z BEN 

2 k j ee IV I, Be \ 

erhält man analog (40) aus (44) für die Massenwiderstände der . N « $ HR RL Be 
M„=ceFkRCL(ER) . En >: 


Sie nehmen mit kR nicht so schnell ab wie bei der Nullstrahlerfläche, da die RER. Gas- 
reine durch die RR EN ‚der Flächenberandung auch Aa bei Fe 


Dom Loorvi RAS, +Hie) 


16) m 2 4 6 8 


Um 0.066 0.0346 0.0253 0.02 IEER 

m: 20 3.526 4.976 6.4 Su 
Au 

a ER 
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Nach diesen Betrachtungen lassen sich die wesentlichen Eigenschaften selbst eines so unüber- 
sichtlichen Strahlergebildes wie es der umlaufende Tragflügel ist, schon an dem Schwingungs- 
modell von Bild 6 wiederfinden. Der Einfluß der Kompressibilität ist dann durch die Wirkung 
einer Kombination von Feder- und Verlust- 
widerstand darstellbar, deren Größen sowohl 
vom Verhalten des Gases als auch vom Cha- 
rakter des Strahlers abhängen. 


kRTm 


0 2 4 6 8 70 12 7 76 7 20 
Bild 11a. Die Massenfunktionen der Dipolfläche Bild 11b. Der Verlauf des Massenwiderstandes der Dipolfläche 


Die Kräfte am Tragflügel 
In die Kraftgleichungen (38) und (59) setzen wir die Geschwindigkeiten der Strahlerflächen 
nach (21) und (27) ein. Für die einzelnen Komponenten erhält man 


d? 
a =—ikgFZKkRSnm 
0 


Ku Bel@noRsch, 
= —ikeoFdnkRSm 
Ky= KadnnoR? Sm 


Die Strahlungskräfte K% und K%sind durch Geschwindigkeit und Strahlungswiderstand 
und die Massenkräfte Kr und Kir durch die Beschleunigung der mitschwingenden Gasmasse 
gegeben. 

Diese Kräfte sind entsprechend ihrer Entstehungsweise als zeitliche und räumliche Funk- 
tionen auf der Strahlerfläche verteilt. Sie laufen, da die Frequenz der festen Elementarstrahler 
gleich der Winkelgeschwindigkeit des Tragflügels ist, synchron mit diesem um. Der durch die 
Tragflügelbewegung erzwungene Erregungsvorgang hat einen Spitzenwert am Tragflügel selbst. 
Daher müssen die Amplituden der Kraftverteilung auf der Fläche gleich den konstanten, am 
Flügel angreifenden Kräften sein, die durch das von allen Flächenpunkten gemeinsam erzeugte 
Schallfeld und damit durch die kompressiblen Eigenschaften des Gases gegeben sind. 

Die Widerstandskräfte K% und K% entstehen durch die Abstrahlung akustischer Energie; 
denn jede rotierende Bewegung führt im kompressiblen Gas zur Ausbreitung akustischer Wellen. 

Aus dem Betrag dieser Kräfte läßt sich je ein Widerstandsbeiwert für die Strahlungs- 
verluste definieren: 


" 


Sie verlaufen gleich den Strahlungsfunktionen 8%,(k R) und S4(kR) (Bild 5a und 9a). 
Die Massenkräfte Ky zu beiden Seiten der Strahlerebene heben sich am Flügel auf; denn 


_ die Beschleunigung gleicher Gasmassen nach entgegengesetzten Richtungen, wie dies die symme- 
 trische Anströmung bewirkt, kann keine Rückwirkung liefern. Ihre physikalische Bedeutung ist 


aus dem Schwingungsmodell zu ersehen. Der hierdurch bedingte Massenwiderstand ist bestim- 
mend für den Strahlungsverlust. Die dazugehörigen Massenfunktionen ©, mit deren Hilfe die 
Berechtigung des Ersatzmodelles Bild 6 nachgeprüft wurde, enthält Bild 5b. 

Bei der Besprechung der Massenkraft Kr halten wir uns an Bild 12. Wir gehen von der 


örtlichen Zirkulationsverteilung um den Tragflügel aus, welche durch das Potential ©! oder 
durch die Auslenkung d;, der Elementardipole darstellbar ist. Da die Beschleunigung um die 
Phase x gegenüber der Auslenkung gedreht ist, verläuft d„, gleich —Bt. 


E 
re j 3 
f ey > % » 4 
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Wie die Strahlungskräfte haben die Massenkräfte ihren Höchstwert immer am Ort der 
größten Auslenkung, d.h. am Tragflügel. Dies gilt auch für alle harmonischen Teilwellen der 
Zirkulation. Diese müssen am Tragflügel gleichgerichtete Kraftkomponenten ergeben, da ihre 

5 Strömungen ihn als Kern und Kraftangriffs- 

MR punkt immer im gleichen Sinn umschließen, 

x ET Die Reaktion besteht in der Beschleunigung 
der Gasmassen im Raum zwischen den Flügeln. 

Dies zeigt der unterste Kurvenzug des Bildes 12. 

Die zeitlich-örtliche Kraftverteilung ist hier bei- 


v— IR. spielsweise als Resultierende von zehn gleich 
großen Teilwellen entstanden, welche am Trag- 
SOFNIIF< flügelort ohne Phasenunterschied, wie im oberen 
Teilbild angedeutet, einsetzen. Die kompressible 
Massenkraft am Tragflügel, der Schub der ge- 
samten Anordnung, ist daher unmittelbar gleich 
der algebraischen Summe der harmonischen 
Kraftkomponenten der Strahlerfläche. In Wirk- 
lichkeit besteht das Spektrum der Kraftvertei- 
lung aus Komponenten, die mit steigender Ord- 
nung der Teilschwingung abfallen. Es ist aus 
f der inkompressiblen Zirkulationsverteilung zu er- 

Bild 12. Zusammensetzung der am Tragflügel angreifenden mitteln. 
a a ee Da der ‚„inkompressible‘“ Auftrieb für eine 

Flügellänge R 
nn Ca 6b 

A= ou Rl= ER RRgeS RER ONE Sa A (62) 


ist, kann man nach (60) für die Schubkomponente schreiben: 
Ky=-- Bu ro BO SAT OH ee (63). 


Sie setzt sich aus der durch die geometrische Anstellung gegebenen Größe A und einem vom 
Strahlungsfeld abhängigen el ? 


Km — 6 fm Sm er; 


zusammen. 
Den Schubbeiwert definieren wir wieder als 


= 24% ee ee AH 


Sein Verlauf ist nach Bild 11a gleich dem der mitschwingenden Gasmasse!”). 


„Bei kleinen Geschwindigkeiten mündet der Schubbeiwert in den unveränderlichen Wert 
“für die inkompressible Strömung. Wie die Massenkraft am Modell in Bild 6 mit & wächst, so 
steigt auch der Schub im inkompressiblen Bereich mit dem Quadrat der Geschwindigkeit. Sein 
weiteres Anwachsen bis zu einem Maximum und der dann anschließende Abfall kann nach den 
vorausgegangenen Überlegungen als ähnlich dem Resonanzeffekt aufgefaßt werden. Der Verlauf 
der Strahlungs- und Massenwiderstände zeigt, daß bei zunehmender Winkelgeschwindigkeit 
schließlich die gesamte vom Antrieb gelieferte Energie in Form akustischer Strahlung ausgesandt 


wird. 
Zusammenfassung 


Die Widerstände, gegen die der rotierende Tragflügel im kompressiblen Gas arbeitet, 
kommen durch die Überlagerung aller wellenförmigen Strömungen zustande, die von den Ele- 


menten der befahrenen Kreisfläche ausgehen. Die Wirkung der Kompressibilität auf die Gas- 


massenbewegung kann daher im Auftreten einer akustischen Strahlung gesehen werden. 


1?) Das Argument der Strahlungs- und Massenfunktionen Yann auch geschrieben werden: 


wR DIT 
u 
C 
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Aus der Energie des Schallfeldes werden die Kräfte bestimmt, die der Tragflügel erfährt. 
Es sind Massenkräfte, denen die kinetische Energie der Strömung und Widerstandskräfte, denen 
die Strahlungsenergie entspricht. Jede dieser Kräfte besteht aus zwei Anteilen, die durch die 
räumliche Ausdehnung des Tragflügelkörpers und durch seinen Anströmwinkel gegeben sind. 

Alle wesentlichen Eigenschaften der untersuchten Anordnung lassen sich auch an einem 
mechanischen Schwingungsmodell wiederfinden. Dies führt auf eine Näherungsdarstellung für 
die am Tragflügel angreifenden Kräfte und damit zu einer Deutung ihres resonanzartigen Ver- 
laufes. 

Der Auftrieb geht mit wachsender Geschwindigkeit gegen Null, da das Gas der aufgezwun- 
genen Bewegung schließlich elastisch ausweicht. Der Tragflügel könnte dann dem Gas Energie 
in Form akustischer Wellen zuführen. 

Eingegangen 1. November 1949. 


Die Berechnung der Grundschwingungszahlen von 


Spiralfedern 
Von A. Weigand in Eichwalde bei Berlin 


Die Drehmomentfedern elektrischer Meßgeräte, deren Mittellinie mit ausreichender Genauigkeit durch 
eine Archimedische Spirale (s. Bild 2) ersetzt werden kann, will man häufig so bemessen, daß ihre Bigen- 
schwingungszahlen möglichst hoch liegen. Dabei ist vorausgesetzt, daß bei einem bestimmten Drehmoment 
(z. B. 100 mg » cm) ein bestimmter Drehwinkel (z. B. 90°) und eine bestimmte Biegenennspannung (2. B. 
5. kg/mm?) auftritt. Um diese Aufgabe lösen zu können, muß man wissen, wie die Schwingungszahlen der 
Feder in und senkrecht zu ihrer Ebene von den Abmessungen abhängen. Im folgenden wird gezeigt, wie 
wenigstens die Grundschwingungszahlen ausreichend genau berechnet werden können. 

The torgue springs of electric measuring instruments are often designed to keep their natural frequencies 
as high as possible. Within sufficient limits, the centre line of the springs is an Archimedean screw. It is 
required that a fived torque (e. g. 100 mg - cm) causes a fixed torsion angle (e. g. 90°) and a fixed nominal 
bending stress (e. g. 5 kg/mm?). To solve the problem stated above, it is necessary to know the way in which 
the natural frequencies with regard to oscillations in. the direction of the plane of the spring and normal 
to it depend’on the dimensions ofthe spring. T'he paper shows how to compute only the first natural frequencies 
with sufficient accuracy. 

On a souvent l’intention de proportionner les ressoris de torsion d’instruments de mesure electriques de 
maniere, que lours frequences naturelles soient les plus hautes possibles. La ligne mediane des ressorts 
peut Eire prösentde, avec une ewactitude suffisante, par une spirale archimedienne. On y a suppose qu’un 
moment de torsion defini (p. ex. 100 mg: cm) correspondait & un angle defini de torsion (p. ex. 90°) et 
a une tension nominale de flexion d£finie (p. ex. 5 kg|mm?). Pour la solution de ce probl&me il faut connaötre 
la dependance entre les fröquences naturelles, dans le plan et vertical au plan des ressorts, et les proportions 
des ressorts. 

UacTo sKkelaTelIbH0 TaR ONpeneKATb PA3MEPEI KPYTHJBHBIX IPY3KHH 3WIERKTPHYECKUX H3MEPH- 
TeIbHBIX IPHÖOPOB, YTOOB UK COÖCTBEHHBIE YACTOTEI ÖBITH KAK MOIKHO ÖoJlee BeJIuKH. CpenHne 
AMHUM TAKUX IIPy3KUH MOryT ÖBITb C NOCTATOYHOH TOYHOCTBIO IIPeNCTABIeHBI CIIUPA.IBIO 
Apxumena (cm. dur. 1). Ilpı rom Npennonaraetca, YTO HEKOTOPOMy OIPejeIeHHOMy 
MOMeHTy KpyueHusa (Hamp. 100 mr CM) CooTBercTByer OfpenesteHHBIÄ YTOA IOBOPOT&a 
(sarp. 90°) u HekoTopoe ompepereHHoe Harndamınee Hanmpstskenme (Haup. 5&2/mm?). YUroöst 
PeIıuHTb 3Ty 3anayy, HeOOXONUMO 3HATb, KAK 3ABUCHT COÖCTBEHHLIE YACTOTEI TIONEePEeYHBIX 


MH INIOCKUX KONeÖannuf IIPy3KHHBE OT ee P&3MepoR. 


I. Schwingungen in der Ebene der Feder 


1. Zusammenhang zwischen den Verschiebungs- und 
Verzerrungsgrößen 
Wir betrachten einen Punkt -? der undeformierten Federmittellinie mit den Polarkoordi- 
naten r,@ (s. Bild 1). Während der Schwingung verschiebt sich der Punkt P nach P,. Die 
Komponenten der Verschiebung seien u (p, £) in Richtung von r und v (p, t) senkrecht zu r. Ein 
Nachbarpunkt P’ mit PP’=ds geht während der Schwingung in einen Punkt P/ über; 
P,P; =ds, ist das Linienelement der deformierten Stabachse, deren Dehnung durch 


ds 
gegeben ist. Außer der Dehnung tritt eine Krümmungsänderung Ax auf; ist x die Krümmung 
der undeformierten Federmittellinie und x, die der deformierten, so ist 
ER My 2 N RER Ba rd) > 
Die Verzerrungsgrößen r und Ax sind durch die Verschiebungen und ihre Ableitungen auszu- 


drücken, wobei ausdrücklich vorausgesetzt wird, daß die Verschiebungen klein gegenüber dem 
3r 


€ 


« \ 


Radiusvektor r sind. Unter dieser Voraussetzung gilt für die POlATEROEUSERE des Verse | ger | $ 


Punktes P/: Fa 
| 1 erw: RN re ee ARELe | 
he v Be 

: Zr 

Bezeichnet man die Abteilung nach g durch einen Strich, so gilt für die Bogenelemente ds und ds: 

0 (EP) el we ua er ae To EEE 
Der N Far ojtldgs ee: BEST i (6). RL 2 


E% Aus (5) und (6) folgt: | 
Re ru Sera +rV —r'v W-, 
2 (De 
Be Für die Dehnung der Feder-- 
% mittellinie erhält man demnach 


_rWtr)Hr(W —v) D. 


2 2 
OP"r+dr 4 ae R r 
A EN Für den Sonderfall der Archi- 


medischen Spirale mit der Polarglei- 
chung r =agist » 


pure) tW—d 7 Bart 
ap) se 


Setzt man in ())r = R = Const, so 

erhält man den bekannten Ausdruck 

Bild1. Zur Ableitung des Zusammenhangs zwischen Verschiebungs- und für die Dehnung eines Kreisbogens 
Verzerrungsgrößen bei den Schwingungen in der Ebene der Feder bei Verformun 8 der Mittellinie. 


Für die Krümmung der undeformierten bzw. der deformierten Federmittellinie gilt: NEL 


Bi | el Se 
A} ( KR do do? tel eh 8) x 
, “= FPAGIEIE ET Te ee A N Pr 
+) 
a ‚GP ET E22 
bzw. es EN | £ 
| dr,\? der i a 5 
142 (72) Lg: { 2 . "er { Fi % ) x BR 
Be % ET T ee at, 3; a er . (8a). PY ” « 
Kar sn 
\ do, er k er ch, mn > a NE 
; CH Y ? ee, 


Gl. (8a) ist mittels der Gl. (3) und (4) so umzuformen, daß statt Hop Ableitungen Bi ch ER 
2y 3 


die Krane 


dp’ dpi und die Komponenten der Verschiebung auftreten. Es ar 
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- 2 
Für den zweiten Differentialquotiente z a 
Yı 
| der, _ RE, ie r" —)® = 2) r' 
Re AN Fr TE da FR " y) 
dgi +u TE 2 v+2\; „)® ER TE! 
Für den Sonderfall der Archimedischen Spirale lauten die Ableitungen: 
dr, RS 
— =a+U+4- ——- ......2.2.2.0. RE MEHR LIE 
dp, ep Rn 
| d?r, ENT 2 v" 
=u— 0 +— U — ....... NER a WO 
dp; VER RR; e) 
Demnach ergibt sich als Ausdruck für die Krümmungsänderung der Mittellinie der Spiral- 


feder: 
An rt Dur (pPtpu —30+3pV  pHl)v" (m) 
a?(p? nn nr J 
Nachdem durch die Gleichungen Ia und (II) die geometrische Voraufgabe der Berechnung 
der Verzerrungsgrößen aus den Verschiebungen gelöst ist, kann die durch die Deformation in der 
Feder entstehende Formänderungsarbeit (d.h. die potentielle Energie) berechnet werden. 


2. Berechnung der potentiellen und kinetischen Energie 


Während der Schwingung treten in der Feder Normal- und Schubspannungen auf. Da 
die Feder ein gekrümmter Stab von kleinem Querschnitt ist, können die Schubspannungen ver- 
nachlässigt werden. 

Die Normalspannungen entstehen durch die Biegung der Feder und durch die Dehnung 
der Mittellinie. Zwischen dem Biegemoment M, der Krümmungsänderung Ax und der Biege- 
steifigkeit Z.J besteht der Zusammenhang 


M=EJAx. N 

Demnach erhält man für die Formänderungsarbeit infolge der Biegung bei konstantem 

= Querschnitt 
_#EJI 
zu 2 
A; a 55° "7 |(4w EEE NR RN ILL (13). 

Für die Formänderungsarbeit infolge der BE der Federmittellinie gilt 
; A=& [as BE BR ANETTE an: a an (14). 
% Die Integrale sind über die ganze Länge der Felder zu erstrecken; da das Bogenelement 
; der Spiralfeder durch 


ds=ayl-+y?dy 


gegeben ist, ergibt sich für die Formänderungsarbeit A nach (Ia) und (II) die Formel: 


eine —2)w + (g+p)u" —30+ 33V" ev? | 
2a? (9? + 1)9/2 Fo ze 
EF f. [p(u+ v’)+ w — v]? 


| | Bu 77 Be ea ee 


_ Die Integrationsgrenzen y, und 9, sind die Polarwinkel, die den Enden der Feder ent- 
sprechen; F ist der Querschnitt der Feder. 
| 'Bezeichnet man mit u die Dichte des Federwerkstoffes, so gilt für die kinetische Energie 


Be): +) Tatenıap. "N 


1 


(11). 


3: Dehnungslose Schwingungen; Ableitung einer Näherungs- 

KHormel für die Grundschwingungszahl mit Hilfe des Verfahrens 
Fr vonRayleighundRitz 

2 Die Spiralfeder setzt der Biegung einen wesentlich geringeren Widerstand entgegen als 

E einer Dehnung der Mittellinie. Bei den Schwingungen in der Federebene handelt es sich also 

Ri Jin der er um reine Biegeschwingungen; die mit einer Dehnung der Stabachse ver- 
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bundenen Längsschwingungen haben eine viel höhere Frequenz als die Biegeschwingungen und 
können vernachlässigt werden. Wir beschränken uns im folgenden auf dehnungslose Schwin- 
gungen; zwischen den Verschiebungskomponenten vu und v besteht also nach Gl. (I) die Be- 
ziehung: 


d 2) 
’ Pre BAR Er ’ ' Di Fa WETTE 
r(u+V)+Hr(W —ov) =ru+trwW+r vs —H 
bzw. für die Spiralfeder 


o(u+v)+wW —v =pu+u+ 9? @) a ER A 15) 


Für die reinen Biegeschwingungen verschwindet das zweite Integral in Gl. (II]). 

Aus den Energieausdrücken können mittels des Hamilton schen Prinzips unter Berück- 
sichtigung der Nebenbedingung Gl. (15) die Differentialgleichungen abgeleitet werden, denen die 
Verschiebungskomponenten genügen. Dieser Weg soll jedoch nicht beschritten werden, da die 
sich ergebenden Differentialgleichungen kompliziert sind und mit Hilfe von Näherungsmethoden 
behandelt werden müssen, um zahlenmäßige Ergebnisse zu erhalten. Es soll vielmehr mittels 
des Verfahrens von Rayleigh und Ritz eine Näherungsformel für die Grundschwingungs- 
zahl abgeleitet werden. Zur Ermittlung der Eigenschwingungen von der Frequenz » macht man 
den Ansatz 

upt)=f(p)sin2rvt ] 
v(p, t) = fs(p) sin Ar vi j ee BE lee Fon le re (16). R 


Da während der Schwingung keine Energieverlyste auftreten sollen, muß die potentielle 
Energie im Augenblick des größten Ausschlages gleich.der kinetischen. Energie beim Durchgang 
durch die Nullage sein. Aus dieser Bedingung ergibt sich für die Frequenz der Rayleighsche 
Quotient 


Br Tas): Br 
TAnBe ach 
fir+on-e arte Vf —3f+ 39 — (HD 
\ (92+ 1)9'2 
FI: RR 2 BRE RE ARBREN GAR PR RER Er RE ENTE ee} (V). 


In+n (pe+ 1)U2 49 


Der Faktor E/u = .c*ist das Quadrat der Schallgeschwindigkeit c in einem geraden Stab, 
der aus einem Werkstoff vom Elastizitätsmodul Z und der Dichte z besteht. Es ist also 


102 
FAT EEE BRENNT Br 
vu Wäre die Schwingungsform bekannt, d.h. 
hl wären die Eigenfunktionen f, und f, bekannt, so 


könnte man aus Gl.(V) den Eigenwert } berech- 
nen. In Wirklichkeit kennt man aber die Eigen- 
funktionen f, und f, nicht. 

Nach dem Verfahren von Rayleigh und . 
Ritz ersetzt man die unbekannten Eigenfunk- 
tionen durch solche Funktionen, die nur die Rand- 
bedingungen an den Enden der Feder erfüllen; 
man erhält auf diese Weise einen zu großen Nä- 
herungswert für die Grundschwingungszahl. Die 
Feder möge so gelagert sein, daß die Verschiebung « 
in Richtung des Radiusvektors für 9=9, und 
9 =, verschwindet; außerdem soll die Verschie- 
bung senkrecht zum Radiusvektor für g =», ver- 
schwinden. Wir betrachten also folgenden Sonder- 
fall für die Randbedingungen . 


AaW=-09 A@)=I Fe)=0 (18). 


Bild 2. Bezeichnungen an der Spiralfeder Die Formel für den Eigenwert soll noch 
etwas umgeformt werden, Far die Konstruk- 
er Windungen, so ist 


tionsdaten der Feder nach Bild 2 eingeführt werden. Ist n die Zahl d 
%—9ı =2nn ha WEL 
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"Ferner folgt aus , = a9, fa = a9: 


Für 9, und 9, gilt 


9ı PEBEE 9%,=2nn ee EURE (21) 
Durch 
; P=9+a—p) 2 =2nn (a t40) .. 2.2.2... (22) 
mit 
Ti; $ 
Oi ee ua Amor (23) 


führen wir in Gl. (V)eine neue Integrationsveränderliche zein. Bezeichnet man die Differentia- 
tion nach % wieder durch einen Strich und setzt noch zur Abkürzung 
1 


a a el LEN 24), 
7 INT Er 


so geht Gl. (V) nach einigen Umformungen über in 
NZ 22 


m Er ST ee (25), 
= 2 2% 97? i 2 lo 9 
fe 
#7 © 
a, er 
N fe+r )(@+0) (4 ) dx SUR RAS (26). 


In dem Integranden des Zählerintegrals Z ist zur Abkürzung 2 = x + o gesetzt; die Inte- 
grationsveränderliche ist jedoch z, das zwischen den Grenzen O und 1liegt. Für die Eigenfrequenz 
erhält man jetzt aus Gl. (17) wegen Gl. (20) 

1 c2J 


2= _—_ 22 DIREKT Va). 
a RT BEER (Va) 

Für die Eigenfunktion f,(x) machen wir den Ansatz 
ee Ras) re. (27), 


der die Randbedingungen G]. (18) erfüllt und noch von dem Parameter & abhängt. Die Eigen- 
funktion f,(x) kann nicht mehr frei gewählt werden, da sie durch die Bedingung der Deh- 
nungslosigkeit der Schwingungen mit f,(x) verbunden ist. 

Nach Gl. (15), (16) und (22) ist: 


nt fit zittert) - 


f c+ 0 
Die Lösung dieser a Elan die die Randbedingung f,(1) = 0 erfüllt, lautet: 
d 
ei mt 
3 NEWIR).,- 5 
J.(®) = (+ 0) —— di... 2.20. (28). 


r | (2+.0)? 
- % 
‚Setzt man in diese Gleichung den Ansatz für f,(x) aus Gl. (27) ein, so erhält man 


av | 122% 1+ 
Aa) = na (a4 HH) —elitoin tt 
er elta 
29). 
te! = une (23) 
In Ui-+e zu =+o 
Feet = 0@+30)1— aus \ 
E- +el- een , HH BEIN Yahgee 2) | 


Mit (27) und (29) wird der Rayleighsche Quotient (26) eine Funktion der Größen n, 0 
und a; die erste Näherung erhält man, wenn man & = 0 setzt. Eine bessere Näherung ergibt 
sich, wenn man & so bestimmt, daß der Rayleighsche Quotient möglichst klein wird. 

Wegen der komplizierten Gestalt der Formeln (26) und (29) wird die numerische Ermitt- 
lung des ersten Eigenwertes ziemlich umständlich, Eine wesentliche Verein-. 
fachungergibtsich, wenn man die Tatsache berücksichtigt, daß 
die Windungszahl der Spiralfedern fast stets mindestens gleich 
2, häufig sogar mindestens gleich3 ist. Die in dem Zähler und Nenner des 
Rayleighschen Quotienten nach Gl. (26) vorkommende Größe n? nimmt nach Gl. (24) 
folgende Werte an: 

n = 2 Windungen: Mr — 0,00633, 


n = 3 Windungen: n? = 0,00281 . 


Man kann also die Gl. (26) durch folgende NER PUng se TUNER ersetzen: 
1 


Pa ERON | AS 
ge@at N ES TOBEN 
1 . \ 


N*= [SH] (a+o)de 
Außerdem kann man in dem Andatz für die Eigenfunktion f,(z) [Gl. 29)1 das Glied mit 


dem Faktor 5 = gegen das mit dem Faktor 2 na vernachlässigen. Man ersetzt also J.(x) nähe- 
rungsweise durch 1 ONLERE 
(a) = 1 1 1 ii 
@) +o)1—n)—e(l+o)i- 5 
(29a). 
1 
+ + — oa al 


Schließlich kann man noch im Nenner des Raylei ghschen Quotienten f,(x) gegen 
3’(&) vernachlässigen, da fy” nach (29a) den Faktor (An)? enthält. Für den ersten Eigen- 
wert A, der jetzt A*? heißen möge, ergibt sich schließlich folgende einfache, für n> 2 gültige 
Näherungsformel 


le — art ala? — ap 


d (2+ 0)? r GE 5 
RZ — 2. rn (80) 5 
Cum)? IHRE) + a HP (a+ ed ‘ | 
mit Me 
. 1 2 Sr 1 \ x u 
H,(e) = — ——— Ha tod ee = | 2 
SAlsby: | 
1—a2? 1 En RN 
H;(a) = — 7 +2 Be - a N 
Als Funktion von « hat der in Gl. (30) Te Quotient die Form # Ba £ 3. & 
+2B,0-+0,0%2 ET REHEN, - : 
Bra 2a 4ı+2B,%+ ya? ir ie DE A 
RR? a Sl a 
mit | 3 ee 
1 x 2 
. [(2 — 22)? NR er—aye 
= dx, B, = I —  — dr, = 
en RER y ; (@+0)° ie NE BET 2 


berechnet. 


Die erste Näherung ergibt sich für « =0 und ER eich et 
Näherung führt auf die quadratische Gleichung a s 


(B,0, — B,0,)0°+ (A, Ü ER 


er 


A a Re 
2 DERNRER: aa 
RN A ER 

ru + + 


Fe N zZ j - . N 
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von der diejenige Wurzel zu nehmen ist, die den kleinsten Wert des Quotienten liefert. Zahlen- 
beispiele sind im Abschnitt III behandelt. 


II. Schwingungen senkrecht zur Ebene der Feder 
Aus den Gleichgewichtsbedingungen am Federelement und dem Zusammenhang zwischen 
Schnitt- und Verzerrungsgrößen soll die Differentialgleichung abgeleitet werden, der der Schwin- 
gungsausschlag genügt ;für Windungszahlen, die mindestens gleich 2 sind, gelingt die näherungs- 
weise Integration durch Besselsche Funktionen von der Ordnung + 2/3. 
Wir betrachten ein Element ds der Mittellinie der Feder; o sei der Krümmungsradius an 
der betreffenden Stelle (s. Bild 3). Mit 


a En a (34) 
lauten die Gleichgewichtsbedingungen des Federelementes 
EOWL- dM & dt M 
Is g=0, -5+9- 0, vn FREIEN (35). 


In diesän Gleichungen ist Q die Querkraft, M das Biege- und T = Torsionsmoment; g ist 
die (senkrecht zur Federebene wirkende) Belastung der Feder je Längeneinheit. Ist w der 
Bene nsschleg und «’ die Masse der Längeneinheit, so ist 


KR) 
ne } Aare 0): 7 gg0% 
ot M*+dM 
Mit 
ps=ayl-+g9?do 
ü= Au BEER, —=q NR N) T+dT 
"nr? +2r?%—rr' 2-+9° n , 
/, 
nehmen die Gleichgewichtsbedingungen folgende Form an N 
i 4 
d 
+ ayl-+g?g=0 
aM g°+2 
LE RL® 1 =D 
ER T+ayl+pQ = 2 (292). 
. 
N - dT p’+2 4 —( Bild 3. Schnittgrößen am Federelement für 
j do p°+ 1 Schwingungen senkrecht zur Ebene der Feder 
; Hat die Feder mindestens 3 Windungen, so ist der Polarwinkel @ wegen Gl. (21) stets 
wesentlich größer als 1. Bezeichnet man die Ableitung nach @ durch einen Strich, so kann man 
| demnach die Gleichgewichtsbedingungen ausreichend genau schreiben: 
Q’ +apg=0 
MEI ea Re (35b) 
T’+M = 
4 Wir betrachten nun den Zusammenhang zwischen der Verformung der Feder und den 


Schnittgrößen M und T. Es sei y der kleine Winkel, den die Tangente an ein aus der Feder- 
ebene ausgelenktes Element mit der Horizontalebene bildet. ® sei der Verdrehwinkel. 
Mit der Vertikalverschiebung w eines Federelementes hängt y gemäß 


TA 
- BISedgE Yretr2 
zusammen; für die Spiralfeder ist . 
w' w' 


afitp ap 
Für die Änderungen der Winkel y und @ gilt (s. Bi&zeno-Grammel, Technische 
+ ei, S.339), wenn in den dort für den Kreis abgeleiteten Formeln r durch o und dydurch dx 


ersetzt wird: 
ed: € Bien, 
dy = ueerr - (er Br te ) dp, 


h r-+r 12 
Tody 


72-272? —rr” 


IT aan ) 
rn) a. 


TE en CH 5, 
Ex I N 
r En ai 
$ 


» 
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Für. die Spiralfeder erhält man: 


aM 92-2 aM 
Da ' ii nee (39) 
REN DI 2,00 Nodh 
und daraus 
ee, 
er SR EN (40). 
T = — (W-+y) 
ap 


Eliminiert man die Querkraft aus den Gleichgewichtsbedingungen (35b),: so erhält man 
—oM"’+-M+oT— T+ayg=0. 


Aus dieser Beziehung, aus der Gleichgewichtsbedingung 7’ +M =0 und aus Gl. (40) 
ergeben sich zwei simultane Differentialgleichungen für die Verformungsgrößen y und ®: 


3 3 2 
BI— 9) +" 9); P—W)l+@J410"+Yy —- +yp)]+arptg=0 

p p F (41). 
EJ(y —d)+GJd [9 + y’ 7," —0 


Ersetzt man y gemäß Gl. (38) durch die Auslenkung w(9, £) des Federelementes, ferner q 


durch die Trägheitskraft — ib und schließlich noch w durch y(p) sin ot bzw. 9 durch 9(p) - sin wi, 
so erhält man schließlich folgendes homogene lineare Differentialgleichungssystem für y(@) 


und (p): 


EJ 6EJ Be ER, er wi: 
Be a Vi RE ja „ ae "+ato w* 
m ae PR gp° oO gi op? Yy p?um?y 
+(EI+0J0) a ad er ad =0 . (42). 
BEER RL, y.eja gr sap ARE 
p p p p p 


Konz man die in 5 (22), (23) und (24) eingeführten Bezeichnungen, so geht (42) in 
15EJ GJd vr fI52J 3@Jd : 
(BR ae fit ey 


ze 
zZ 2 N 22 


5 + 
a 14 Has (ra —n) 9’ — 728 EJ +20 Jd) (ra — )® 


i i +n:3EJ(, —n)29 =0 (42a) 
@Jda+EJ „ EJr2GJd ,„ MJdtr, —r)=, G6Jd(r —r) 
Be Wand er, en) 


nz 0 
über; die Ableitung nach z=x-+-.o ist durch einen Strich bezeichnet. 

Es ist unwahrscheinlich, daß das komplizierte Differentialgleichungssystem (42a) durch 
bekannte Funktionen integriert werden kann. Es wurde deshalb folgendes grobe Näherungs- 
verfahren benutzt. 

Man vernachlässigt in (42) alle Glieder, die den Faktor 7? oder 7° haben. Man erhält aus 
der zweiten Gleichung (42a) 


»=0 
und aus der ersten: 
er I H pael &.0.. ee a (43) 
mit } 
9. P) w 2 zz r;) 2 
A= (Zn) I RER NE RE (44) 


 Z. angew. Math, Mech 


Ba. 31 Nr. 1/2 Jan,/Febr. 1951 Weigand, Die Berechnung der Grundschwingungszahlen von Spiralfedern 48 


Durch die Vernachlässigung der Glieder mit. den Faktoren n? bzw. n 5 wird das System (42a) 


- auf die Differentialgleichung zweiter Ordnung (43) reduziert. Infolgedessen entfallen die für 


(42a) zu stellenden dynamischen Randbedingungen, die sich auf das Verhalten der Schnitt- 
größen an den Federenden beziehen; es bleiben nur die geometrischen übrig, die aussagen, daß 
der Ausschlag an den Enden der Feder verschwindet. 

Gl. (44) soll durch Einführung der Dichte « statt der Masse der Längeneinheit «’ noch 
etwas umgeformt werden. 

Ist # der Querschnitt der Feder, so ist « = uF. Fürrechteckigen Querschnitt ist ?=b-h 
und Jd =kbh°®, wobei k ein von dem Verhältnis a/h abhängender Zahlenfaktor ist, der aus 
Tabellen entnommen werden kann; für b/h = 10 ist z.B. k = 0,312. Mit 


ES 
u 
nimmt dann Gl. (44) die Form an 
f —r.\4 
2 (ana ar. U SER TR ALTEN: 


Es wurde darauf verzichtet, die Näherungsrechnung durch eine Störungsrechnung zu ver- 
bessern, da ihr Ergebnis hinreichend genau mit experimentell gewonnenen Werten überein- 
stimmt. 

Aus (43) erhält man leicht den Rayleighschen. Quotienten für den Eigenwert 42. 
Man kann'nämlich Gl. (45) mit z2= x +0 auch schreiben 


(on) +re@toy=0 De A 93 (43a), 


Man multipliziert diese Differentialgleichung mit %, integriert zwischen 0 und 1 und formt 
das erste Glied durch Teilintegration um, wobei zu beachten ist, daß der Ausschlag an den Enden 
der Feder (e=0 und x =1) verschwinden muß. 

Man erhält so: 


Y 
[ (ro 
0 


; . (45). 
/ (et o)y’de 
Setzt man in (45) für y einen die Randbedingungen erfüllenden Ansatz ein, z.B. 
y=a(ll—xa)(l+ax), 
so erhält man mit & = 0 die erste Näherung für den ersten Eigenwert, die man verbessern kann, 
indem man & so bestimmt, daß der Rayleighsche Quotient möglichst klein wird. 
Auf Grund der Zusammenstellung in Jahnke-Emde, Tafeln höherer Funktionen, 
4. Aufl., S. 150 kann die Differentialgleichung (43) durch Besselsche Funktionen integriert 
werden. Man hat in der ersten der dort angegebenen Differentialgleichung zu setzen: 
A 2 
a“=2, Bp=2, y=3, p=t7. 

Die allgemeine Lösung lautet also, wenn mit A und B die Integrationskonstanten be- 

zeichnet werden: 


y-al4rnl #)+Byln(b: ). DENE N Ur 347% (A6), 


Aus den Randbedingungen y(o) = 0 und y(1 +.) = 0 ergibt sich die Gleichung für den 
Eigenwert: 


Jays (222). re (+ e)| = Ianlee ul +0») Ba 8 (AT) 


Da keine Tafel für J 497 (x) zur Verfügung stand, wurden bei den Zahlenbeispielen zum Teil 
die Potenzreihen, z. T. die asymptotischen Reihen von Hankel verwendet. 

Die Rechnung wurde auf die Ermittlung des ersten Eigenwertes beschränkt. Aus Gil. (47) 
erhält man für ihn ein genaueres Ergebnis als aus dem Rayleighschen Quotienten, da in 
Gl. (47) nur der Fehler steckt, dessen Ursache die Annahme großer Windungszahlen ist. 


III. Beispiele 
An einigen Beispielen soll gezeigt werden, daß die Grundschwingungszahlen von Spiral- 


federn auf Grund der hier durchgeführten Berechnungen ausreichend genau ermittelt werden 
können. Es handelt sich um Drehmomentfedern aus, Bronze elektrischer Meßgeräte, deren Eigen- 


er." 


FOREN 4 
S 


ur 


. x see a Tg Nano, Main Me 
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schwingungszahlen durch Aufnahme der Resonanzkurven auf dem Schütteltisch experimentell 


ermittelt wurden. 
Die Materialkonstanten ce und c’ haben die Werte: 


c = 3,59 - 10° mm/sec 
c’ = 2,27 - 10° mmj/sec. 


1.Schwingungeninder Federebene. Br, 
Wir betrachten zunächst eine Feder mit den Abmessungen 
b=0,5 mm, h = 0,05 mm, 1,8 mm, 3020; n = 5 Windungen 
Daraus berechnet man 
Y; 
Ya—t; 
Wegen Gl. (Va) und Gl. (30) lautet die Formel für die en 


=.0,778: 


DI 


ee — 22+ a(2? — 2°)? 
ya es N 0 eb ; 
auf Tg 1 S 
TEN Fna)+aHsle)a+ 0)rde 
Aus den gegebenen Zahlenwerten findet man 
1 ch 
872 Y3 n(ra — ri)? h 
Für die bei der Berechnung des Rayleighschen Quotienten vorkommenden Integrale. 
erhält man aus Gl. (32): 
A, =, 0184, B.:=.0,00757,,.,0, =.0,003773 
As = 0,00647, =..B,:= 0,00353,2.0, —=0,00195. 
Die quadratische Gleichung für & lautet nach Gl. (33) 
1,45 &2+11,5%+16,1=0. 


Die Wurzel, die das Minimum des Rayleighschen Quotienten liefert, ist« = —6,12. 
Für den Minimalwert findet man 1,845, demnach ist die Grundschwingungszahl 


v = 12,97 y1,845 = 17,6 Hz. 


In erster Näherung (& = 0) ergibt sich v = 12,97 \ we = 21,9 Hz. Gemessen wurde 


eine Grundfrequenz von rd. 17 Hz; sie stimmt praktisch mit der zweiten Näherung überein. 
Als zweites Beispiel betrachten wir folgende Feder: 


b=06mm, h=0,06mm,  », =5,25mm,: 1, —15mm, 'n = 7,25 Windungen. 


=12,97 ec; 


Daraus folgt: 
0= 0,4, 


Für die Grundschwingungszahl erhält man aus der im ersten Beispiel angegebenen Formel 


1547. 


falls mit Z bzw. N Zähler bzw. Nenner des Rayleigh schen Quotienten bezeichnet werden. 
Durch Berechnung der Integrale findet man 
Z 0,0621 + 0,0458 x + 0,0106 «? 
N 0,0033 +0,00384& +0,01 02° E 
Für « = —5,3 nimmt dieser Quotient den Minimalwert 7,5 an; die zweite Näherung e era. 
gibt also 2 
v — 15,4.775 = 42,5 Hz. 
In erster Näherung (« =0) erhält man 
0,0621 


y—= 15,4 0,0035 7 Hm 
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| - Gemessen wurde eine Grundfrequenz von rd. 50 Hz; die zweite Näherung liegt alsotiefer 
. als der gemessene Wert. Der Grund für diese Erscheinung liegt in folgendem. Bei & = 0 ist 
- die Drehmomentfeder mit der Achse des beweglichen Organs des Meßgerätes verbunden; dort 
muß also eine dynamische Randbedingung erfüllt sein, die aussagt, daß das von der Feder auf 


das bewegliche Organ wirkende Moment mit dem Massenmoment im Gleichgewicht stehen muß, 
das bei der Schwingung des beweglichen Organs entsteht. Diese Bedingung, bei deren Formu- 
lierung der Eigenwert selbst vorkommt, wird von dem hier benutzten Ansatz nicht erfüllt; es 
war für die Eigenfunktion f,(&) bei x = 0 keine Bedingung vorgeschrieben. Einen oberen Grenz- 
wert für die Eigenfrequenz wird man erhalten, wenn man den Ansatz so wählt, daß f,(0) =0 
wird. Man erhält in diesem Fall in zweiter Näherung den Wert » = 70,6 Hz für die Grund- 
schwingungszahl, der wesentlich höher liegt als der gemessene. 
2. Schwingungensenkrecht zur Federebene 
Zunächst wurde für mehrere Federn die Grundschwingungszahl in erster und zweiter 
Näherung nach Rayleigh ermittelt. Für alle Federn war 
1, =8mm, rn =3omm, 7 0= 0,778. 
Aus dieesn Werten und aus dem Näherungsansatz y = z(1— x)(l1 +xx) für die erste 
Eigenfunktion erhält man für den ersten Eigenwert nach Gl. (45) 
92 _ 0,277 + 0,1542 & + 0,0483 a? 
— 0,0426 +0,05 & +0,013402° | 
In erster Näherung ist also 
\ 0,277 
Pe ER ek >= 
4 0.0126 6,5, 1=2,5. 
Die Berechnung der zweiten Näherung führt auf die quadratische Gleichung 
0,57 &® — 16,5 — 29,5 =. 
Die in Frage kommende Wurzel ist & = 30,7; für diesen Wert erhält man 
IESBR  2 FER 
Aus (44a) erhält man für die Grundschwingungszahl 
1  yYkch h 
= — ——41=2840 Yk -A. 
? Am: (Va —r)?n 4 Yk 7 2 
In der nachfolgenden Zahlentafel sind die berechneten Werte mit den gemessenen ver- 
glichen. 
n ® 2) e 
b h fin. t eite 
ß Win- erste zwei 
, dungs- Näherung | Näherung | Messen 
% mm mm au Hz Hz Hz 
B.\ 2 0,05 0,333 = 41,6 31 30 
| 0,99 0,099 0,312 8 49,7 37 34 
re 0,84 0,084 0,312 7 48,7 36,2 32 
L 0,93 0,093 0,312 \8 46,8 34,8 3l 
; 0,76 0,076 0,312 - 7 44,1 32,8 30,5 
® 0,68 0,068 0,312 6 45,8 34,2 30,5 
” 0,37 0,037 0,312 4 37,4 27,9 25 
RX 0,61 0,061 0,312 6 41,2 30,7 28 


\ Für die Feder mit den Abmessungen 
U b=06mm k=006mm »,—525mm n=1ömm n= 17,25 Windungen 
_ wurden außer den Schwingungen in der Federebene auch die Schwingungen senkrecht zur Feder- 


> ebene untersucht. Man erhält 


ae 


1.-Vkch 4 
ee m in 3 a ’ == 18,9 ng . 
ER ERTSEEN 18,9 sec v A 


« EEE ER a N ie a ee at un 
EB N ae Sa 
/ u Br Zr ü es Pe 
\ Sr . z fr - 
Pe, 3 Du, I 
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findet man aus Gl. (45) 1,452, + 0,624 «+ 0,163 a? 


DEN eh. 
A 0,030 + 0,0324 «+ 0,00977 a? 
In erster Näherung erhält man Be ZA, 
1,452 ST = 
2 — Des = — | ‚ 
4 0,030 48,4, ) = 6,96 


Die Berechnung der zweiten Näherung führt auf x = 12,5 und ergibt 
12 21954 4 =4,05. 

Demnach ist die erste Näherung für die Grundschwingungszahl 
v = 18,9 - 6,96 = 131 Hz 


v= 18,9 -4,05 = 76,6 Hz. 


"Gemessen wurde eine Frequenz von rd. 65 Hz, die noch um 18% von der berechneten 
abweicht. Schließlich wurde für diese Feder der erste Eigenwert A mittels der Besselschen 
Funktionen aus Gl. (47) berechnet. Es ergab sich 


1=3,10. 
Mit diesem Wert erhält man für die Grundschwingungszahl 
ER ER Zr 
Dieser Wert weicht nur noch um rd. 8% von dem gemessenen ab. Schließlich betrachten 
wir noch eine Feder mit den Abmessungen 
b = 0,6 mm, :h =0,06mm, k = 0,312 
„=8mm, n=6mm 0 > gi 


Mit PER 
Bis _Vkeh — 160sec-1 
4: n(r, —r;)? 


und die zweite 


erhält man 
v= 1604. 
Aus Gl. (45) ergibt sich 


22 0,00839 + 0,00694 & + 0,00263 a? 


0,1167 + 0,1190 & + 0,0346 0 
% = 1,5 liefert das Minimum des Ra yleighschen Quotienten, für den man den Wert 


42 = 0,0662 
erhält; die erste Näherung ergibt sich aus 
0,00839 
Pe Zn 
0,1167 Nach 


Für die erste Näherung erhält man 


v = 160 Y0,072 = 42,8 Hz 
und für die zweite / 
v = 160 Y0,0662 = 41,2Hz. 
Berechnet man den ersten Eigenwert nach Gl. (47), so erhält man A = 0,256 und 
v = 160 - 0,256 = 41 Hz. Alle berechneten Werte unterscheiden sich nicht wesentlich vonein- 
ander und von dem gemessenen Wert, der 35 Hz betrug. 


IV. Zusammenfassung 

Für Federn, deren Mittellinie die Form einer Archimedischen Spirale hat, wird die Grund- 
schwingungszahl für Schwingungen in und senkrecht zur Federebene mittels des Verfahrens 
von Rayleigh näherungsweise berechnet. Es wird dabei vorausgesetzt, daß die Feder 
mindestens 3 Windungen hat. Für die praktisch besonders wichtigen Schwingungen senkrecht 
zur Federebene wird außerdem die Differentialgleichung für die Eigenschwingungen aufgestellt, 
die sich für große Windungszahlen durch Besselsche Funktionen von der Ordnung +2/3 inte- 
grieren läßt. 

An einer Reihe von Beispielen wird durch Vergleich mit gemessenen Werten gekeigt, daB Gr 
sich die Gr mittels der entwickelten Gleichungen ausreichend genau 2 
berechnen lassen. j 
Eingegangen am 24. 10. 1949, 


x 
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Eine neue Interpolationsmethode zur Ermittlung der gesamten 
Zerstreuungsfigur zentrierter optischer Systeme mit unendlich 
fernem Bild 


Von H. Köhler und @. Pradel in Jena 


Es wird ein Rechenverfahren mitgeteilt, das zur Ermittlung der gesamten Zerstreuungsfigur der aus einem 
optischen Instrument für subjektive Beobachtung austretenden- Strahlenbüschel dient. Die Anwendung des 
Verfahrens hat nur die irigonometrische Durchrechnung der im Meridionalschnitt verlaufenden „Koma- 
strahlen“ zur Vorausseizung. Aus ihren Abweichungen gegen den Hauptstrahl wird die Abweichung eines 
beliebigen Strahls innerhalb der Pupille, also auch der Sagittalstrahlen, durch Interpolation gefunden. Das 
Verfahren beruht auf einer Anwendung der Theorie der Bildfehler fünfter Ordnung. 


A caleulating melhod is described, which is used for finding the entire diverging figure of ray pencils 
emergent from an optical instrument for subjective observation. The use of ihe method is based merely on 
the Irigonomelric caleulation of Ihe ”coma rays‘‘' passing along the meridional section. T’heir deviations from 
the principal ray permit to find by interpolation the deviation of any ray within Ihe pupil, and therefore also 
of the sagittal rays. T’he method is based on an application of the theory of image errors of the fifth order. 


On a deerit un procede de calculation qui s’applique a la determination de Ventiere figure de divergence 
d’un faisceau de rayons Emergents d’un instrument d’optique pour observalion subjective. L’application du 
procede est basee simplement sur la calculation trigonometrique des ”rayons coma‘‘ passent par la section 
me£ridionale. Leurs deviations relatives au rayons principal permettent de trouwver par interpolation celle 
d’un rayon' quelcongue dans la pupille et, par consequent Egalement celle des rayons sagittau®. Le procede est 

basee sur l’application de la theorie des erreurs d’image du cinquieme ordre. 

OnnesipaeTca TEeopeTuyecknä MeToA ONpeHeneHuA burypsr paccesung IyyRa Iyyei, 
BEIXONAIHUX 3 ONTNYEeCKOTO HPHOOPA MIA HEeIMOCPeACTÜCHHEIX (CYOBEKTUBHEIX) HAÖMIONeHUH. 
Ilpumeneune 3TorO MeTona OCHOBEIBAETCA TOABKO HA TPHTOHOMETPHYeCKOM BEIYHCIHEHHH 
ayyeii „‚KoMa‘‘, UMYIUMX BAONb MEPHNMOHANBHOTO ceyeHus. Ms uXx OTKIOHCHHUA OT TIABHOTO 
ayya OlpenenaertcH, Ipu IOMOINH UHTEPHONANHN, OTKIOHeHHE AIWÖOTO AyYa B Hpenenax 
3payka, T. e. TaK3Ke H CATUHTAJIBHOTO JIyYa. ITOT METOA OCHOBBIBAETCAH Ha TPHMEHECHHH TLOpHH 
HeTOynocTei M30ÖpaskeHnä IATOTO MOPAAKA. 


I. Einleitung 
Nach neueren Untersuchungen, die im Zeißwerk [1] durchgeführt wurden, hat die Kenntnis 
der vollständigen Zerstreuungsfigur bei optischen Instrumenten zur subjektiven Beobachtung 
entscheidende Bedeutung für die Beurteilung der Bildgüte. Das gilt besonders für Feldstecher 
mit Gesichtsfeldern über 70°. Wegen der großen Einfallswinkel an den einzelnen brechenden 
Flächen lassen sich bei diesen modernen Systemen im allgemeinen keine exakten Schlüsse aus 
der bisher üblichen Durchbrechung im Hauptschnitt auf die gesamte Zerstreuungsfigur, also 
besonders auf ihre sagittalen Abweichungen ziehen. Auf die Notwendigkeit, die gesamte Zer- 
streuungsfigur für die Beurteilung heranzuziehen, hat im übrigen schon M. Herzberger [2] 
hingewiesen. 
Zur Ermittlung der vollständigen Zerstreuungsfigur ist normalerweise die Durchrechnung 
mehrerer windschiefer Strahlen erforderlich. Diese stellt einen immerhin beachtlichen Arbeits- 
| aufwand dar. Es wurde deshalb von H. Köhler [3] bereits ein Interpolationsverfahren aus den 
Seidelschen Bildfehlern 3. Ordnung abgeleitet, das jedoch nicht in allen Fällen anwendbar ist 
und im wesentlichen zur schnellen Orientierung während der Korrektionsarbeit gedacht ist. Das 
Bedürfnis, aus den Durchrechnungswerten des Hauptschnittes und mit einem Minimum von 
Rechenaufwand die gesamte Zerstreuungsfigur mit einer Annäherung von wenigstens 10 bis 20% 
zu ermitteln, führte zu dem im folgenden beschriebenen Interpolationsverfahren. Dieses Ver- 
fahren ist zunächst für Flächenfolgen mit unendlich fernem Bildpunkt, also im wesentlichen für 
Geräte zur subjektiven Betrachtung, entwickelt worden. Die Abänderung der Formeln für 
Geräte mit endlich fernem Bildpunkt bleibt einer späteren Arbeit vorbehalten. 


II. Grundgedanke des Interpolationsverfahrens 
Das nunmehr im einzelnen beschriebene Verfahren besteht aus vier Abschnitten und zwar: 
i 1. Aufstellung einer Beziehung zwischen den meridionalen und sagittalen Abweichungen 
des austretenden Strahles gegenüber dem Hauptstrahl in Abhängigkeit von den Koordinaten der 
_ Austrittspupille, wobei die Beziehung eine Potenzreihe mit Gliedern bis einschließlich der 5. Ord- 
nung darstelle. 
2. Bestimmung der zunächst unbekannten Koeffizienten dieser Beziehung aus den Rechen- 
ergebnissen im Hauptschnitt. 
3. Ermittlung der zunächst unbekannten Koordinaten der Austrittspupille in Abhängigkeit 
# von denen der Eintrittspupille, die im allgemeinen gegeben sind. 


Te a 
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4. Die eigentliche Berechnung der Zerstreuungsfigur mit den nunmehr ermittelten Be- 
stimmungsstücken der unter 1. aufgestellten Beziehungen. 


III. Aufstellung der Fehlergleichungen 
Bei der Aufstellung der Fehlergleichungen kommt es darauf an, von vornherein diejenigen 
Potenzprodukte der Pupillenkoordinaten auszuscheiden, die auf das Zustandekommen der Bild- 
fehler ohne Einfluß sind. Man gewinnt eine solche Beziehung am elegantesten aus dem Eikonal. 
Bezeichnet man dieses mit V, ferner die meridionale Koordinate der Austrittspupille mit m’, 
die sagittale Koordinate mit M’, entsprechend die Koordinaten des Objektpunktes mit y und z, 
dann ist das Eikonal wegen der Drehsymmetrie des Systems darstellbar in der Form 


V=V (2%, ym -F2zM, m2 + U) NEN EN SEES: 


Legt man den Objektpunkt in die Ebene z=0, dann vereinfachen sich die unabhängigen Variablen 


zu 9°, ym’, m’®-- M’?. Von ihnen sind innerhalb der 6. Ordnung folgende Kombinationen 
vorhanden: 

y  ym’ m’®+ M'’? 

vom em) gem ya mer) (me Ma 

yon yma Me) Yin ma) ya aa | 3 © = ©) 

ym®  yPm?(m2+M2)  ym’(m®4M2)  (m’®4+ M%): 
Nimmt man das von nun an feste y in die unbekannten Zahlenfaktoren der Potenzprodukte 
auf, dann läßt sich Gleichung (1) in Annäherung 5. Ordnung so ansetzen (n'—Brechungsindex): 
1 = Vor Pam’ + cos w|S mat m + Bm + Im’ mt 04T mi 

. (3). 

ma 4 MB + Sm ma 4 22) + Sm’ (me + MEHR (me+ My 
Cos w’ und die Zahlenfaktoren sind wegen späterer Bequemlichkeit eingeführt. w’ist die Austritts- 


neigung des Hauptstrahls. Um vom Eikonal auf die Abweichung des austretenden Strahles zu 
kommen, benutzt man die Brunsschen Beziehungen 


Richtg.der m-Achse aV 


EI MED. 
Richtg. der M-Achse Eyes a a (4), 
Strahlrichtung 
oV 
=n Nr ti 
Bun: 6). 


wobei P’ und @’ die auf die Systemachse be- 

Richtg. der Systemachse zogenen Richtungscosinus des austretenden 

Strahls’ darstellen. Diese Richtungscosinus 

drücken wir nun aus durch «‘; bzw. u}, die 

Winkel, die der in die Ebenen M'=0 bzw. m’—=0 projizierte Strahl mit der Systemachse 
bildet (vgl. Bild 1, hierin ist sing = P’ und siny=@)). 

Diese rein geometrische Betrachtung ergibt!) 


Bill ı 


PN —-smwsimw=sinwWcow..:. 2... 00. (6), 


O1 —sinweinwW=sinucosW :.. 0... (7). 
Da wir eine unendlich ferne Gaußsche Bildebene voraussetzen, sind die Größen u, — w' und u), 
klein von mindestens 3. Ordnung in bezug auf die Systemachse. Infolgedessen Ben mit Ge- 
nauigkeit 5. Ordnung die Taylorschen Entwicklungen 
sin u, = sin w' + (u) — w') cos w — Hs 
cos uU, = cos w' — (u, — w') sin w' — *** 


sinyw=Ww— re Sala (9.8 


csmw=1— '': 
1 —sinwsi®,=1— 


1) Zuder an dieser Stelle gewählten Darstellung würden wir durch eine Bemerkung von Dr. Boe gohol d.4% 


angeregt, der den Entwurf zu dieser Arbeit freundlicherweise einer Korrektion USERADBER, hatte, 
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Setzt man (8) in (6) und (7) ein und läßt wieder Größen höherer als 5. Ordnung weg, dann erhält 
man 


E Fan ls —Ww)cosWN. Ai (9), 
EN De (10). 
| Führt man nun die Differentiation nach (4) und (5) aus, dann entsteht 
= Vıtcosw Am’ +— (m®4 9) +Om®+-- PN A BEN 
10V 
a 5 EN SE ER  l  ee aan —4 4 ER Ra EEE (12). 
Für m’—= M’=0 folgt aus (11) und (9) 
N TE NR ee Re (13). 


Setzt man schließlich (9) und (13)in (11), (10) in (12) ein, kürzt mit cos w’, schreibt m ?+ M?—p'®, 
geht vom Bogen- zum Gradmaß über, multipliziert mit 100 und nimmt rechts die davon 
herrührenden Zahlenfaktoren in die noch unbekannten Koeffizienten A, ...., H auf, dann ergeben 
sich die Endformeln 


| 100 (u — 0) Am’ +9 p°+ (C+ DJm? + (B+ mp? 
. (14) 


4 
+ Fm®+GE + @m’?p'?-+ Hm’p’* 


N und 


100 u, = (B+ Om’ a NER LM REN. (15) 


die sich nach Einführung der Abkürzung 
K=Cm'+Ep?®+Fm®?+Gmp? Hp“ ..... 2.2.0. (16) 
schließlich so schreiben lassen 


wm -u+ mw roR + Dme tr 2 RT eva 
| und 
F- | 1004 —=(B+K)M' | EU BF RAN SER, (18). 


‘ Damit ist der erste Teil der Aufgabe gelöst. 


Die Beziehungen zwischen den meridionalen und sagittalen Winkelabweichungen des aus- 
tretenden Strahl sind als Funktion der Pupillenkoordinaten gefunden. 


IV. Berechnung der Konstanten der Fehlergleicehungen 


‚Die Koeffizienten der Gleichungen (14) und (15) sollen aufgabengemäß aus den Ergebnissen 
der Durchrechnung im Hauptschnitt gewonnen werden. Als Bestimmungsstücke stehen einmal 
die Strahlabweichungen gegenüber dem Hauptstrahl im Meridianschnitt zur Verfügung. Man 

_ erhält die dazugehörige Beziehung, wenn man in (14) M’ = 0 setzt. Es ergibt sich dann 


100 (u, — w) = Am’ + © 04 D) m®+ (E-+2F)ms+ 2 Omt-+ Hms .. (19). 


Als weitere Bestimmungsstücke benutzen wir die Sagittal- und Tangentialschnittweiten die so- 
wohl auf dem Hauptstrahl als auch auf den einzelnen Komastrahlen trigonometrisch errechnet 

werden können, wobei man den jeweiligen Komastrahl als Leitstrahl auffaßt. Den Zusammenhang 
zwischen den sagittalen und tangentialen Schnittweiten und unseren Fehlergleichungen gewinnt 
Br man in folgender Weise: 


' Man bildet aus (14) und (15) 


). 


a). Br em HE+ me Gm + Hm RETET (20) 


role 


50 


und 


100 ER 


om. 


) — A+(30+2D)m +(3E+6F)m? + Soma 5Hm“ (21). 
De) i 


Aus Bild 2 folgt 
ww + Aw) tgwm _ 1 


Am’ -. (-+At)cosw 


=> Für Am’—0 ergibt sich daraus | ) ; 


- . x 


ee) a ea DERESEST * 
om’ /w=o cosu,\om/i=-0o +tcosw, 


und weiter wegen w’ — const 


(u, — w') _ 57,3 cos u, 
om’ FEN t 


20 


" Merıdianebene 
(M=0) 


(#+At)cosı u 


t« Schnittwerte meridionaler Nachbar- 2) 
strahlen auf dem Leitstrahl ER leitstrahl } (S+Af) :0s u, 3 | = 


? Bild 2 i Bild 3, Projektion auf die Ebene m’=0 
/ Schnittweite sagittaler Nachbarstrahlen auf dem Leitstrahl 
Aus Bild 3 folgt on RN N 
Bun ; | Atg u, N = a 


AM (f-+Af)cosw Ka Aw Bw 


und daraus für M’ = 0 
u Tre, I u, 
it — en RN ee 
| M’=0 . EN 3: 


om 2 [cos w; 

‘ Setzt man in (22) bzw. (23) ein 
{a8 | SDR Aura a 
t 1000’ ‚fr 210903. a un no Se 

dann folgt aus (20) und (22) SER; a Sale ee Re 
Ey: 5:7 DI HE a en a re ee 
j ag a ar br are .. (24) Eu 
und aus (21) und (23) R | ; N n. Kar Er 
5,73 Dt. cos u’ = 44 8042D)m’ + BE+6Mw2 Gm +5 Hm Be 0 

Man hat nun an Hand der Beziehungen (19), (24) und (25) eine mehrfache Auswahl An B 
mungsstücken für A,..., H. Man braucht dazu die sagittalen Schnittweiten auf ind 
drei Strahlen. Nach Belieben können aber auch vier oder fünf sagittale Schi 
werden. -Die fehlenden drei bis fünf Bestimmungsstücke kann man belie s de 


abweichungen im Hauptschnitt und den tangentialen Schnittweiten auf H wupts‘ I ahlb 
strahlen auswählen. Benutzt man die de von Yen: Komas 3 


ne Rechnungsgang: 
Gegeben sind 


w Neigung des Hauptstrahls, 
Ws... %, Neigung der Komastrahlen n E- ” 
MM, Austrittspupille der Koe sranlal D 


5 ange. at, en 1os1 Köhleru.Pradel, Eine neue Interpolationsmethode i al 


Dsp, Ds,, Ds, sagittale Schnittweiten auf dem Hauptrahl und den Komastrahlen 2 und 3, 
gemessen in Dioptrien von der Austrittspupille aus, 


Din - tangentiale Schnittweite auf dem Hauptstrahl, ebenfalls i in Dioptrien von der 
Austrittspupille aus gemessen. 


Daraus berechnet man der Reihe nach: 


A = 5,73 Dty cos w' 2 2 Den 
cos w' 
10 (w,—w) 4 N 
fi EEE W.=1,2,3,4) 
An —m, m, Pı=m, + m, 
Y Ya ont Pi 1 Ayı Par 
3 SENSE. e: „_mhomf, 
D‘ Mn ie < 4 ’ Ab ==, ’ ’ 
| m, — m, m — m, 
a a SE EN NR ee NE ns 
% j) m, fe m,S; 
BR re a 
ed Pe 7 B = Pa— Pr 
EN A = Ay As % 
__ AdarA— öx y_"ZA, 
 5ay— ßö  ay— ßÖ 
TR 5 3 5 
X. E+2F—erg— gßn@—YaH SOHDE As 4% @ — 6. H 
| 5 5 
ENT 4 Paı G—YaıH Fr — An 4a G—öyH. 


Hier liefert jeweils die Formel mit den Indizes 32 den genaueren Wert. Die en Aus- 
drücke können als Kontrolle benutzt werden. 


f, = — ai 1 B) — Om’? — Hm‘; Y 4, genau entsprechend 
4 Ms —Me Wan Er m, Hy — m, Us | 2-6 = 
Hrn nn. D=-(50+D2)-50 
an En a BSEHNE 


Abschließend sei hier noch bemerkt, daß es nur bei geeigneter Auswahl der in den Ansatz (3) 
aufzunehmenden Potenzen von m’ und M’ möglich ist, die Koeffizienten von (14) und (15) aus 
Größen nullter und erster Ordnung von Strahlen des Meridianschnitts zu bestimmen. 


Bi er r \ V. Berechnung der deformierten Austrittspupille 


- Man benutzt für die Berechnung der deformierten Austrittspupille ein Annäherungsver- 
u, fahren 3. Ordnung, das man u.a. bei M. Herzberger [1] findet. Bezeichnen m und M die 
Koordinaten der BRLSSSEnpIle, dann kann ‚man die Banzcnaten der Austrittspupille wie folgt 
2. “en Bir, 


m im +[' ie a +bamtem+e! sta. 


NER Baer . (26) 
ga wir Heumkanret ey 
uf + Hm ten + Be Mm 

ar (27) 
i ne +em+: M)+e 14 ss 

en ar " Be, . 4* 


2 


AT n: 


Br Ir 
ya) 
1, h 
N 

er 
A # 
x 
ER 
5 
B, 
Pr 
N 
W 
h 
EL 
er. 
a 

: 
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Betrachtet man wieder den En, allein interessierenden Fall z—= 0 und y= const, dann geht (26) 
und (27) über in 


ni (BEE) m m m a + Emm... 
und 


Mm —(B+%) M+Dymit+ 5 Mm + M9) EEE (29). 


Da es im folgenden nur auf den Zusammenhang zwischen m’ und M’ ankommt, setzt man 
m —= mR, wobei R den Radius der vollen Eintrittspupille bedeutet, und nimmt die Konstanten R 
und y in die vorläufig noch unbestimmten Koeffizienten auf, Somit lauten die endgültigen 
Bestimmungsgleichungen für die Austrittspupille 


| m am tem dm HP) tem) | - > - > - - @0) 
und 


M'=bM+cmM+eMm®+M)|---:...:....(). 


Diese Beziehungen gelten streng zunächst nur für die Gaußsche Austrittspupille. Bekannter- 
maßen unterscheidet sich die Lage der Austrittspupille von der der Gaußschen bei den hier 
interessierenden Systemen immerhin um Zehntel der Okularbrennweite. Man kann jedoch zeigen, 
daß die allgemeine Form der Gleichungen (30) und (31) auch für die wirkliche Austrittspupille 
eine Näherung darstellt, wenn man f= 0 setzt. Auf die Wiedergabe dieser Untersuchungen wird 


hier verzichtet. Bei der Bestimmung der fünf Konstanten a, b, c, d, egeht man in ähnlicher Weise er 
wie im vorhergehenden Abschnitt zu Werke. Als Des DOIDEupE ES stehen aus der Durch- 
rechnung im Hauptschnitt zur Verfügung: 

1. sagittale und meridionale Pupillenvergrößerung auf dem Hauptstrahl; 

2. sagittale Pupillenvergrößerung auf mindestens zwei Komastrahlen, 

3. wenigstens zwei Werte der Austrittspupille im Hauptschnitt selbst. 

Um Formeln zu erhalten, die zu dem auf S.8 beschriebenen Rechnungsgang passen, be- 
nutzen wir m,, m, m,, m, und vg, %, ds, die sagittalen Vergrößerungsverhältnisse zwischen 
Austritts- und Eintrittspupille auf dem Hauptstrahl und den Komastrahlen 2 und 3. 

Bezeichnen BR f ‘x die auf die Fläche x des Systems bezogene sagittale Schnittweite vor 
und nach der Brechung an der Fläche, / gr die auf die Eintrittspupille bezogene Anfangsschnitt- 
weite, J ‚ap die auf die Austrittspupille bezogene Endschnittweite, dann ist für den jeweils be- 
trachteten Strahl bei einem %-flächigen System 

_ Jr. Js... Ser, FR 
Ser Se } 
Für die Koeffizienten a, ..., e erhält man bei x 
m —1 My = =L N Kir m—=—1 | 
1 ’ 2 y2 ’ 3 y2 ’ 4 2 
die Formeln h2 

b=ü%yR Zi e=(%-+ v;) R—2b 

1 / 1 BRER \? 

a= Im, — m; +2 (m, — m;)—3e] Im tm tm tHm)—e| 


Dabei sind b, c, e allein aus den sagittalen Vergrößerungsverhältnissen bestimmt. aundd 
sind mittels der Methode der kleinsten Quadrate so ausgeglichen, den My m, m, m, mit IORnE 
lichst kleinen relativen Fehlern wiedergegeben werden. Brit 

Es darf nicht verschwiegen werden, daß in manchen Fällen die Formeln (32) nicht aus- ve 
reichen. Die nötige Genauigkeit läßt sich dann wahrscheinlich mit anderen Verfahren erreichen. 
Untersuchungen darüber sind im Gange. - zer 32 Kost Cy er: 


VI. Die Berechnung der Zerstreuungsfigur & 


Er Set > 
Nachdem im Abschnitt IV die Konstanten A, B,C,D,E,F,G,H ermittelt En ar 'her 
gehenden Abschnitt die Koordinaten m’ und M’ bzw. p' aus den ee Werten d 
| hi 
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Eintrittspupille mittels der Formeln (16), (17) 
berechnen. | 
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Eintrittspupille bestimmt worden sind, lassen sich nunmehr für jeden beliebigen Punkt der 


und (18) die resultierenden Winkelabweichungen 


VI. Schlußbemerkung 


Variable auf den rechten Seiten die Koordinaten der Eintrittspupille einführen würde. Die Praxis 


| Die Beziehungen (16) bis (18) würden theoretisch befriedigender, wenn man als unabhängige 
j 


bessere Annäherung liefern. 


hat jedoch bisher gezeigt, daß die auf die Austrittspupille bezogenen Reihenentwicklungen eine 


j Literatur 

| [1] H. Köhler: Ein neues Verfahren zur quantitativen Beurteilung der Bildgüte optischer Instrumente zur 
j subjektiven Beobachtung. (Veröffentlichung demnächst geplant.) 

| [2] M. Herzberger: Theory of image errors of the fifth order in rotationally symmetrical systems. J. opt. 


Soc. Amer. 23 (1939), S. 395. 


ad 
. 


Mech, 30 (1950) 8. 226. 
Eingegangen am 3. Dezember 1949. 


7 Orthogonalisierung (Unitarisierung) von Ma- 
trizen nach E. Schmidt!) und ihre prak- 
tische Durchführung. 

Im n-dimensionalen Raum können n gegebene linear 
unabhängige reelle Vektoren a, orthogonalisiert und 
normiert werden. Bilden die gegebenen Vektoren die 
Zeilen einer quadratischen Matrix VW, so wird diese 
durch den Orthogonalisierungsprozeß in eine Matrix 
mit im allgemeinen nicht normierten, also nur zeilen- 
weise orthogonalen, Vektoren transformiert. Werden 
j diese Vektoren noch normiert, dann erhält man eine 
H orthogonale Matrix 9 mit 8 = VD =E. 

Handelt es sich um komplexe ‘Vektoren, dann gelten 
die Aussagen immer noch, wenn an Stelle der Ortho- 
gonalisierung die Unitarisierung tritt. Die komplexe 
Matrix X wird dann in eine unitäre Matrix transfor- 


miert. 

In beiden Fällen wird die Matrix X mit einer Drei- 
ecksmatrix vormultipliziert, wobei im komplexen Falle 
h auch komplexe Koeffizienten auftreten. Die Ermitt- 
i lung dieser Koeffizienten ist im allgemeinen recht un- 
übersichtlich. Im folgenden soll gezeigt werden, daß 
r bei der Orthogonalisierung bzw. Unitarisierung nur 
das Gleichungssystem AW = € bzw. WAWX=€ 
betrachtet zu werden braucht. Führt man das ent- 


Algorithmus in ein gestaffeltes Gleichungssystem über, 
i dann erhält man direkt die gesuchten Koeffizienten. 
Bekanntlich kann die Auflösung linearer Gleichungs- 
systeme auch durch Orthogonalisierung (Unitarisie- 
rung) vorgenommen werden, was hier ebenfalls kurz 
besprochen werden soll. 

Zunächst wird nur mit reellen Vektoren gearbeitet, 

l. Orthogonalisierung von Vekto- 
ren. Um die Matrix XV mit den Zeilenvektoren a; zu 
orthogonalisieren, macht man bekanntlich den Ansatz: 


“—1 
zug 2 Bi 
= Aus der Forderung, daß r, zu t,, ty ...,%;_, 

-  gona] sein soll, also 
en -0 für k=12% ..,i-1. 
erhält man zur Bestimmung der Koeffizienten die Glei- 


Er ei u + fü i=1,2...,r (). 


ortho- 


3). 


mid t: Rend.Cire. Mat. Palermo 25 (1908); A. Wint- 
theorie der unendlichen Matrizen, Leipzig 1929, $6. 


sprechende Gleichungssystem mit dem Gaußschen .. 


I 


[3] H. Köhler: Ein einfaches Verfahren zur Ermittlung der gesamten Zerstreuungsfigur optischer Instrumente 
aus der meridionalen Durchrechnung auf Grund der Bildfehlertheorie dritter Ordnung. Z. angew. Math. 
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DieGleichung (1) kann man auch in der Form schreiben 


x, a, Pb 01:04:7,,050 
= fi Ra A 
2 = mit = a (4), 
: L sa: IN 0 
In 17} Ini In2 In,n-ı 1 
woraus folgt: 
a\ 
x i x 
. 2 
:|n=m. (5). 


(5) 


Die Dreiecksmatrix, mit: der X vormultipliziert werden 
muß, um eine zeilenweise orthogonale Matrix R zu er- 
halten, ist also 2-1, 


Führen wir noch die Diagonalmatrix ® ein, diein der 
Diagonale die Beträge-der orthogonalen Vektoren ent- 


hält, 
B Ere 


®-|: - | mit 4,- + 09), 
N 


1) 


so können die Vektoren in # normiert werden und man 
erhält schließlich eine orthogonale Matrix ©: 
DAILIY-DIR=OD 


2. Auflösung linearer Gleichungs« 
systeme?). Auch zur Auflösung eines linearen 
Gleichungssystems 


kann man das Orthogonalisierungsverfahren heran- 
ziehen. Man multipliziert (8) mit D-1 2-1 vor und er- 
hält ein Gleichungssystem mit einer orthogonalen 


Matrix: ‚ 
Dr=-D181H9 (9) 


Multipliziert man (9) mit DO’ vor, so erhält man unter 
Berücksichtigung von 9 DO =€:) 


t=- VID Hy =D DIT. . (10). 


a FR et hr, 


2) Vgl.E. Bode wig: Bericht über die verschiedenen Methoden 
zur Lösung eines Systems linearer Gleichungen mit reellen Koeffi- 
“zienten. Proc. Akad. Wet. Amsterdam 1947, insb. 8,1105°—1110. 


Sr, 
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Will man die Kehrmatrix berechnen, dann ist für 9 
in (10) die Einheitsmatrix zu setzen: 


AA-WEADELI=-HDIII.... 


Die Auflösung von Gleichungssystemen durch Ortho- 
gonalisierung kann allerdings wie schon E. Bode- 
wig?) festgestellt hat, praktisch kaum empfohlen 
weıden. Denn gegenüber dem G a ußschen Algorith- 
mus liegt die Zahl der Rechengänge wesentlich höher. 


3.Bestimmung der Matrizen & und 1, 


a) Man kann zunächst die li; in bekannter Weise 
unter Benutzung von (3) bestimmen. Die Matrix R 
kann dann nach (4) Zeile für Zeile berechnet werden. 
2-1 wird als Kehrmatrix einer Dreiecksmatrix entweder 
durch das Produkt 


E.0,2.28205, ryl%. 0 
01. o 0o\/o. 
g-1 - 
1.8200 
000 0) a 0 
0 0 En — Inn-1 1 0% x Rue In, n-2 0 1 (13) 
0 0 
el 0) 


ln 0a 8 > 1) 


oder übersichtlicher durch Aufrechnen des gestaffelten 
Gleichurgssystems 2X = € bestimmt. 

b) Eine übersichtlichere Anordnung, wenn auch unter 
einer geringfügigen Erhöhung der Rechengänge, erhält 
man durch folgende Betrachtung. Aus (11) folgt durch 
Vormultiplikation mit W’-1: 


AW)1=L19224 AV=-LDW.. a3. 
Diese Beziehungen legen den Gedanken nahe, die 
Matrix Adurch Heı anziehungdesGleichungssystems mit 


der symmetrischen Matrix A W’ zu orthogonalisieren. 
Führt man das Gleichungssystem 


ARE SER. eh (14) 


mit dem abgekürzten Verfahren von Ga u ß°))5) in ein 
gestaffeltessüber, spaltet also die Matrix AW’ in die 
zwei Dreiecksmatrizen 2 und D° 2’ auf und multipli- 
ziert mit: 2-1 vor, so erhält man auf der rechten Seite 
eine zeilenweise orthogonale Matrix. 


VYH-MA-R ....(B) 


Der Zusammenhang mit dem Gleichungssystem wird 
weiter nicht benötigt, da das. Rückwärtseinsetzen — 
das Aufrechnen des gestaffelten Systems — unter- 
bleibt. Die Matrix &, wird nicht bestimmt. 

In der Hauptdiagonale der Dreiecksmatrix D? 2%’ 
stehen die Normen der orthogonalen Vektoren. 

Die gesuchte orthogonale Matrix erhält man durch 
Vormultiplikation mit D-1, 


Din ne, 


Verwendet man das abgekürzte Verfahren von Gauß 
in der Anordnung von Cholesky ?) 5), so erhält man 


aus (14): 
DV/’H=- DILL AS IIFERAN, 


womit sofort die normierte orthogonale Matrix erhalten 
wird, 


®)E.Andersen, Solution of great Systems of Normal Equa- 
tions together with an Investigation of Andrae’s Dot-Figure. Geod. 
Inst. Skr. (M&m. Inst. geod. Danemark) (3) 11 (1947). 

*) Benoit: Sur une möthode de r&solution des &quations nor- 
males provenant de l’application de la möthode de moindres carr6s 
ä un systeme d’&quations lineaires en nombre inferieur ä celui des 
inconnues. — Application de la methode ä la resolution d’un sy- 
steme defini d’@quations lineaires (Proced& du Commandant Cho- 
lesky). Bull. g&od. 1924, Nr.2, p.67—77; vgl. auch: Tryggve- 
Rubin, Ett nytt sätt att lösa normalekvationer. Svensk Lant 
mäterj-Tidskr. Lidköping 1926, Häft1, 3—9. 

°) R, Zurmühl: Matrizen, Berlin 1950, $8.248—262, 
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‘Will man die Auflösung von Gleichungssystemen 
durch Orthogonalisierung vornehmen, so geht man von 
dem. Gleichungssystem 


(U) Er re (18) 


aus und erhält dann nach der Überführung in ein 
gestaffeltes System auf der rechten Seite zunächst-2-1 
(biCholesky®-1272). Man bildet dann 2-2 A=R. 
Die Lösung des Gleichungssystems nach (10) bzw. die 
Bestimmung der Kehrmatrix nach (11) erfordert nur 
noch weitere Matrizenmultiplikationen. 


4, Unitarisierung von Matrizen mit 
komplexen Elementen, 

Durch die komplexen Vektoren a; ergeben sich im 
allgemeinen komplexe Vektoren rt; und komplexe Ko- 
effizienten I,,. Die Bedingungen (2) und (3) lauten 
jetzt: Be 

el ae ae 

Ik 7 he TI a I 

ı fa 15T 
(Querstriche bedeuten konjugiert komplexe Werte.) 
Weiter sind in der Diagonalmatrix (6) die positiven 
Quadratwurzeln der Normen t,r; der komplexen Vek- 
toren einzusetzen. - 

D1214Y=MU ist dann eine unitäre Matrix mit 
DU = : 

Bei der Behandlung eines Gleichungssystems erhält 
man an Stelle von (10) 


Er FDA Ty2R 

während sich für die‘ Kehrmatrix ergibt: 
KA=-NETISeR 1, ;... (11a). 
Aus (11a) folgt durch Vormultiplikation mit W!: 
(AW)I=- FD bzw. AU -LDV (13a). 
Im Falle einer komplexen Matrix A empfiehlt es sich 

also, das Gleichungssystem 
(AU)L-A bw. (AW)R,=- € 


mit der Hermiteschen Matrix AW-zu behandeln. 
Führt man durch Anwendung des Ga u ßschen Alge 
rithmus das Gleichungssystem in ein gestaffeltes über, 


. 


. '.(108), 


‚so erhält man die gesuchte komplexe Matrix 2-12 A 


bzw. die Dreiecksmatrix 2-1, 


5. Beispiel zur Orthogonalisierung 


314 26 18 11 
Die Matrix A = (& 2 1) mit AV = (15 21 1) 

Ze 11.2138 
ist zu orthogonalisieren. Mit dem Verfahren von Cho- 
lesky°) erhält man durch Überführung des Glei- 
chungssystems AW’ X = Win ein gestaffeltes: DL’X 
—-DI8R1YA=D. (In dem Schema steht rechts vom 
Doppelstrich die orthogonalisierte Matrix. 


Darmstadt 


Ein Nomogramm erleichtert die Bestimmung 


eines Funktionsverlaufes[. 


In seinem Buch ‚‚Notions &l&mentaires de Math&“ ‘a 
matiques“ [2] zeigt L. Brillouin folgende Formel, 


deren ungefährer Verlauf dort nur angedeutet inte) 


ER GEN | | 
BI N 100 2 I . . et ve 
RN wo 7 n Ge 2 » v& cz 
PN 


Z.ang6w.Math.Mech. 


} 
e 
i 
9 


Z.angew. Math, Mech, 
 Bd.31 Nr.1/2 Jan./Febr. 1951 


Um den Verlauf dieser an und für sich einfachen For- 
mel zeichnerisch darzustellen, bedarf es umfangrei- 
cher Rechnungen, selbst dann, wenn die Funktion e 
bis zu # = 1000 tabuliert vorliegen sollte. Auch die 
logarithmische Rechnung erfordert einen großen Zah- 
len- und Rechenaufwand. Selbst die Zuhilfenahme 
4 der Höheren Mathematik entbindet uns nicht von den 
E sich anschließenden Berechnungen. 

Hier erweist sich nun die Anwendung nomogra- 
{ Phischer Methoden als wertvolle Hilfe, wie im folgen- 
den gezeigt wird. Zunächst logarithmieren wir (1) 
und erhalten 


lgy=s-lge— MWlge. (2) 
Durch weitere Umformung erhalten wir schließlich 


ı_ les + 100 1g x 


0,4343 © RE Tg 


j, Bild 1, 


Nach Schwerdt [3] können wir (3) der folgenden 
Grundform zuordnen: 


% F(B) + Fa) 

3 j F = 207 nn ,,..(4 
E.:: IF NUERERT) 

2 Das danach entwickelte Nomogramm zeigt Bild 1. 


Der Verlauf der Funktion kann nunmehr unter 
Zwischenschaltung beliebig vieler Werte leicht be- 


Bei 


er 7000 


— 


A END EL ERNRTEN. Dr U 
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stimmt und ohne große Mühe sofort gezeichnet wer- 
den (Bild 2). 

Das Minimum wird aus dem Nomogramm durch 
Ziehen der Tangente an die Kurve für x gefunden, 
es liegt für z = 100 bei y = 10-156,6, 

Da die Ablesegerade immer durch den Punkt P 
gehen muß, der also ein Drehpunkt ist, lag es nahe, 
die Funktion in Kreisform darzustellen. Begrenzt 
man die y-Leiter auf die Werte von 10”160 pjs 10+200 
und läßt diese Strecke zur Peripherie eines Kreises 
werden, dann erhält das Nomogramm folgendes Aus- 
sehen ( Bild 3). 


r/) 700 


Bild 3. 


Weitere aufschlußreiche Zusammenhänge enthüllt 


die Funktion aber erst, wenn man den Exponenten u 
ebenfalls variabel macht. Statt (3) erhalten wir nach 


Umformung 
Igy+u-lgx 4 
1= BET Sy RE (5). 

-150 

70 

50 
> et 
on” 
u 

I 

oF—7 
s nn” 
y N 0° 


Se 


Bild 4, 


Als Bereiche für u wurden die Werte zwischen + 100 
und — 100 vom Verfasser als ausreichend angesehen. 
Um das Nomogramm anschaulich darzustellen, wur- 


a Ze 


; 
f 
k 
w 
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den die Parameter 1= — 1/100 und m= +1 ge- stelligen Logarithmentafeln errechnen kann, dann ist 


setzt [4] Das für die allgemeine Funktion 


konstruierte Nomogramm zeigt Abb. 4. Die Über- 

sichtlichkeit des Nomogramms erlaubt es, die Lage be- 

sonderer Werte sofort abzulesen: 

1. Alle Werte x, bei denen u = 0 ist, liegen auf der 
Verbindungsgeraden vom Koordinatennullpunkt 

zum Werte y= 1. 

2. Der Drehpunkt P mit y = 1 verbunden liefert alle 
x- und u-Werte, bei denen sich als Ergebnis der 
Wert y = 1 ablesen läßt. 

3. 2 = u sind Tangentenpunkte der von P ausgehen- 
den Ablesegeraden und zugleich 

4. die Minimumwerte für den Verlauf der e-Funktion 
bei bestimmten «. 

Interessant ist der Verlauf der y-Werte bei #» = win 

der Nähe von 1. Man findet bestätigt, was man durch 

Einsetzen in die Formel natürlich ebenfalls findet: 


Füru=0 wirdy=]1 
u) ms ME 
er ulN Bildr53} 


Bemerkenswert ist bei diesem Nomogramm, daß man 
recht schnell einen Überblick über die zu erwartende 
Größenordnung des Ergebnisses bekommt. Wenn 
man sich vergegenwärtigt, daß man dies sonst nur 
unter Zuhilfenahme von Rechenmaschinen und viel- 


der Vorteil der Anwendung eines Nomogramms in 
diesem Falle besonders zu erkennen. 


Bild 5 
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F. Mühlig, Der 24m-Interferenzkomparator 
des Geodätischen Institutes in Potsdam. 
(Deutsche Akademie der Wissenschaften zu Berlin. 
Veröffentlichungen des Geodätischen Institutes in 
Potsdam, Nr. 2.) 508. mit 30 Abb. Berlin 1949. 
Akademie-Verlag. Preis brosch. 7,— DM. 

Da der direkten Interferenzmessung nur kurze 
Strecken bis höchstens 25 cm zugänglich sind, stellen 
sich ihrer Anwendung in der Geodäsie zur Zeit noch 
erhebliche Schwierigkeiten entgegen. Immerhin gibt 
es doch bereits einige auch zur Messung längerer 
Strecken geeignete mit Lichtinterferenzen arbeitende 
Apparate: z.B. den eine Weiterentwicklung des 
Fabry-Perotschen Interferometers vorstellenden 
Apparat der Japaner Watanabe und Imaizumi 
mit einem Meßbereich bis 500 m, die Interferenz- 
Kömparatoren für Invar-Basisdrähte von Väisälä in 
Finnland und von Warlich und Schwarz in der 
Sowjetunion. Einen weiteren sehr wesentlichen Fort- 
schritt in dieser Richtung bedeutet der neuerdings 
vom Geodätischen Institut in Potsdam entwickelte 
24 m-Interferenzkomparator. 

Bei der Konstruktion dieses neuen Komparators 
ist man besonders darauf bedacht gewesen es so einzu- 
richten, daß alle zu einer Eichung gehörenden Ope- 
rationen in möglichst rascher Aufeinanderfolge ge- 
schehen können, damit der Apparat auch unter feld- 
mäßigen Bedingungen verwendbar sei. Grundsätz- 
lich hat der Apparat von Watanabe und 
Imaizumi als Vorbild gedient; geringfügige 
konstruktive Abänderungen gestatten aber auch 
das Arbeiten nach dem Prinzip von Väisälä, 
Der bereits angefertigte Apparat ist lediglich ein 
noch mit mancherlei technischen Mängeln behaf- 
tetes Versuchsmodell; die damit gemachten Erfah- 
rungen sollen bei der endgültigen Konstruktion 
nutzbringend verwertet werden, 

An die Beschreibung der Einrichtung und Hand- 
habung des Komparators und ausführliche Angaben 


über die ausgeführten Versuchsmessungen reihen sich 
in besonderen Kapiteln gebrachte wertvolle Abhand- 
lungen über mehrere Einzelprobleme der Interferenz- 
messung: Theorie des Meßkeils, Einfluß der Ab- 
weichung vom Parallelismus in der Einstellung der 
Etalonplatten, Einfluß einer prismatischen Gestalt 
der mittleren, beiden Etalons angehörigen Platte, 
Theorie des Vergleichs des Ausgangsetalons mit der 
12,5 cm-Strecke, Einfluß der Abblendung eines Teil- 
bündels, Einfluß der Dispersion des Phasensprungs 
auf die Messungen. Den Abschluß bildet eine kurze 
Darstellung der wichtigsten Interferenzerscheinungen, 

Alles in allem genommen: eine als außerordentlich 
wertvolle Bereicherung der geodätischen Literatur 
einzuschätzende Veröffentlichung. 


Dresden, A. Buchholtz. 


. 
Dr, Walter Weizel (o. Professor der Physik an der 
Universität Bonn), Lehrbuch der Theore- 
tischen Physik. Erster Band: Physik der 


Vorgänge. Bewegung, Elektrizität, L’cht, Wärme. . 


XV + 771 8. mit: 270 Abb. Berlin/Göttingen/Heidel- 
berg 1949. Springer-Verlag. Preis geh. 53,— DM, 
geb. 56.90 DM. 


Das neue Lehrbuch der Physik umfaßt zwei Bände, 


von denen der erste, jetzt vorliegende, den Untertitel 
„Physik der Vorgänge“ trägt und die phänomenolo- 
gische Theorie enthält, während der zweite Teil die 
Atomtheorie behandeln soll. 


Inhalt des ersten Bandes: Mechanik der Massen- 
punkte und starren Körper, Mechanik der Kontinua, _ 


Elektrodynamik, Optik, Elektrodynamik bewegter 
Körper, Thermodynamik. Der Verfasser sagt zwarim 
Vorwort, daß sein Werk nicht den Anspruch erhebt, 


völlig neue Wege bei der Darstellung des Lehrstoffes 
zu beschreiten. Dies mag im Gesamtbild zutreffen, 
der Leser wird jedach feststellen, daßzahlreichecharak- 
teristische Einzelzüge das Buch vor anderen auszeich- > 
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nen. Der Inhalt geht in vielen Punkten sehr weit über 
das übliche Maß der kurzgefaßten Kompendien hinaus, 
ohne daß ein ungebührlich großer Umfang von vorn- 
herein abschreckt. Die Darstellung ist so ausführlich, 
daß ein Student der mittleren Semester leicht folgen 
kann. Durch kurze Inhaltsangaben und eine besondere 
Zusammenstellung der jeweils neu eingeführten Größen 
wird das Durcharbeiten der einzelnen Abschnitte er- 
leichtert. Alle Gleichungen sind im Giorgischen 
Maßsystem geschrieben; dadurch wird dem Elektro- 
techniker die Darstellungsweise vertrauter. Man darf 
dem Verfasser zu dieser ausgezeichneten Arbeit Glück 
wünschen und für das Buch eine verdiente, große Ver- 
breitung erwarten, 


Dresden, A. Recknagel. 


v 

K. Marguerre (o. Prof. a.d. Techn. Hochschule Darm- 
stadt), Neuere Festigkeitsprobleme 
desIngenieurs. Ausgewählte Kapitel aus der 
Elastomechanik. Von Prof. Dr.-Ing. W. Flügge, 
Stanford USA, Prof. Dr.-Ing. R. Grammel, 
Stuttgart, Prof. Dr.-Ing. K. Klotter, Karlsruhe, 
Prof. Dr.-Ing K.Marguerre, Darmstadt, Prof. 
Dr. G. Mesmer, Darmstadt. Herausgegeben von 
K.Marguerre. VIIL+ 253 S. mit 120 Abb. Ber- 
lin-Göttingen-Heidelberg 1950. Springer-Verlag. 
Preis geb. 25,50 DM. 

Anläßlich einer Vortragsreihe, die vom damaligen 
Vorsitzenden des VDE-Berlin angeregt worden war 


. und im Winter 1941 vor wissenschaftlich interessierten 


Ingenieuren der großen Berliner Betriebe gehalten 
wurd£, entstand eine zusammenfassende Darstellung 
des seinerzeitigen Standes der wissenschaftlichen 
Festigkeitslehre, soweit sie vorwiegend für die tech- 
nische Praxis von unmittelbarem Interesse ist. Das 
‚vorliegende Buch, das eine Zusammenstellung der 
Manuskripte jener Vorträge darstellt, gibt zugleich 
‘einen Einblick in die bis zum Jahre 1941 erfolgten 
Weiterentwicklungen der Festigkeitsforschung, und 
zwar auf folgenden Teilgebieten: Abschn. I erstreckt 
sich auf die experimentellen Methoden der Spannungs- 
bestimmung (Dehnungsmessung, Reißlackverfahren, 
Anbohrverfahren, Röntgenverfahren, Spannungsop- 
tisches Verfahren). Die Abschn. II und III be- 
ziehen sich auf Spannungsverteilungsprobleme der 
Elastizitätslehre. Hierbei wird nach einem kur- 
zen Rückblick auf Stab und Platte die Schale 
als ein “neuerdings wichtig gewordenes Bauglied 
ausführlicher dargestellt. Abschn. IV bis VI han- 
deln von den Verfahren zur Lösung technischer 
Eigenwertprobleme. Zunächst werden Schwingungs- 
.erscheinungen im Bau- und Maschinenwesen be- 
handelt, dann Drehschwingungen von Reihen- 
motoren, Kritische Drehzahlen, Torsionsschwingungen, 
Biege- und Flatterschwingungen. Schließlich werden 
Knick- ‚Beul- und Durchschlagprobleme für Stab und 
Platte verhältnismäßig ausführlich behandelt. Da die 
Verfasser der einzelnen Kapitel zu den führenden 
Fachleuten auf den entsprechenden Spezialgebieten 
ehören, vermittelt das Buch einen sachlich sehr guten 
Überblick über den derzeitigen Stand der Festigkeits 
forschung. 


Dresden. H. Neuber. 


Ernst Hölder, Überdie Variationsprin- 
zipeder Mechanikder Kontinua. (Ber. 
ü. d. Verh. d. Sächs. Akad. d. Wiss., math.-nat. 
Klasse. Bd. 97, H.2); 13 S. Berlin 1950. Akademie- 
Verlag. Preis brosch, 2,50 DM. 

Verfasser benutzt das Hamilton’sche Prinzip zur 
Formulierung eines Variationsproblems für das all- 
gemeine Kontinuum, mit Anwendung auf die statio- 


_ näre Strömung eines Gases im Raum. Die Arbeit 
stellt eine Verallgemeinerung des Variationsprinzips 
von ÖÜlebsch dar [Über eine allgemeine Transfor- 
_ mation der hydrodynamischen Gleichung, J. f. reine 


F 


und angew. Mathematik 54 (1857), S. 293-312, 
ferner 56 (1859)]. Für den ebenen und rotations- 
symmetrischen Fall hatte Verf. das entsprechende 
Variationsprinzip auf der DMV-Tagung in Jena 1941 
bzw. in Bd.1 Seite 449—450 der meteorologischen 
Rundschau 1948 mitgeteilt. 


Dresden. H. Neuber, 


Prof, Dr. Friedrich Tölke, Mechanikdefor- 
mierbarerKörper. Erster Band: Der Punkt- 
förmige Körper. VIII + 388 S. m. 339 Abb. Berlin- 
Göttingen-Heidelberg 1949. Springer-Verlag. Preis 
geb. 45,— DM. 

Der Verf. hat sich die gewagte Aufgabe gestellt, das 
bereits umfangreiche Schrifttum auf dem Gebiete 
der Mechanik .durch eine neue Buchreihe zu ergänzen, 
welche der Mechanik der deformierbaren Körper ge- 
widmet ist und damit eigentlich fast das gesamte 
Gebiet der Mechanik überhaupt umfassen müßte. Das 
geplante Werk sollsich vorwiegend an die technischen 
Anwendungen der Mechanik anpassen und vor allem 
die durch die Weiterentwicklung der Technik (und 
zwar vornehmlich des Maschinenbaues) bedingten 
Hauptprobleme behandeln. Der Verf. will sich dabei 
an folgende Aufteilung halten: Bd.I: der Punkt- 
förmige Körper; Bd.II: der statisch beanspruchte 
feste Körper; Bd. III: der dynamisch beanspruchte 
feste Körper; Bd. IV: der thermisch beanspruchte 
feste Körper; Bd.V: Flüssigkeiten und Gase. Es 
liegt zunächst der 1. Band vor, welcher sich trotz des 
einschränkenden Titels nicht nur auf den Massen- 
punkt, sondern bereits auch auf materielle Systeme 
(Punkthaufen) erstreckt und in folgende Kapitel glie- 
dert: Der geradlinig bewegte, punktförmig idealisierte 
Körper; vektorielle, geometrische und kinematische 
Grundlagen; mechanische Grundlagen; Bewegungen 
in zentralen Potentialfeldern; Mechanik der Raum- 
und Relativbewegungen; Massenmittelpunkt des 
Haufensystems; Mechanik des Haufensystems; die 
gekoppelten harmonischen Schwingungen in Verbin- 
dung mit erzwungenen Schwingungen; die gedämpf- 
ten Schwingungen. Neben den ausführlich darge- 
stellten allgemeinen Grundlagen der Mechanik, wie 
z.B. den Newtonschen Axiomen und den zuge- 
hörigen kinematischen Definitionen, bringt Verf. ver- 
hältnismäßig ausführlich eine Zusammenstellung 
mathematischer Grundlagen, insbesondere der Vek- 
torreehnung. Dadurch, wie auch in technischer Hin- 
sicht, entsteht eine Überschneidung des Buchinhaltes 
mit anderen Werken der Fachliteratur. Dieser Ein- 
druck wird durch die zu abstrakte Wahl des Buch- 
titels, sowie der Kapitelüberschriften noch verstärkt. 
Die ziemlich umfangreiche Behandlung der Punkt- 
haufendynamik läßt dieEinführungder Lagrange- 
schen Gleichungen und der Hamiltonschen 
Theorie vermissen, welche die Darstellung vielfach 
wesentlich verkürzt und erleichtert hätten. Das Buch 
enthält jedoch manche technisch interessanten Bei 
spiele (insgesamt 55), welche den lernenden Leser 
anregen und zum Teil auch dem Fachmann Neues 
bieten werden. Dennoch bleibt die Frage sehr um- 
stritten, ob bei der "bereits bestehenden umfang- 
reichen Fachliteratur der Druck eines so umfang- 
reichen Werkes, das jainfolge der hohen Anschaffungs- 
kosten kaum einem großen Leserkreis zugänglich sein 
wird, wirklich vertretbar ist, oder ob die Abfassung 


“derartiger Werke, welche umfangmäßig die Rolle von 


Standardwerken zu übernehmen hätten, in Zukunft 
nicht besser als Gemeinschaftsarbeit mehreren Fach- 
gelehrten übertragen werden müßte, welche den In- 
halt den Bedürfnissen des Ingenieurs in den verschie- 
denen Richtungen anzupassen in der Lage sind, so 
daß sich dann wenigstens alle einschlägigen Biblio- 
theken zur Anschaffung entschließen könnten. 
Schließlich hätte auch eine geschicktere Wahl des 
Buchtitels den Absatz erleichtert. 


Dresden. H. Neuber. 
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Dr, Heinrich Blasius, Mechanik, Physika- 
lische Grundlagen vom technischen Stand- 
punkt. Dritter Teil: Kinematik, Dynamik, Hydrau- 
lik. Dritte neubearbeitete Auflage. 351 S. mit 255 
Figuren und 111 Aufgaben. Hamburg 1950. Boysen- 
Maasch-Verlag. Preis: kart. 14,— DM, geb.17,—DM, 


Nachdem in dieser Zeitschrift bereits auf die Neu- 
auflagen des 1. und 2. Teiles der Blasiusschen Me- 
chanik hingewiesen wurde, [Bd.29 (1949) S.128; 
Bd. 30 (1950) S. 127] soll nunmehr auch die kürzlich 
erschienene dritte Auflage des 3. Teiles besprochen 
werden. Die 1. Auflage des 3. Teils, die bereits 1935 
erschien, und damals noch nicht ausführlich bespro- 
chen wurde, enthält, wie die ersten beiden Teile, vor- 
wiegend eine Reihe ingenieurmäßiger Aufgaben, und 
zwar aus dem Gebiete der Kinematik, Dynamik und 
Hydraulik. Der Verf. bevorzugt‘ auch hier das Ver- 
fahren, beim Leser an Hand praktischer Beispiele die 
Erkenntnismethodik und damit die Umdeutung des 
technischen Sachverhaltes in einen der Mechanik zu- 
gänglichen Modellvorgang zu schulen, wenn auch für 
die rein logisch und abstrakt-mathematisch durchzu- 
führende weitere Behandlung kein höherer Standpunkt 
eingenommen wird (das Buch wendet sich vorwiegend 
an den nichtakademischen Ingenieur). Dadurch ge- 
winnt der Leser eine gewisse Lösungstechnik, die durch 
die Mannigfaltigkeit des dargebotenen Beispielma- 
terials gefördert wird. Das Herauskommen der 3. Auf- 
lage, die gegenüber der 1. und 2. lediglich hinsichtlich 
des Wortlautes noch einmal überarbeitet ist, wobei 
gleichzeitig kleinere : inhaltliche Fehler ausgemerzt 
wurden, läßt erkennen, daß das Buch bereits Anklang 
gefunden hat. 


Dresden. H. Neuber. 


Adalbert Duschek, Vorlesungen über hö- 
here Mathematik, I.Band. X- 3958. 
mit 167 Abb. Wien 1949. Springer-Verlag. Preis 
$ 7,80, geb. $ 8,70. 

Wenn der Verfasser hier den ersten Band einer auf 
vier Bände berechneten Darstellung der Mathematik 
in dem Umfange, wie Techniker und Physiker sie für 
ihre wissenschaftliche Arbeit bedürfen, vorlegt, so 
glaubt er, eine gewisse Rechtfertigung schuldig zu 
sein. Er sieht sie darin, daß die Anforderungen, die 
die moderne Teghnik an den Ingenieur in den For- 
schungsabteilungen und in den Berechnungsbüros 
stellt, in den letzten Jahren gewaltig gestiegen sind 
und immer mehr steigen werden, und daß der wissen- 
schaftlich arbeitende Techniker ein guter Mathe- 
matiker sein muß, wenn er diesen Anforderungen 
entsprechen will. Die Wiener Technische Hochschule 
habe inihrer neuen Studienordnung diesem Umstande 
Rechnung getragen, da die Industrie vielfach eher auf 
praktische Ausbildung verzichten könne als auf die 
theoretische, und sehe eine starke Vermehrung der 
mathematischen Vorlesungs- und Übungsstunden vor. 
Wir können dem Verfasser dankbar sein,. daß er sich 
durch diesen äußeren Anlaß bestimmen ließ, seine 
Vorlesungen in Buchform herauszugeben. Mit allem 
Nachdruck bekennt er sich zu der Forderung, daß 
dem Naturwissenschaftler und Techniker zwar eine 
nach seinen Bedürfnissen ausgewählte Mathematik, 
diese aber mit der nötigen Strenge ohne Umgehung 
grundsätzlicher Schwierigkeiten geboten werden muß, 
So erstrebt die Darstellung in erster Linie wirkliches 
Verständnis der Begriffe und Methoden und lehnt 
jede nur rezeptmäßige Darbietung ab. Dieses hohe 
Ziel hat der Verfasser, der ja schon in seinen früheren 
Veröffentlichungen sein großes didaktisches Geschick 
bewiesen hat, in vollem Maße erreicht. Da bisweilen 


‚ cine gewisse Breite nicht gescheut wird, da gute Ab- 


bildungen, Aufgaben und Anwendungsbeispiele (mit 
durchgeführten Lösungen) den "Text beleben, ist die 
Darstellung allenthalben von ausgezeichneter Klar- 
heit und Verständlichkeit. 


N a N ir, 
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Systematisch entwickelt wird allein die Analysis; 
Algebra und Geometrie erscheinen nur als Anwen- 
dungsbeispiele oder Hilfsmittel. Auch die wichtigsten 
Verfahren der praktischen Mathematik finden Be- 
handlung. Besonders zu begrüßen ist eine moderne 
Darstellung der Grundbegriffe und Grundoperationen 
der Wahrscheinlichkeitsrechnung bis zu den Funda- 
mentalsätzen. 


Aus dem Inhalt: I. Zahlen und Zahlenfolgen. II. 
Der Funktionsbegriff. III. Das Integral und die Ab- 
leitung. IV. Die elementaren transzendenten Funk- 
tionen. V. Ergänzungen zur Differential- und Inte- 
gralrechnung. VI. Polynome, algebraische Glei- 
chungen und rationale Funktionen. Anhang: Lö- 
sungen der Aufgaben. | 


Der zweite Band soll die Integration und Differen- { 

tiation der Funktionen von mehreren Veränderlichen, 
den Abschluß der Wahrscheinlichkeitsrechnung, die 
Fehlertheorie, Ausgleichungsrechnung, lineare Glei- 
chungen und Determinanten bringen. Der dritte 
Band wird gewöhnliche und partielle Differential- 
gleichungen, NReihenentwicklungen, Integralglei- 
chungen und Variationsrechnung, der vierte Band 
schließlich die Grundzüge der Funktionentheorie 
behandeln. ! 


Auf Grund seiner großen Vorzüge darf man dem 
Werk, das bis Ende 1951 vollständig vorliegen soll, 
einen großen Benutzerkreis unter den Studierenden 
nicht nur der Technischen Hochschulen, sondern auch 
der Universitäten wünschen. Druck und Ausstattung 
sind mustergültig. 4 
Aachen. Günther Schulz, 


ee er 


Dr.-Ing. Franz Ollendorf (Prof. am Hebrew Tech- 
nical College Haifa), DieWeltderVektoren. 
Einführung in Theorie und Anwendung der Vektoren. 
Tensoren und Operatoren. VIII -+ 470 S. m. 68 Abb. 

Wien 1950. Springer-Verlag. Preis 37,50 DM, geb. | 
40,— DM. 

Der in diesem Buch gebotene Stoff geht weit: über 
den in den üblichen Büchern über Vektorrechnung 
behandelten hinaus. Am Anfang des Buches gibt der 
Verfasser selbst ein Programm dessen, was ihm bei 
der Abfassung des Buches vorgeschwebt hat. Einige 
Sätze daraus seien hier zitiert: ‚In ihrer bisher ent- 
wickelten Form ist die Vektorrechnung auf den An- 
schauungen der klassischen Physik basiert.‘ „Die 
Anwendungen der Vektorrechnung bilden... einen 
untrennbaren Teil ihrer selbst“. Die Willkür in der 
Wahl des Bezugssystems gilt es „aus den ‚„‚Gesetzen““ 
der Geometrie wie der Erfahrung zu eliminieren oder 
mit anderen Worten: „Invarianten aufzusuchen, 
welche ihren Wert bei beliebigen Transformationen 
des Bezugssystems wahren.‘ ‚‚Die Forderung der 
allgemeinen Invarianz führt von der klassischen Vek- 
torrechnung des dreidimensionalen Raumes über die 
wesentlich noch Euklidisch gebundene, be- 
schränkte Relativitätstheorie zur allgemeinen Rela- 
tivitätstheorie, die somit zur Disziplin der Vektor- 
rechnung zugehört.‘‘ „Statistische Aussagen ent- 
springen ihrem Wesen nach einer grenzenlos ge- 
steigerten Beobachtungszahl, und die aus ihnen u 
erschließenden Gesetzmäßigkeiten sind der vekto- 

‚riellen Darstellung nur in einem abstrakten Raume 
zugänglich, dessen Dimensionszahl keiner Beschrän- 
kung unterliegt: Mit der Konzeption dieses Hil- 
bertschen Raumes vollendet sich die Welt der 
Vektoren“, } ae ee 

Diesem Programm entsprechend ist das Buch auf- 
gebaut. Es gliedert den Stoff in die folgenden ach 
Kapitel: 1. Skalare und Vektoren, 2. Vektorfelder, 
3. Vektorrechnung in affinen Koordinaten, 4. Algebra 
der Tensoren, 5. Tensoranalysis im affinen Ra 
6. Der Minkowskische Raum, 7. Der Ri isch 
Raum, 8. Der Hilbertsche Raum. Auf 


ar 


} 
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einzugehen würde zu weit führen. Es kann nur jedem, 
der sich für den hier behandelten Stoff interessiert 
empfohlen werden, das außerordentlich inhaltsreiche 
Buch selbst in die Hand zu nehmen, 


Dresden, Willer.s. 


Dr, phil. A. Duschek (o. Prof. a.d. Techn. Hoch- 
schule Wien) und Dr, techn, A. Hochrainer, Grund - 
züge der Tensorreehnungin analy- 
tischer Darstellung. I. Teil: Tensor- 
analysis. VII-- 338 S. m. 64 Abb. Wien 1950. 
Springer-Verlag. Preis 24,80 DM. 

Der vorliegende zweite Band des dreibändigen Wer- 
kes behandelt alles das, was mit der Differentiatioh 
und Integration variabler Vektoren und Tensoren zu- 
sammenhängt. Da die Anwendungen in Physik und 
Technik dem dritten Bande vorbehalten bleiben sollen, 
kommen für den vorliegenden vor allem zwei Gebiete 
in Frage: die Differentialgeometrie und die Theorie 
der Felder. Das Buch beginnt daher mit der Theorie 
der Raumkurven und Flächen im Euklidischen 
Raume, die im üblichen Umfang behandelt wird, 
fügt aber am Ende noch einige Abschnitte über die 
Geometrie im Riemannschen Raume hinzu, in denen 
es besonders auf Parallelverschiebung und auf den 
.Krümmungstensor eingeht. Von der Theorie der 
Felder wird wie üblich zunächst das quellen- und 
wirbelfreie, dann das wirbelfreie und schließlich das 
quellenfreie Feld behandelt. Darüber hinaus sind 
zwei Abschnitte dem ebenen Feld gewidmet. In 
Anhängen werden die Lösungen der Aufgaben des 
ersten und zweiten Bandes gegeben. Das Buch gibt 
eine gute Einführung in die Tensoranalysis, die über 
das, was gewöhnlich in den zahlreichen kürzeren 
Büchern über Vektorrechnung geboten wird, hinaus- 

eht. Die Darstellung ist didaktisch geschickt, so 

aß das Buch durchaus empfohlen werden kann. 


Dresden. :Willers. 


\ 


Non-linear Problemsin Mechanics 
of Continua. Proceedings of Symposia in App- 
lied Mathematics. Vol.I. VII-+ 2198. New York 
1949. Published by the American Mathematical 
Society. Preis $ 5,25. 

Im Jahr 1946 faßte eine Gruppe von Mitgliedern 
der Amerikanischen Mathematischen Gesellschaft den 
Entschluß, alljährlich eine besondere Tagung über 
Probleme der Angewandten Mathematik zu veran- 
stalten, in der Erwägung, daß dieser Zweig der mathe- 

‚  matischen Wissenschaften bis dahin in Amerika nicht 
die Beachtung gefunden hätte, die ihm zukäme, Die 
erste dieser geplanten Tagungen hat als ‚Symposium 
in Applied Mathematics‘‘ vom 2. bis 4. August 1945 

an der Brown University stattgefunden und war den 
nichtlinearen Problemen der Mechanik der Kontinua 

_ gewidmet. Die Vorträge sind dankenswerterweise in 
einem gut ausgestatteten Bande zusammengefaßt, 
nur wenige davon in bloßen Auszügen. 

Fast drei Viertel des Bandes beziehen sich auf 
e hydrodynamische Probleme. Wenn man von’ Stö- 
_ .  rungs- oder Stabilitätproblemen und von den ein- 

x fachsten Aufgaben über Potentialströmungen inkom- 

. pressibler Flüssigkeiten absieht, so haben ja in der 
Tat alle hydrodynamischen Fragestellungen nicht- 
‘linearen Charakter, so daß das starke Hervortreten 
der Hydrodynamik in diesem Bande kein Wunder 

ist. Die hydrodynamischen Vorträge wiederum be- 
handeln zum überwiegenden Teil Themen der jetzt 
aktuellen Gasdynamik, die hier nicht einzeln auf- 
geführt werden können, J. J. Stoker weist in seinem 

Vortrag auf die schon oft studierte Analogie der Gas- 

E = .strd n mit den Wellen in seichtem Wasser hin. 
 _W.R. Sears behandelt Grenzschichten in einerinkom- 

ER er Flüssigkeit an schräg angeströmten Zy- 
N 


e- Beitrag zur Betz-Festschrift getan hat. Die Probleme, 
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ndern, ähnlich wie es L. Prandt] 1947 in seinem 


welche die Theorie der Freistrahlgrenzen von Poten- 
tialströmungen inkompressibler Flüssigkeiten nach 
den Arbeiten von Helmholtz, Kirchhoff, Levi-Civita 
und Villat noch bietet, werden in ihrer neueren Ent- 


„wicklung verfolgt in zwei Vorträgen von Alexander 


Weinstein und von Garrett Birkhoff. 

Im letzten Viertel des Bandes kommt die Theorie 
der Elastizität und der Plastizität zur Sprache. Die 
Berücksichtigung endlicher Formänderungen bei 
elastischen Platten führt auf grenzschichtartige 
Phänomene an den Plattenrändern, die in den Vor- 
trägen von K.O. Friedrichs und von Eric Reißner 
neben anderen Fragen der Plattentheorie studiert 
werden. Eine Methode zur Lösung der Kärmänschen 
Differentialgleichungen für die stark durchgebogene 
dünne rechteckige Platte wird in einem Vortrag von 
Samuel Levy mit Ergebnissen mitgeteilt. C. F. Car- 
rier stellt allgemeine Erwägungen an über die Stabili- 
tät elastischer Bauelemente, die kleinen zeitlich ver- 
änderlichen Oberflächenkräften unterworfen werden; 
mit ähnlichen Problemen hat sich in Deutschland 
E. Mettler beschäftigt. F.D. Murnagham erweitert 
die Elastizitätstheorie auf sehr starke Beanspru- 


' chungen, wobei insbesondere die Kompression eines 


Mediums durch extremen hydrostatischen Druck und 
die Dehnung eines elastischen Zylinders durch starken 
Zug betrachtet werden. Plastisches Verhalten von 
Stoffen mit Verfestigung erörtert der Vortrag von 
D.C. Drucker. W. Prager betrachtet Unstetigkeiten 
der Spannungen bei statisch bestimmten plastischen 
Problemen, Unstetigkeiten, die etwa den Verdich- 
tungsstößen der Gasdynamik entsprechen. 4 

Aus dieser notgedrungen unvollständigen Über- 
sicht dürfte bereits zur Genüge erhellen, daß dies 
erste Symposium in Applied Mathematics in seinem 
dokumentarischen Niederschlag die aufmerksame 
Beachtung der Fachgenossen verdient. 


Göttingen. W. Tollmien, 


Dr.-Ing. W. Meyer zur Capellen. Integraltafeln: 
Sammlung unbestimmter Integrale elementarer Funk- 
tionen. VIIT-+ 292 8. Berlin-Göttingen-Heidelberg. 
Springer 1950. Preis geb. 36,— DM. 

In der Tafel sind etwa 3000 unbestimmte Integrale 
elementarer Funktionen zusammengestellt, auch 
solcher, die auf vertafelt vorliegende sog. höhere 
Funktionen führen oder deren Integrale durch Reihen 
leicht zu bestimmen sind. Dabei ergeben sich ge- 
schlossene Ausdrücke oder auch Rekursionsformeln, 
aus denen in einigen Fällen auch das Schlußergebnis 
berechnet ist. Einteilung und Aufbau der Tafeln sind 
sehr geschickt. Weitgehende Unterteilung und ein 
ausführliches Inhaltsverzeichnis erleichtern das Auf- 
finden eines gesuchten Integrals. Eine Zusammen- 
stellung der wichtigsten Integrationsregeln und Hin- 
weise auf den Gebrauch der Tafeln leiten das Buch 
ein, eine Zusammenstellung wichtiger Foımeln, ins- 

‘besondere für die vorkommenden höheren Funktionen 
und von vielgebrauchten Konstanten, sowie ein 
Schriftenverzeichnis schließen es ab. Satz, Druck und 
Ausstattung des Buches sind musterhaft. 

Dresden. Willers. 

Adalbert Prey, Einführungindie sphä- 
rischeAstronomie. VII-+ 316 S. mit 123 Abb. 
Wien 1949. Springer-Verlag. Preis 22,— DM., geb. 
24,— DM. 

Das vorliegende Lehrbuch ist nicht nur für Stu- 
dierende der Astronomie gedacht, sondern darüber 
hinaus für alle diejenigen, die auf Nachbargebieten 
arbeiten. Es bringt eine knappe, wohlerwogene Aus- 
wahl aus dem Gesamtgebiet der Astronomie unter 
Ausschluß der Astrophysik. Der erste Abschnitt be- 
handelt u. a. die sphärische Trigonometrie, die Koor- 
dinatensysteme, die scheinbaren Bewegungen der 
Himmelskörper, die Zeitrechnung, Sonnenuhren, 


% 
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Finsternisse, Präzession, Nutation, Aberration, Paral- 
laxe, Refraktion. Im zweiten Abschnitt wird die 
Theorie der astronomischen Instrumente und Hilfs- 
apparate gebracht (u.a. Spiegelsextant, Universale, 


Meridiankreis- und Passageinstrument sowie das, 


Äquatoreal). Der dritte Abschnitt schließlich ist der 
geographischen Ortsbestimmung gewidmet. Ein An- 
hang bringt das Wichtigste aus der Ausgleichungs- 
rechnung. 

Die Darstellung verrät allenthalben die große Lehr- 
erfahrung des Verfassers und zeichnet sich durch 
Klarheit und gute Verständlichkeit aus. Wie der 
Verfasser selbst im Vorwort bemerkt, haftet einem 
Buch über sphärische Astronomie leicht etwas Alt- 
modisches an, da im Verlauf der letzten 100 bis 150 
Jahre nichts grundsätzlich Neues hinzugekommen sei. 
Es entsteht in der Tat bisweilen der Eindruck, als ob 
die rechnende Astronomie, eines der ältesten An- 
wendungsgebiete der Mathematik, an der stürmischen 
und erfolgreichen Entwicklung der praktischen Mathe- 
matik allzd wenig Anteil genommen habe. Bedauer- 
lich ist, daß die meisten Abbildungen von Kugeln mit 
ihren Großkreisen und Dreiecken im Sinne der Dar- 
stellenden Geometrie falsch gezeichnet sind, obwohl 
richtig konstruierte Bilder erfahrungsgemäß von 
größerer Anschaulichkeit sind. 


Aachen, Günther Schulz. 


Heinrich Rödel, Mechanik für Ingenieure, 
Bd. 4: Dynamik einschließl. Schwingungslehre, 
(Westermanns Fachbücher der Ingenieur-Fach- 
kunde.) 198 8. m. 242 Abb. Braunschweig-Berlin- 
Hamburg 1949. Georg Westermann-Verlag. Preis 
kart. 9,20 DM, geb. 10,40 DM. 

Ähnlich wie bei der Blasiusschen Mechanik be- 
vorzugt hier der Verf. die anschauliche Methode und 
stellt dem Leser ein reiches Aufgabenmaterial zur 
Verfügung, gegliedert in Beispiele, die vollständig 
durchgerechnet sind, und Aufgaben, die der Leser 
sebst lösen soll und deren Lösungen am Schluß des 
Buches angegeben sind. Die Theorie beschränkt sich 
auf die allereinfachsten Ableitungen, da das Buch 
hauptsächlich für den nichtakademischen Ingenieur 
bestimmt ist. Hinsichtlich des Aufbaues wäre zu be- 
merken, daß die Anwendung des Impulssatzes auf 
Stoßvorgänge erst am Schluß des Buches gebracht 
wird, während vorher schon ausgiebig vom Energie- 
satz Gebrauch gemacht wurde. U. a. wird im Beispiel 
48 die Bewegung eines schweren Klotzes behandelt, 
der auf eine abgefederte Unterlage herabfällt. Verf. 
vernachlässigt hierbei die Masse der Unterlage, ohne 
es besonders zu erwähnen. Bei ihrer Berücksichti- 
gung würde eih Stoßverlust auftreten, der die An- 
wendung des Energiesatzes in der einfachen Form 
nicht mehr gestattet. 


Dresden. 15 Neuber.. 


Dr. E. Madelung (o. Prof. a. d. Univ. Frankfurt 
a.M.), Die mathematischen Hilfsmittel des 
Physikers.. (Mathematische Wissenschaften in 
Einzeldarstellungen mit besonderer Berücksichtigung 
der Anwendungsgebiete Bd. IV.) Vierte, vermehrte 
und verbesserte Aufl. XX + 5318. mit 29 Abb. 
Berlin-Göttingen-Heidelberg 1950. Springer-Verlag. 
Preis brosch. 47,— DM, geb. 49,70 DM. 

Die dritte Auflage dieses Buches war 1942 ver- 
griffen. In sechsjähriger Arbeit hat es der Verfasser 
dann einer gründlichen Überarbeitung unterzogen und 
dabei manches hinzugefügt, was heute für den Phy- 
siker von Wichtigkeit ist. So hat er z.B. dem 
mathematischen Teil, der etwa drei Viertel des Um- 
fanges ausmacht, ein Kapitel „Zahlen, Funktionen, 
Operatoren‘ vorangestellt. Im einzelnen alle Ande- 
rungen aufzuzählen ist unmöglich; denn es findet sich 
kein Abschnitt, der gänzlich unverändert geblieben 


’ 
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ist. Als Beispiel sei nur erwähnt, daß in der Gruppen- 
theorie ein Abschnitt über Symmetriegruppen (Kri- 
stalle) hinzugefügt, die Abschnitte über Vektor- 
analysis, spezielle Koordinaten usw. wesentlich er- 
weitert, der über Störungstheorie neu geschrieben 
wurde usw. Im physikalischen Teil wurden die Ab- 
schnitte ‚‚Relativitätstheorie, Quantentheorie und 
Thermodynamik“ neu abgefaßt. Der Umfang des 
Buches hat sich so um etwa 150 Seiten vermehrt. 

Über die Art der Darbietung sagt der Verfasser 
im ‘Vorwort: ‚Ich weiß, daß mein Buch in vielen 
Punkten nicht die Zustimmung aller Mathematiker 
findet. Sie tadeln die Leichtfertigkeit, mit der ich 
mit ihrem geheiligten Besitz umgehe. Aber vieles, 
was jenen besonders wichtig und interessant erscheint, 
ist für uns ohne Bedeutung. Wir sind oft gezwungen, 
die Gefahren des Leichtsinns auf uns zu nehmen und 
können die Maschine, mit der wir arbeiten wollen, 
nicht mit gar zuviel Sicherungen und Warnungs-, 
schildern umgeben, ohne uns empfindlich zu hindern. 
Gegenüber dem Mathematiker ist der Physiker immer 
ein Praktiker.‘ 

Das Buch hat sich in der Hand der Physiker seit 
langem bewährt und wird auch in der neuen Form 
allen Physikern von außerordentlichem Nutzen sein. 


Dresden. Willers. 


Dr.-Ing. Winfried Oppelt, Stetige Regelvor- 
gänge. (Bücher der Technik. Herausgeber: Dr.-Ing. 
Alfred Kuhlenkamp.) 144 S. mit 42 Abb. Han- 
nover-Wolfenbüttel 1949. Wissenschaftl. Verlags- 
anstalt K. G. Hannover und Wolfenbütteler Verlags- 
anstalt G.m.b.H. Preis brosch. 7,— DM., geb. 
7,80 DM. 

Verfasser gibt anschließend an sein Buch ‚‚Grund- 
gesetze der Regelung‘‘ (Besprechung in Z. angew. 
Math. Mech. Bd. 29 [1949], S. 62) eine systematische 
Übersicht über stetige Regelvorgänge. Obwohl als 
Einführung gedacht, ist an zahlreichen Stellen eine 
handbuchartige kurze Form gewählt, die allerdings 
durch zahlreiche (etwa 100) geschickt gewählte 
Skizzen und graphische Darstellungen wirkungsvoll 
ergänzt wird. Auf die zahlenmäßige Durchrechnung 
von Beispielen glaubt der Verfasser verzichten zu 
können. (Ref. ist allerdings der Ansicht, daß dadurch 
dem. Anfänger die Einarbeitung auf dem Gebiet der 
Regelung erschwert wird, erst Zahlen geben den 
richtigen Einblickin die Wirksamkeit einzelner Regel- 
glieder.) 

Der Aufbau von Regelsystemen erfolgt in dm ı\ 
Buch synthetisch aus einfachsten Einzelbausteinen, 
die sich in der Praxis jeweils leicht verwirklichen 
lassen. Die Wirkung von Zusammenschaltungen wird 
vorzugsweise in schematischer Form graphisch ver- 
folgt. Überall wird die Teilung der Anlagen in Regel- i 
strecke und Regler weitgehend ausgenutzt. Der NG 
systematische Aufbau reicht mit einer kurzen Be- | 
merkung bis zur Behandlung von Mehrfachregelungen. IcuE 

Dem Buch ist eine Fehlerberichtigung zu W. Ir 
Oppelt, Grundgesetze der Regelung beigelegt. = 

Dresden, N. J. Lehmann. 


Prof, Dr.-Ing. A. Leonhard, Die selbsttätige 
Regelung. (Theoretische Grundlagen mit pak- 
tischen Beispielen). IX + 284 S. mit 254 Abb. Ber- : 
lin-Göttingen-Heidelberg 1949. Springer-Verlag. Preis 
brosch. 24,— DM, geb. 27,— DM. ah Er 

Das vorliegende Buch stellt eine wertvolle Be- Br a 
reicherung der Literatur über die Regelungstechnik 
dar. Obwohl zunächst als Einführung gedacht, führt 
es bis zu den modernsten Arbeiten üb r beispielsweise 


die optimale Wahl von Regelungskonstanten und ei t 


P 


damit auch dem erfahrenen Fachmann neue 
regungen. Die Darstellung ist diesem Zw: 
paßt. Mit einer etwas Be € 


wird das gesamte Gebiet der se bsttätig 
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Regelung in leicht verständlicher Form aufgebaut. 
Zahlreiche instruktive Beispiele, Skizzen und gra- 
phische Darstellungen sowie viele auch zahlenmäßig 
durchgerechnete Aufgaben machen den Leser schnell 
mit den Problemen und Methoden der Regelungs- 
technik vertraut. Sie vermitteln das auf diesem 
h . Gebiet besonders nötige „Fingerspitzengefühl“. Die 
g erforderlichen mathematischen Grundlagen — ins- 
» besondere für die Stabilitätsuntersuchung — sind 
soweit dargestellt, daß bei der Verwendung kaum 
Schwierigkeiten auftreten können; wie überall, so 
steht auch hier besonders die Anwendung im Vorder- 
grund. 

B Einen Einblick in den reichhaltigen Inhalt des 
- Werkes vermittelt die Aufzählung einiger charakte- 
ristischer Kapitel und Abschnittsüberschriften: Allge- 
meines über Regelungen, Verhalten von Meßwerken 
und Einzelgliedern des Verstellsystems, Ermittlung 
R- des Regelvorganges (auf dem Wege über die Diffe- 
rentialgleichung, Heayviside-Kalkül, Laplace Trans- 
formation, graphische Methoden), die Stabilität der 
E- Regelung (nach der Methode der kleinen und der 
. selbsterregten Schwingungen), Kontrolle der Dämp- 
; fung, Verbesserung der Stabilität, Festlegung frei 
| wählbarer Regelkonstanten (nach verschiedenen Me- 

thoden und Gesichtspunkten). 


Dresden, 


N.J. Lehmann. 


Dr. Ing. E. Schmidt (o. Prof. a. d. Techn. Hoch- 
schule Braunschweig), Einführung in die 
technische Thermodynamik undiin 
die Grundlagen der chemischen 
Thermodynamik. Vierte erweiterte und über- 
arbeitete Auflage. XVI -+ 520 S. m. 244 Abb., 69 Ta- 
bellen u. 3 Tafeln. Berlin-Göttingen-Heidelberg 1950. 
Springer-Verlag. Preis geb. 30 DM. 

Die Thermodynamik vonE.Schmidt zählt, seit 
sie vor rund 20 Jahren die Fortsetzung des Schüle- 
schen Werkes übernahm, zu den anerkannt besten 

- Lehrbüchern des Fachgebietes und iistin der Standard-' 
Literatur ein nicht mehr hinwegzudenkender Begriff 
geworden. Eine grundsätzliche Würdigung des 
Buches darf sich daher auf die erneute Kennzeichnung 
seines Wertes beschränken, der vor allem in der aus- 


heit und Folgerichtigkeit der thermodynamischen 
Entwicklungen und in der durch instruktive Zahlen- 
beispiele wirksam unterstützten Anwendung der ge- 
wonnenen theoretischen Erkenntnisse auf die wichtig- 
sten technischen Maschinen und Geräte zu erblicken 
ist. - 

Von den Erweiterungen, die die neue Auflage be- 
reichern, scheint zwar das Kapitel Strömungsmaschi- 
nen des Guten etwas zu viel zu tun, da es sich mit 


dichterbaues beschäftigt, die für den daran beteiligten 
Konstrukteur nur Fragmente bedeuten, für den Ther- 
'  modynamiker aber zu abseits liegen. Dagegen ist 
außerordentlich dankenswert die ausführliche Be- 
handlung der Thermodynamik des Raketen- und des 
 Luftstrahlantriebes, die sich in geschlossener Dar- 
. stellung sowohl mit den verbrennungs- und strömungs- 
technischen Vorgängen bei der Energieumsetzung als 
- auch mit den mechanischen Bedingungen des Raketen- 
fluges beschäftigt. 
Über den Wert der letzten beiden, ebenfalls neuen 
Abschnitte aus der chemischen Technik für das vor- 
liegende Werk kann man geteilter Meinung sein, Sie 
_ enthalten nur kleine Teilprobleme aus dem großen 
Fachgebiet, nämlich die Wärmetönung, das Massen- 
wirkungsgesetz, das Nernstsche Wärmetheorem und 
Anwendung auf Verbrennung, Vergasung und 
ziation. Für den Wärmeingenieur sind diese 
obleme nur von geringem Interesse. Den Che- 
enieur können sie allein nicht befriedigen. Sie 
[olgedessen die Geschlossenheit des Buches 


führlichen Erörterung der Grundbegriffe, in der Klar- 


reinen Spezialproblemen des Turbinen- und Ver- 


und sollten, wie der Verfasser selbst im Vorwort an- 
deutet, besser einem besonderen, der chemischen 
Thermodynamik gewidmeten Werke vorbehalten 
bleiben, 

Dafür wäre in dem Abschnitt Wärmeübertragung 
eine ausführlichere Behandlung der Vorgänge beim 
Kondensieren und Verdampfen, sowie ein Eingehen 
auf die Prandtlsche Grenzschichttheorie, der die 
Mehrzahl der im Kapitel 109 angeführten Gleichungen 
für den Wärmeübergang entstammen, und auf die nur 
angedeutete Flammenstrahlung, zu empfehlen. 


Dresden, Br:Rraltın, 


Elecetromagnetie Theory (Proceedings of the second 
Symposium in Applied Mathematics ot the American 
Mathematical Society.) 91 8. mit.8 Abb. New York 
an American Mathematical Society. Preis geb. 
38. 

Von einer 1948 von den amerikanischen Institutes 
of electrical Engineers, of Physics und of Radio 
Engineers veranstalteten Tagung sind 17 Arbeiten, 
meıst auszugsweise wiedergegeben, Hervorgehoben 
seien: 

1. Eine neue Quanten-Elektrodynamik (H. Fesh- 
bach). Für die in letzter Zeit für das „Vakuum“ 
entworfene, interessante Schwankungstheorie wird 
eine kovariante, maßstabunabhängige Form für die 
Wechselwirkung zwischen Strahlung und Masse auf- 
gestellt. 2. W. H. Watson gibt in ‚„Diskontinuitäten 
ın der elektromagnetischen ‘[heorie‘‘ die Zusätze zu 
den Maxwellschen Gleichungen an, die diskonti- 
nuierliche Elektronenbewegungen zu beschreiben ge- 
statten, z. B. die Bahnsprünge, 3. L. Infeld erweitert 
die Technik der von Dirac stammenden ‚‚factori- 
zation method“, die die Berechnung von Eigenwert- 
problemen vereinfacht. 4. J. L. Synge behandelt 
ın „Elektromagnetismus ohne Metrik‘“ 2 schiefsym- 
metrische Tensoren, um zu den Grundlagen für eine 
Vereinigung der Gravitations- und elektrodynami- 
schen Erscheinungen beizutragen, 

Von den weiteren, oft nur 1 Seite umfassenden 
Referaten, seien erwähnt: Eine Methode der analy- 
tischen Fortsetzung bei Eigenwert- und Streupro- 
blemen der Quantentheorie (EB. Feenberg); Systeme 
von Wiener-Hopf Integralgleichungen und ihre An- 
wendung auf Randwertprobleme des Elektromagnetis- 


- mus (A. Heins); Nichtlineare elektrische Netzwerke 


(J. Duffin); Reflexionen an Krümmungen und 
Ecken in Hohlrohrleitern (Rice); Wellenausbreitung 
in Hornstrahlern (F, Stevenson); Kompensation 


sphärischer Elektronenlinsenfehler durch Spiegel 
(Ramberg). 
Dresden. H. Schönfeld, 


Iris Runge, Carl Runge und sein wissen- 
schaftliches Werk. (Sonderdruck aus den Ab- 
handlungen der Akademie der Wissenschaften in 
Göttingen. Mathematisch-physikalische Klasse.) 
2148. mit 2 Abb. Göttingen 1949. Vandenhoek & 
Ruprecht. Preis geb. 13,— DM. 

Mit Carl Runges Namen ist die Entwicklung der 
angewandten Mathematik in Deutschland auf das 
Engste verbunden. Schon bevor er 1904 den Lehr- 
stuhl für angewandte Mathematik an der Universität 
Göttingen, den ersten dieser Artin Deutschland, über- 
nahm, hatte er sich bemüht die mathematischen 
Methoden den Bedürfnissen der Techniker anzu- 
passen und bequem anzuwendende Verfahren zu 
entwickeln, die bis zu zahlenmäßigen Ergebnissen 
führten. Dazu befähigten ihn seine hervorragenden 
mathematischen Kenntnisse (er war Schüler von 
Weierstraß, Kronecker usw.), seine große 


Energie, mit der er die bei Ingenieuren, Vermessungs- 


technikern usw. gebräuchlichen speziellen Methoden 
sich aneignete und so durcharbeitete und verall- 


j 
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gemeinerte, daß sie allgemein verwendbar wurden, 
und sein Blick für die Bedürfnisse der Praxis. Die 
von ihm entwickelten Verfahren bilden heute einen 
wesentlichen Teil aller Vorlesungen über graphische 
und numerische Methoden und über praktische 
Analysis, wie sie heute an unseren Universitäten und 
Hochschulen gehalten werden. 

Das vorliegende Buch, von Runges ältester 
Tochter geschrieben, zeichnet liebevoll seinen 
Lebenslauf, schildert Runge als harmonisch aus- 
geglichene Persönlichkeit und geht auf seine wissen- 
schaftlichen Leistungen vor allem auf dem Gebiet 
der Spektroskopie ein. Daneben entwirft sie ein an- 
schauliches Bild der damals führenden Männer der 
mathematischen und naturwissenschaftlichen For- 
schung und des Lebens in diesen. Kreisen, insbeson- 
dere des Lebens in Göttingen, das in jener Zeit die 
„Hochburg der Mathematik‘ war. Jedem, der Runge 
gekannt hat, wird die Lektüre des Bandes große 
Freude bereiten. 


Dresden. Willers. . 


Ernst Lindelöf, EinführungindieHöhere 
Analysis. Zum Selbststudium und für Studieren- 
de der ersten Semester. Deutsch von Egon Ulb- 
rich (Prof. a. d. Universität Gießen). IX+ 
526 S. m. 84 Abb. 2. Aufl., Leipzig 1950. B. G. Teub- 
ner Verlagsgesellschaft. Preis geb. 14,80 DM. 

Das vorliegende Buch hat den Zweck, „einen 
stetigen Übergang von den Schulstudien zu den 
Universitätsstudien zu vermitteln“, und ist außerdem 
besonders auf die Bedürfnisse eines Studenten des 
höheren Lehramts zugeschnitten. Das Ergebnis dieser 
Zielsetzung ist ein Werk, das sich in seinem Aufbau 
weitgehend von den übrigen im Gebrauch befind- 
lichen Büchern ähnlichen Inhalts unterscheidet. Eine 


kurze Inhaltsübersicht vermittelt am besten einen 
Eindruck von der Anordnung und dem Charakter 
des Buches: 

Ein ausführliches Einleitungskapitel bekandelt die 
elementaren Funktionen, untermauert so den Schul- 
stoff und schafft eine solide Grundlage für das weitere 
Studium. Nach dem Rechnen mit Näherungswerten, 
das mit einer Strenge behandelt wird, die man diesem 
Gebiet sonst nicht zu widmen pflegt, werden Ketten- 
brüche und Grenzwerte (einschließlich Reihen mit 
konstanten Gliedern) gebracht. Es folgt ein Kapitel 
über Ableitungen von Funktionen, eine Differential- 
rechnung bis zum Mittelwertsatz und der Bestimmung 
von Extremwerten. Ein Zwischenkapitel, das der 
exakten Einführung der Begriffe Länge, Flächeninhalt 
und Volumen gewidmet ist, leitet über zur Integral- 
rechnung, dem Kernstück des Buches. Bis hierher 
werden gewisse Sätze aua der Irrationalzahltheorie 
sowie über Grenzwert und Stetigkeit als bekannt 
oder evident vorausgesetzt. Diese Sätze werden erst 
jetzt, bei hinreichender Reife .des Lesers in voller 
Exaktheit bewiesen und dürften hier auf besseres 
Verständnis stoßen, als dies bei dem üblichen syn- 
thetischen Aufbau der Fall ist. Das Buch wird be- 
schlossen mit einem Kapitel über komplexe Zahlen 
und drei Anhängen über lineare Gleichungssysteme 
bis zu 3 Unbekannten und über Determinanten zur 
Inhaltsberechnung von Polygonen und Polyedern. 
Zahlreiche Aufgaben (allerdings ohne Lösungen) bieten 
genügend Übungsstoff. Die Darstellung ist sehr breit, 
lebendig und außerordentlich leicht lesbar und macht 
das Buch zur einem ausgezeichneten Hilfsmittel für 
Lehrer und Lernende, nicht zuletzt auch ein Ver- 
dienst: des Übersetzers. 


Dresden. ha G. Opitz. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER Rn 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten) 


B.L, van der Waerden, (Prof. a. d. Univ. Amster- 
dam), Moderne Algebra. Unter Benutzung 
von Vorlesungen von E. Artin und E. Noether. 
Erster Teil. (Die Grundlehren der mathematischen 
Wissenschaften in Einzeldarstellung mit besonderer 
Berücksichtigung der Anwendungsgebiete. Band 
XXXIII). 3. Aufl. VIII -+ 292 Seiten. Berlin-Göt- 
tingen-Heidelberg 1950. Springer-Verlag. Preis.brosch . 
24,— DM, geb. 27,— DM. 


Dr. R. Grammel, (Prof. a. d. Techn. Hochschule 
Stuttgart), Der Kreisel, seine Theorie 
und Anwendung. Zweite neubearbeitete Aufl. 
ErsterBand: Die Theorie des Kreisels.. XI -+ 2818. 
mit 137 Abb. Berlin-Göttingen-Heidelberg 1950. Sprin- 
ger-Verlag. Preis brosch. 30,— DM, geb. 33,— DM, 


Dr. Wilh, Maak, (apl. Prof. a. d. Univ. Hamburg), 
Differential-undIntegralrechnung 
(Bücher der Mathematik und Naturwissenschaften. 
Herausgeg. von Dr. Henry Poltz.) 235 S. mit 17 Abb. 
Wolfenbüttel-Hannover 1949. Wolfenbütteler Ver- 
lagsanstalt. 


Dr. H, Hasse, (o. Prof. a. d. Humboldt-Univ. Berlin), 
VorlesungenüberZahlentheorie. (Die 
Grundlehren der Mathematischen Wissenschaften in 
Einzeldarstellungen mit besonderer Berücksichtigung 
der Anwendungsgebiete.) Band LIX. XII-+ 4748. 
mit 28 Abb. Berlin-Göttingen-Heidelberg 1950. 
Springer- Verlag. Preis 42,— DM, geb. 45,— DM. 


Dr. W. Maak, (Prof. a. d. Univ. Hamburg), Fast- 
periodische Funktionen. (Die Grund- 
lehren der mathem. Wissenschaften in Einzeldarstel- 
lung mit besonderer Berücksichtigung der Anwen-. 
dungsgebiete. Band LXIL) VIIT+ 240 8. Berlin- 
Göttingen-Heidelberg 1950. Springer-Verlag. Preis 
21,60 DM. geb. 24,60 DM. 


H.Krup, Bestimmungder allgemeinen 
Lösung der Schrödinger-Gleichung 
für Coulomb-Potential. (Berichte über 
die Verhandlungen der Sächs. Akademie der Wissen- 
schaften zu Leipzig. Mathematisch-naturwissensch. 
Klasse. Band 97, Heft 8.) 28 S. Berlin 1950. Aka- 
demie-Verlag. 


H. Schubert, Übereinelinearelntegro- 
differentialgleichung mit Zusatz- 
kern. (Berichte über die Verhandlungen der Sächs, 
Akademie der Wissenschaften zu Leipzig. Mathe- 
matisch-naturwissenschaftl. Klasse. Band 97, Heft 7). 
52 S. Berlin 1950. Akaderierfeias, Preis Rn DM. 


Dr. Fr. Hand; (0.\Prot a ANCDr oder Bin... 


führungin "die theoretische Physik, 


Band 5. Atom- und Quantentheorie. (Meyers kleine 


Handbücher, Band 56/57.) 314 S. mit 48 Abb. Leip- 


zig 1950. Bibliographisches Institut. Preis 5,80.DM, 2 
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ZUSCHRIFITEN AN DEN HERAUSGEBER 


Zu K.Ludwig: „Die Biegung der Rechteckplatte 
ohne die Bernoullischen oder andere Annah- 
men“. Z.angew. Math. Mech.29 (1949), S.18—19. 


Die Verschiebungen £, n, & sind als Polynome von ? 
unter Beschränkung auf die nullte, erste, zweite und 
dritte Potenz angesetzt.- Unter der vereinfachenden 
Annahme, daß es sich um ein statisches Problem han- 
delt, verschwinden die Diff. Quotienten nach {, und 
die &,, &, usw. sind nur noch Funktionen von x und 
y. Es stellen sich dann alle Spannungs- und Ver- 
schiebungsgrößen als Polynome von z bis zum dritten 
Grade dar. 

Damit enthalten die Gleiehgewichts-Gleichungen (Be- 
wegungs-Gleichungen)sebenfalls solche Polynome, z. B. 
lautet, die Gleichgewichtsgleichung der x-Richtung: 


ih 
alarm sr $) 
Kt [4 Be en & a : =) 
el re) 
su Hei) 


Es müssen hierin die Koeffizienten für jede der vier 
Potenzen von z verschwinden. Der Verfasser setzt 
aber nur die Kosffizienten von 2° und 2! null, gewinnt 
damit 2 Bestimmungsgleichungen für &, und £,; so 
erhält er, weil er nicht auch die Koeffizienten von 2? 
und 2° verschwinden läßt, freie Verfügung über die 
Funktionen &, und &£, und verwendet diese dann zur 
Befriedigung von Randbedingungen. Entsprechend 
verfährt er in den Gleichgewichts-Gleichungen der y- 
und z-Richtung, womit er die Verfügung über die 
Funktionen n, und n, bzw. £, und Z, behält. 

Damit sind jedoch die Gleichgewichtsbedingungen 
nicht erfüllt. 


Darmstadt. A. Mehmel, 


Erwiderung 
Die Verschiebungskomponenten £&, 7 und £ sind bis 
z° entwickelt. Die Längsdehnungen 


KT 3% und 


‚sind auch bis 2° entwickelt. Die Querdehnung 
; 0% 


NR dr 


ist dagegen nur bis 2? entwickelt. Daher sind die 
Normalspannungen auch nur.bis 2? entwickelt. Der 
Scherungswinkel | 


also auch die Schubspannungskomponente T;y sind 
bis 2? entwickelt. Die Scherungswinkel 


207 0 _mn,% 
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also auch die Schubspannungskomponenten ry,; und 


Tr; sind bis 2? entwickelt. Die Komponenten 
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Da die Summanden mit 2? und 2? doch nicht voll- 
ständig sind, wäre ihre Beseitigung sinnlos. In dieser 
1. Näherung werden die Summanden bis z beseitigt, 
in 2. Näherung die mit z?, in 3. die mit 2°,... Wenn 
die Bewegungsgleichung nicht angenähert sondern 
streng erfüllt würde, hätte man die Lösung des räum- 
lichen Problems. Dann aber brauchte man die Platten- 
theorie nicht mehr, 


Hannover. K, Ludwig. 


Entgegnung 


Nur ‚durch das Differenzieren nach 2 werden die 
Exponenten von 2 kleiner. Damit bleibt die dritte 


\ ; IMG: 
Potenz von z sowohl in e, und o; wie auch in —: er- 


0x 
halten. Sie kommt demnach auch in den Bewegungs- 
gleichungen vor, wie ja schon in meinem Schreiben 
angedeutet war. 

Die Gleichgewichtsbedingungen sind also tatsächlich 
nicht erfüllt. Daß Sie die höheren Potenzen von z in 
einem Iterationsverfahren beseitigen wollen, ist mir 
an Hand Ihres Aufsatzes nicht klar geworden. Viel- 
mehr hat mich auch die Überschrift zu der Annahme 
verleitet, daß es sich hier um die strenge Lösung des 
räumlichen Plattenproblems handelt. 

Ich glaube auch nicht, daß durch die Hinzunahme 
weiterer Potenzen von 2 und schrittweise vorgehende 
Ausschaltung der niederen Potenzen das Ziel zu er- 
reichen ist. Durch den Ansatz weiterer Potenzen von 
2 bekommt man immer genau so viel frei verfügbare 
Koeffizienten &£ (x, y), nk (8, y), &k (x, y) wie man 
Bestimmungsgleichungen hat. Will man alle Gleich- 
gewichtsgleichungen befriedigen, so bleiben keine 
Koeffizienten für die Erfüllung der Randbedingungen 
übrig. Im umgekehrten Fall ist die Befriedigung der 
Gleichgewichtsbedingungen nicht möglich, es sei denn 
daß der Nachweis für die Bedeutungslosigkeit der, 
Glieder mit hohen Potenzen von z erbracht wird. 

Ich möchte noch auf etwas anderes hinweisen. Ver- 
einfacht man das Problem auf die Weise, daß man die 
Belastung zur Hälfte oben als Druck, zur anderen 
Hälfte. unten als Zug angreifen läßt, und legt den 
Koordinatenanfang in die Mittelebene der Platte, 
so muß x 

€ und nin z ungerade, 
3 Zin 2 gerade sein. 
Für diesen Fall ist also 


=s>0 
Die hier fehlenden Glieder stellen also den Einfluß 
einer Belastung dar, die oben und unten als Druck 
wirkt und mithin für sich im Gleichgewicht ist. Für 
p= const. z. B. ist der Fallelementar. Dieser in z sym- 
metrische Belastungsfall ist wirklich von untergeordne- 
ter Bedeutung, und man kann nicht durch seine Hinzu- 
nahme die Gleichgewichtsgleichungen des obigen An- 
satzes in Ordnung bringen. 


Darmstadt. A.M e hmel. 


Zu O.Emersleben, „Die Schwingungsdauer 
eines umlaufenden Pendels als Analogon zum 
Potential eines Kreises‘. Z. angew. Math. u 


- Mech. 29 (1949) S. 279-282 und 30 (1950) $. 64° 


Der Verfasser gibt das Potential u einer auf einem 
Kreis gleichförmig verteilten Ladung E im Auf- 
punkt P mit den kleinsten und größten Abständen 
0, bzw. o, des Aufpunktes vom Kreis mit 
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an. Dabei ist K(k) das vollständige elliptische Inte- 
gral 1. Gattung in der Legendreschen Normal- 
form mit dem Modul k. Das Potential liegt zwischen 
den beiden Werten, die sich ergeben, wenn man sich 
E nacheinander in den Punkten des Kreises konzen- 
triert denkt, die vom Aufpunkt die Abstände o, und 
0, haben. Die Bestimmung dieses Zwischenwertes 
leistet nun das ‚„arithmetisch-geometrische Mittel‘ 
eleganter als K(k). Man erhält das arithmetisch- 
geometrische Mittel M., z. B. von o, und o,, indem 
man nacheinander von 0, und o, arithmetisches und 
geometrisches Mittel bildet, von den erhaltenen 
Werten wieder bis beide Werte binreichend gleich 
sind. Die Konvergenz ist sehr stark. 


Nach Gauß!) ist nun 
n/2 


Ki ——— 
Mas (1,/1 — R2) eo Ni 


Aus (1) folgt dann sofort 


E 
U ET RES 
M as (01, 03) (9) 
Es ist das eine Lösung, die bei technischen Pro- 
blemen die Auswertung mit dem Rechenschieber 
gestattet. 

Selbstverständlich bleibt die vom Verfasser ge- 
fundene Analogie zur Schwingungsdauer eines um- 

laufenden Pendels gewahrt. Hierfür gilt 


RE REN 


Legt man auf die formale Analogie weniger Wert, so 
erhält man mit o&, = 2z/T, und @, = 2n/T, 

> 5 
Ua UT) 


m 


oder 
Be NEN 
Mas(T, T,) 
Auch alle diese Lösungen sind anschaulicher und 


bestätigen erneut?) die vielseitige Anwendbarkeit 
des May. 


Stuttgart. 


Ki Ppotthott 


Schrifttum. 

[11 K.Kommerell: Das Grenzgebiet der elemen- 
taren und höheren Mathematik. F. K. Koehler 
Verlag, Leipzig 1936. 

[2] K. Potthoff: Die Berechnung der Gegen- 
induktivität koaxialer kreisförmiger Windungen. 
2. f. Elektrotechnik 2 (1949) S. 20 und-S$. 139. 


Erwiderung 


Gleichungspaar (0) meiner kleinen Mitteilung, auf 
das Herr Potthoff vorstehend unter (1) zurück- 
kommt, ist aus einer zitierten älteren Arbeit (1925) 
des Verfassers übernommen, wo es die Grundlage für 
eine Ausgestaltung des Potentials einer kreisförmigen 
Ladung gebildet hat!) und wo u.a, gezeigt wurde, 


‘) Zs. f. Hochfrequenztechnik 26 (1925) 61— 73. 


daß sich auf solche Ladungen insbesondere ein für 
Punktladungen bekanntes Spiegelungsverfahren an- 
wenden läßt. Die Formeln haben sich für die Durch- 
führung numerischer Rechnungen praktisch bewährt, 
vgl. z. B. die Bilder 5 und 7 meiner damaligen Ver- 
öffentlichung!). Für die dabei vorkommenden 
elliptischen Integrale wurde von den dafür vor- 
handenen bequem zugänglichen Tafeln Gebrauch ge- 
macht. — Wenn ich dabei statt jedes Tabellenwertes 
ein arithmetisch-geometrisches Mittel hätte berechnen 
sollen, wäre selbst angesichts der bekannt schnellen 
Konvergenz dieser Zahlenfolgen und auch bei durch- 
gehender Beschränkung auf Rechenschiebergenauig- 
keit die Arbeit eine vielfache geworden. 


Herr Potthoff betont gelegentlich — Zeit- 
schrift für Elektrotechnik 2 (1949) 21 — als Vorteil 
für die Verwendung des arithmetisch-geometrischen 
Mittels, daß sich die Rechnungohne Gebrauch 
von Hilfstafeln mit dem Rechenschieber 
durchführen lasse, Da jedoch die Beschaffung der 
hierbei benötigten Tafeln für vollständige elliptische 
Integrale auch jetzt kein Problem ist — z. B. „Jahnke- 
Emde‘“ — liegt für den Praktiker in der Regel kein 
Bedürfnis vor, aufihren Gebrauch zu verzichten. 


Der von Herrn Potthoff in der Einleitung 
seiner zitierten Arbeit vertetene Standpunkt, eine 
Formelserunanschaulich, wenn man die Dar- 
stellung mif Hilfe elliptischer Integrale wähle, sie 
werde anschaulich, wenn man den Gauß- 
schen Algorithmus des arithmetisch-geometrischen 
Mittels einführe, dürfte Ansichtssache sein. Ähnlich 
subjektiv verhält es sich wohl mit dem Begriff „ele- 
gant‘‘, wobei man unterscheiden kann zwischen der 
Eleganz der äußeren Form — z. B. Formel (3) in der 
vorstehenden Zuschrift — und der Eleganz in der 
Handhabungsmöglichkeit einer Gleichung. Der Unter- 
schied, daß das vollständige elliptische Integral nur 
von einer Größe, dem Modul, abhängt, das arith- 
metisch-geometrische Mittel zweier Größen jedoch — 
und zwar symmetrisch — von diesen beiden, führt 
dazu, daß einigermaßen vollständige Tabellen dieser 
Funktion unhandlicher sein müßten. Sowohl die 
Tabelle von Gauß (Werke III, 403) als auch von 
K. Hayashi (Tafeln der Besselschen, Theta-, 
Kugel- und anderer Funktionen, Berlin 1930, 62—69, 
vgl. auch S. 121) beziehen sich auf den Sonderfall des 
arithmetisch-geometrischen Mittels zwischen der 
festen Zahl1 und einer zweiten Zahl <1. Daher 
wird. dort das arithmetisch-geometrische Mittel als 
Funktion einer Variablen behandelt und somit auf 
die Darstellung der — „eleganten‘‘ — Symmetrie in 
den beiden Argumenten verzichtet. Dann ist es aber, 
wie auch Formel (2) bei Herrn Potthotf zeigt, bisauf . 


7 
den Faktor % nichts anderes als der reziproke 


Wert des komplementären vollständigen elliptischen 
Integrals erster Gattung. E 
Anders liegen die Verhältnisse bei den arithmetisch- 
geometrischen Mitteln von mehr als zwei Größen, 
für die man keine Darstellung durch andere bekannte 
Funktionen besitzt, : 


Berlin. 0. Emersleben. 
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